РАСЧЕТ ПОТЕНЦИАЛЬНЫХ 
И ПРОГРАММНЫХ УПРАВЛЕНИЙ 
ДЛЯ ЛИНЕАРИЗУЕМЫХ СИСТЕМ 



























ГЬсударстьѳняый комитет РСФСР по делам науки и высшей школы 
Уральский ордена Трудового Кра'оного Знамени политехнически! 
институт им.С.М.Кнрова 


РАСЧЕТ ІЮТЕНЦИАДЬНЫХ И ПРОГРАММНЫХ 
УПРАВЛЕНИЙ ДЛЯ ЛИНЕАРИЗУЕМЫХ СИСТЕМ 


Межвузовский сборник научных трудов 


Екатеринбург 

1991 


УДК 62.50 


Предложены эффективный в инженерных приложениях метод рас- 
чета потенциальных управлений с обратной связью и эффективная 
численная реализация методики расчета на основе метода конеч- 
ных элементов. 

Основное содержание составляют статьи, посвященные разра- 
ботке методов расчета управлений с обратной связью для нелиней- 
ных систем и программных управлений. 

Программные управления рассчитываются для линейных и не- 
линейных систем с помощью методов типа принципа максимума. Со- 
ответствующие постановки задач состоят в управлении объектами 
относительно заданной совокупности целевых множеств в фазовом 
пространстве . 

Ряд статей сборника содержит результаты по методам расче- 
та импульсных управлений, а также инженерные рекомендации по 
раочѳту управлений для конкретных объектов. 

Сборник предназначен специалистам в области автоматиза- 
ции управления и студентам старших курсов специальности 2201 
"Автоматизация и управление в технических системах". 

Рецензент кафедра прикладной математики Ураль- 
ского государственного университета им. А.М. Горького 

Утвержден рѳдакциовно-мадатольским советом института 

Редакционная коллегия: 

от. науч. оотр. ,канд. тѳхн.наук А.Н.Оботшш - отв. редактор; проф. , 
д-р тѳхн.наук Ю. И. Алимов (Урал. политехи. ин-т) ; д-р физ.-мат. 
наук С.Т.Завалишин (У$) АН СССР); проф. , д-р тѳхн.наук В.Г.Ла- 
бунец; инж. В.П.Серов ~ отв. секретарь (Урал. политехи. ин-т) ; 
д-р физ-мат. наук А. Г. Ченцов (УЦО АН СССР) 


© 


Уральокий политехнический 
институт им. С. М. Кирова, І99І 


ПРЕДИСЛОВИЕ 


Разработка инженерных методик расчета законов управления 
с обратной связью, обеспечивающих требуемые динамические свой- 
ства, являетоя одной из основных проблем при практическом при- 
менении методов математической теории управления. 

Требования, предъявляемые к системам на практике, являют- 
ся противоречивыми. В математических постановках чаще всего 
формулируются задачи с единственными Цилями синтеза, например 
задачи оптимального управления. В инженерных приложениях обыч- 
но нельзя достаточно обоснованно указать функционал качества 
системы, поскольку при этом нужно уметь оценивать соответствие 
свойств системы и выбранного функционала. Последняя задача на- 
много труднее по сравнению о рѳечетом оптимальных управлений. 

В сборнике представлен ряд статей, направленных на созда- 
ние эффективной методика инженерного синтеза законов управле- 
ния с обратной связью (динамических регуляторов для нелинейных 
систем). Соответствующие управления названы потенциальными. 
Суть подхода заключается в том, что решение задачи синтеза не- 
линейной системы складывается аз решений множества задач ли- 
нейного синтеза. Заметим, что цель синтеза нелинейной оиотѳмы 
формулируется на основе частных целей синтеза к свойствам ли- 
неаризованных замкнутых систем. Такой подход пбзволяет разре- 
шить конфликт противоречий к свойствам системы на этапе линей- 
ного синтеза. Численно потенциальные управления рассчитывают- 
ся с помощью метода конечных злемѳнтов. 

Другая группа статей посвящена математическим вопросам 
расчета оптимальных угтавлѳний, о утверждениями, аналогичными 
принципу максимума. 

На примере задачи управления робо том-манипулятором в сре- 
де показано, что соответствующие управления содержат импульс- 
ные составляющие. Представлены результаты по методам раочета 
импульсных управлений. 

Ряд статей сборника содержит решение прикладных задач уп- 
равления конкретными объектами. 
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РАСЧЕТ ЗАКОНОВ УПРАВЛЕНИЯ ДЛЯ МНЕАГИЗУЕМЫХ 
СИСТЕМ 

Введение 

Основой построения моделей реальных физических объектов 
о помощью дифференциальных уравнений является идея линѳэризв- 
щи. Еоли правые часта дифференциальных уравнений - непрерыв- 
но дифференцируемые функции, то это означает, что "в малом" 
зти функции линейны, т.ѳ. в окрестности каждой из точек опре- 
деления они о точностью до бесконечно малых второго порядка 
описываются линейными функциями. Фактически все основные физи- 
ческие законы в дифференциальной форме записи получены на ос- 
нове перехода от локальных описаний, справедливых для малых 
приращений, к глобальнда описаниям с использованием диффѳрен*- 
іральных уравнений. Однако глобальность описания объекта диф- 
ференциальным уравнением является условной, поскольку практи- 
чески используемые численные методы решения дифференциальных 
уравнений основаны на "сшивании" локальных решений. 

Для многих задач инженерного характера представляется 
естественным задавать желаемые глобальные свойства замкнутым 
оиотемам управления путем ваданкя их локальных свойств. Ко- 
нечно, соответствие локальных и глобальных свойств нелинейных 
систем не является в общем случеѳ простым. Для часто встречаю- 
щихся на практике задач управления, в которых требуется полу- 
чить одинаковые или слабо меняющиеся свойства замкнутых сис- 
тем во всей допустимой области рабочих режимов, обеспечение 
локальных свойотв приводит и к аналогичным глобальным свойст- 
вам. 

Нелинейный объект можно представить в виде множества ли- 
нейных объектов путем линеаризации в каждой точке допустимой 
области. Каждый линеаризованный объект описывает поведение ис- 
ходного нелинейного объекта в некоторой, возможно очень малой, 
окрестности точки линеаризации. "Сшивание" линеаризованных 
объектов определяет исходный нелинейный объект. 
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Надлежащим расчетом законов управления можно получить 
замкнутые линеаризованные системнее необходимыми локальными 
свойствами , а закон управления для нелинейного объекта в целом 
# получить за счет "сшивания" соответствующих законов управления 
для линеаризованных систем. Именно эта идея и развивается в 
данной работе до получения практической инженерной методики 
расчета законов управления. Всю методику расчета невозможно 
изложить в рамках одной статьи, поэтому оснозное внимание бу- 
дет уделено математической постановке задачи расчета потенци- 
ального управления, ее формальному решению и некоторым свой- 
ствам замкнутых систем вблизи точек равновесия. 

Использование термина "потенциальный" подчеркивает, что 
нелинейный закон управления в виде обратной связи получен в ре- 
зультате выделения потэнциального поля из векторного поля, об- 
разованного всевозможными векторами коэффициентов обратной свя- 
зи для линеаризованных систем. 

Постановка задачи 

Пусть объект управления описывается нелинейным дифферен- 
циальным уравнением вида 

х(і)~!(х(і),и(1)), х(0) = х°, (I) 

где х - П -мерный вектор состояния; и - т -мерный вектор 
управления; х° ~ начальное состояние, функция у определе- 
на в некоторой области С х V из произведения пространств 
Я х Я и имеет в этой области непрерывные частные производ- 
ные д//дх и д]Іди , не зависящие от управления и . Символы 
д{!дх л дЦди обозначают матрицы частных производных вектор- 
ной функции у по переменным состояния X и управления 
и . Принятое допущение, что частные производные д^/дх и 
д}Іди не зависят от и , означает независимость свойств 
линеаризованных систем от управления- в данный момент времени. 

Для подавляющего большинства объектов задачи синтеза нели- 
нейного управления в виде обратной связи могут быть поставлены 
таким образом, чтобы удовлетворялись сформулированные выше ог- 
раничения на функцию у . Множество С в пространстве- со- 
стояний определяет область допустимых режимов работы объекта, 
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С математической точки зрения полагаем, что это множество 
является замкнутым и имеет Липшицеву границу Г. 

Пусть и * - некоторое фиксированное значение управле- 

ния их* - некоторое фиксированное значение вектора состоя- 
ния, тогда линѳарэозаннѳя система имеет вид • 

6±(і)~ & (х*)бх(і) + ^ (х*)6и(і)+}(х', и *), ( 2 ) 

где бх(і)- х(і)-х*, ди(і) ~и(()'-и ч , бх(О)- х°-х*° . 

Дифференциальное уравнение (2) описывает свойства нелинейной 
системы (I) в некоторой окрестности относительно состояния 
сг* и управления и* . 

Допустим, что линейные законы управления в виде обратной 
связи для линеаризованных систем, обеспечивающие необходимые 
им динамические свойства, известны и имеют вид [і, 2 ] 

ди(1)~-і(х*')дх(+) , ( 3 ) 

где і - некоторая матрица обратной связи. 

Не нарушая общности рассуждений, полагаем, что управление 
и является скалярным и і представляет собой вектор-стро- 
ку. Случай векторного управления следует непосредственно из 
рассмотрения скалярного случая. 

Поставим задачу: рассчитать закон управления и(Х) в ви- 
де обратной связи для нелинейной системы (I), обеспечивающий 
заданные овойотвв линеаризованных замкнутых систем вблизи каж- 
дого состояния X * . 

Нетрудно показать, что в общем олучвѳ такого закона управ- 
ления не существует. Линеаризованная система (I) с законом уп- 
равления и(х) имеет вид 

6± ( і ) = '№(» + $ (*' (х ЧМ > ( 4 ) 

где и* = и(х*). 

Для того, чтооы линеаризованная система (4) обладала заданны- 
ми свойствами, в соответствии с (3) необходимо выполнение ра- 
венства 
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^(х')~-Цх'). (о) 

Есліі потребовать, чтобы свойства линеаризованных систем обес- 
печивались для кавдого состояния .г * , то (5) преобразуется 

к виду 

^~(х) = -Цх) , хе а, (6) 

т. ѳ. вектор-строка Ь(т) должна быть градиентом некоторой 
скалярной функции. 


Необходимым и достаточным условием существования функции 
и , удовлетворяющей (6), является равенство частных произ- 
водных [3] 

г! г. г)Р- 

^ • • • » Л , (?) 


Ці 

дх- 

/ 


де. 

дХ: 


Сп(х)). 


где принято , что і(х) = ( 6, (х), 

Поскольку элементы матрицы Л определяются только из усло- 
вия получения заданных динамических свойств линеаризованных 
систем, то условие (7) практически никогда выполняться точно 
не будет. 

Задачу синтеза нелинейного закона управления из множества 
законов управления сформулируем в следующем виде. 

Задача расчета потенциального управления . Определить за- 
кон управления в виде обратной связи, обеспечивающий минимум 
функционала 


эе (и) - + і(х)Ц*сіх , (8) 

где Й'І/ обозначает евклидову норму векторов и матриц в ко- 
нечномерных пространствах. 


Смысл функционала (8) на языка теории поля ооотоит в том, 
что из векторного поля, порождаемого множеством векторов обрат- 
ной связи і , выделяется потенциальное поле. 

Введение функционала де (и) вида (8) не является един- 
ственно приемлемым. В зависимости от конкретной физической за- 
дачи в функционал (8) можно, например, ввести весовую функцию, 
позволяющую обеспечить более точное приближение желаемых свойств 
замкнутой системы в заданном подмножестве допустимой области 
С . 


- 7 - 


Важно отметить, что в предлагаемом подходе вид функцио- 
нала (Ь) должен слабо влиять на свойства замкнутой системы, 
которые в основном должны определяться векторной функцией 
Цх ) . Если это не выполняется, то следует изменить желаемые 
свойства линеаризованных систем, приводящие к изменению век- 
торной функции і(х) . 

функционал типа (8) допускает наиболее простое решение 
задачи на минимум, что является очень важнш для практическо- 
го применения предлагаемого метода. 

Решение задачи расчета потенциального управления 

В предыдущем разделе задача раочѳта потенциального управ- 
ления поставлена из обычных физических соображений. Дня кор- 
ректной математической постановки задачи необходимо рассмот- 
реть класо функций, на котором будем искать минимум функцио- 
нала (8). В качестве пространства функций выберем гильбертово 
пространство Ѵ/^ П (С) функций и , имеющих обобщенные 
производные ди/ дх ив пространства і г (С) интегрируе- 
мых о квадратом функций [ 4 ] . 

Конструкция пространства VI ^ (С) достаточно сложная, 
но совершенно необходима для решения вопроса о существовании 
решения задачи раочѳта потенциального управления. Отметим ос- 
новные элементы структуры пространства ѴІ г и) (С) [4] . 

Как известно, пространство і г (С) представляет ообой 
пространство измеримых функций, суммируемых с квадратом, т.е. 
для каждой функции и ив величина [и*(х)Ых 

существует и конечна. Интеграл понимается в смысле Лебега. Ска- 
лярное произведение функций из і г ( С) определяется соотноше- 

нием 

( и і&\ г (а) “ !и(х)#(х)ах . (9) 

& 

Для некоторых функций Ъз (С) можно ввести еще одно ска- 
лярное произведение. При этом нам понадобится понятие обобщен- 
ной производной. 

Отображение ди/дх является обобщенной производной 
первого порядка функции и е , если для каждого эле- 

мента Фе С™ (С) выполняется соотношение 
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(ІО) 


/ л ? Ф 0 (' г ^' 2 = - ] и № 4 х <Х ^ Х * 

а а 

где С “(С) обозначает пространство основных бесконечно диф- 
ференцируемых функций с компактным носителем. Естественно, что 
если функция И дифференцируема в обычном смысле, то 
обобщенная производная совпадает, о обычной производной. 

Используя понятие обобщенной производной, введем для под- 
множества функций из 1*2 (С) , имеющих обобщенные производ- 

ные первого порядка, скалярное произведение 


(Ч^кѵ/Ѵ) * ( и ^\ г (с) + {дх ' с/х)** 


(4) ' 


(II) 


где для векторных функций (зектор-отрок) принято 

' ди\ д&‘ , 


(ди_ д$\ = /7 

\дх 1 с/х/^Ф •/ { 


дх) дх 


с/х 


( 12 ) 


Штрих обозначает операцию транспонирования матриц. 
Подпространство функций из і г (и) со скалярным произве- 
дением (II) образует пространство Ы г <г> (С) , воли граница 

Г области С Дипшицѳва [4] . 

Липшицѳва граница представляет собой іипѳрпов&рхность 
П -мерного пространства , в окрестности кавдой точки которой 
можно ввести локальные координаты и описать участки гиперпо- 
верхностей в этих локальных координатах Липшицев ыми функциями. 
Области, имеющие гладкие или кусочно гладкие границы без осо- 
бых точек .заведомо являются Лиігаицѳвши. 

В соответствии с определенна! выше скалярным произвѳдѳни- 
(II) метрика в ( С) определяется равенством 


ем 


И и И уМф “ ( и » и) ѵу/б (е) 


(13) 


Пространство содержит плотное подпространство С°°(С) 
непрерывных и дифференцируемых функций и поэтому каждый элемент 
из Ѵѵ^ ( 5) может быть с произвольной степенью точности в мет- 
рике (13) аппроксимирован бесконечно дифференцируемыми функция- 


ми. 
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Рассмотрим задачу отыскания элемента и °е , до- 

ставляпцѳго минимум функционалу (8). Для удобства дальнейшего 
наложения представим функционал (8) в виде 

Щи) 2 1& (*)• С ’Шх +//і(х)/‘4х . (14) 

Используя формулу интегрирования по частям, получим 


л+^ифІ(4М$+І/и(х)/ л с/х,(І5) 


где V - внешняя нормаль к Г. 

В [4] показано, что минимум функционала (15) на простран- 
стве VI® (С) достигается на некотором элементе и" . Любые 
два минимизирующих управления и 9 отличаются на постоянную 
почти всюду на С . Функции и° являются слабыми решения- 
ми уравнения Пуассона 


і,і от} Т?і дх, ’ 


(16) 


с краевым условием на границе 

Ц- ѵ« - и , х « г. (і7) 

Слабое решение существует тогда и только тогда, когда 

(18) 




Управление и 9 может быть приближенно получено с любой точ- 
ностью по методу Ритцв и, следовательно, о помощью метода ко- 
нечных элементов (МКЭ). Это открывает возможность эффективно- 
го численного решения задачи раочѳта потенциального управления 
о помощью МКЭ. 


Ѵ-зщощ дасйсчаа чийтом. д дахаддиашм. люавдвнш 


Как уже отмечалось выше, потенциальный закон управления 
и(х) решает вадачу "сшивания" линейных законов управления, 
рассчитанных для каждой точки линеаризации из области С 
Таким образом, необходимые динамические свойства нелинейных 
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систем обеспечиваются на этане линейного синтеза. Наиболее 
предпочтительной для практических вычислений является методи- 
ка линейного оинтѳза, основанная на задании желаемых собствен- 
ных значений и векторов матриц динамики [і , 2] . 

Поскольку не существует ваконов управления в виде обрат- 
ной связи, обеспечивающих произвольные заданные овойства лине- 
аризованным системам, то потенціальный аакон управления обес- 
печивает приближенное выполнение заданных овойств в каждой точ- 
ке линеаризации замкнутой системы из области О . Точные ко- 
личественные оценки степени приближения завиоят от сиорооти 
изменения овойств иоходннх линеаризованных систем и от качест- 
ва рассчитанных для них линейных законов управления. 

В качестве примера рассмотрим устойчивость замкнутой оис- 
тѳыы с потенциальным управлением. Замкнутая система о потенци- 
альным управлением и (X) описывается уравнением 

±~1(Х,и(&). (19) 

Функцию Ляпунова для системы (19) запишем в виде 

іѴх) » !'(х, и(х)) Ѵ(х) ! ( х , и(&) , (20) 

где Ѵ(х) - некоторая оишѳтрячная матрица. 

Полная производная функции & ■ в силу системы (19) имеет 

вид 

‘(х, и(Х))(А'(х) Ѵгх) * \/(х) А(Х)+ $)/(х,и ($) , (21) 

где 

мо. Щ“? ± ♦ Щит. , (гг) 

дх ди ох 

(23) 

Матрица дѴ(х)/дх является прямоугольной и- имеет размер- 
ность пх п* и получается заменой каждого элемента Ул* мат- 
рицы V на отроку дѴ^ /дх . Знак " 9 " обозначает тензор- 
ное произведение матриц, /„ - единичная матрица размернос- 

ти п*п . 

Матріцу У(Х) выберем с помощью матричного уравнения 
Ляпунова 
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А'(х)Ѵ(х)+ Ѵ(х) А (х)--4, (24) 

ідѳ В - некоторая положительно определенная матрица. 

йоли потенциальный закон управления обеспечивает прием- 
лемые динамичѳокиѳ свойства замкнутым линеаризованным систе- 
мам (т.ѳ. заданные расположения в левой полуплоскости собствен- 
ных значений и необходимые направления собственных векторов 
матриц динамики А(х) для каждого х из В ) , то У(х) 
будет положительно определенной матрицей Ухе В 

В случае, когда нелинейная система (19) имеет одну точку 
равновеоия в области В , траектории оиотѳмы будут стремить- 
ся к этой точке равновесия, если матрица 

М(х)-АІ'х)ѴГх)+Ѵ/х)А(2)+-2^ ) (і п в{(х, и(х)))--<3+ а (25) 

*(!„•/ (х,и(*і) 

является отрицательно определенной в В . Производное ото- 
бражение дѴ(х)!дх можно определить путем дифференцирования 
(24). В результате получим систему матричных уравнений 

А(Х) +А'& 2Ш « - (щ Ѵ(х)) л. (26) 

' I 

Анализ уравнений (26) показывает, что для получения уо- 
тойчивой системы с потенциальнмі управлением желательно обѳо- 
печить медленное изменение овойотв замкнутой линеаризованной 
сиотѳмы в пределах области 0 . 

Заключение 

В работе предложен метод раочѳта нелинейных систем путем 
"ошивания" законов управления для линеаризованных систем. Чис- 
ленное отыскание нелинейных потенциальных эаконов управления 
осуществляется с помощью метода конечных элементов. Пример 
расчета потенциального управления методом конечных элементов 
представлен в [5] . Предложенный метод без существенных изме- 
нений пригоден для раочѳта нестационарных и дискретных систем 
управления. 


дт 

ъ: 
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УДК 62.50 А.Н.0ботнин, Л.П.Алеоѳнко 

(Уральский политехничес- 
кий институт) 

РАСЧЕТ ПОТЕНЦИАЛЫЮГО УПРАВЛЕНИЯ МЕТОДОМ 
КОНЕЧНЫХ ЭЛЕМЕНТОВ 

Введете 

Метод расчета нелинейных законов управления в виде обрат- 
ных овязей для класса нелинейных систем предложен в [і] . 

Этот метод позволяет получать потенциальные законы управления 
для нелинейных сиотем через законы управления для линеаризо- 
ванных систем. 

Пусть уравнения динамики нелинейной системы имеют вид 

х(0)-х°, іе[0,°°), (I) 

где хе Ос Я п , и е , функция I непрерывно дифферен- 
цируема и матрицы частных производных д^дх , д^/ди не за- 
висят от управ летя и(і) 

Полагаем, что для каждой точки линеаризации х* в мно- 
жестве & и для фиксированного и * известен закон уп- 
равления для линеаризованной системы 
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(2) 


6±(і)-& (х')6х + -& ди , 


имеющий вид 


6и(і)т-1(х*)6и(і), 


( 3 ) 


дде дх •* х-х щ , ди^и-и*, 1(х*) - вектор обратной овявя для 
линеаризованной системы. 

Потенциальный,- вакон управления минимизирует выпуклый 
функционал 



(4) 


на пространстве (І, 2] , и является олабж решени- 

ем уравнения в частных производных типа Пуассона [і] 



(б) 


} ( V - внешняя 


нормаль к дС ). 

Хотя потенциальное управление. и(х) является решением 
уравнения в чаотных производных типа Пуассона (5), использо- 
вать его для численного нахождения решения рассматриваемой за- 
дачи управления (I )— <4) неэффективно. Наиболее предпочта тель- 
ным является непосредственное нахождение минимума функционала 
(4) с помощью метода конечных элементов (МКЭ). Суть метода 
заключается в разбиении области С на конечные элементы, 
т. ѳ. на множества О с , С *■ и &■ , в пределах которых задают- 
ся интерполяционные полиномы' и‘(Т) , определяемые узловы- 

ми значениями, с последующей минимизацией функционала по уз- 
ловым значениям. 


Численный алгоритм расч ета управления 


Согласно требованиям к МКЭ из [з] , представим замкнутую 
область 6 с непрерывной по Липшицу границей в виде объе- 
динения замкнутых с непустой внутренностью нѳпересѳкающихся 
множеств 6 І , т. ѳ. 
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( 6 ) 


где к - количество конечных элементов; 0 ; - некоторые 

замкнутые области, называемые в дальнейшем конечными элемен- 
тами. Пусть в пределах каждого конечного элемента О- с /Р" 
искомая функщя и представлена интерполяционным много- 

членом и 1 (х) , і = 1, к . Каждый интерполяционный полином 
является функцией у еловых значений: 

(и (х, и(х р ))\ (7) 

где х п ...,Хр -узловые значения; и(х, и(Хр) - ис- 
комые узловые значения; р - количество узловых значений. 

Интерполяционный полином в пределах элемента &■ запи- 
сываем в виде 

и 1 (х) -А/;ГХ)Ц’, Хе &; , і- V* , (8) 

где Л/; (X) » (М (і (X),..., /Ѵ. р (х)) - функция формы, определяемая 

видом интерполяционного полинома. Любая функция формы облада- 
ет свойством 

- 

где бу - вектор с единичным у -м элементом и нулевыми 
остальными элементами. 

В результате разбиения области (г на конечные эле- 
менты и введения интерполяционных многочленов функционал (4) 
примет вид 

эе(и)- ^ (іо) 

где 

(и 1 ) - ]/ - у і'Шх +!! ЦХ)І*СІХ . (II ) 

С учетом (8) функционал зе (и) будет функцией увловых значе- 
ний 

эе(и) -с эе(Ѵ)~ Е я. («0 , ( 12 ) 

«»/ * 

где г 

ДО- / ЦіГ'^ЦсІТ+2 ’У^ - + іЦих)Ц‘сіх . 

Ъ 1 а і I 


0 } Х 4 & * 


(9) 


( 13 ) 


Дифференцируя (12) но вектору 
мое и достаточное условие минимума 

* <• ! Яч\! Яч\ ' к 


II получаем неоЗходи- 


[.5І(Щ'Н "--ь !(■§)"*»*■ (І4) 

1 в, о.' 

Соотношение (14) лредставляет собой систему линейных уравне- 
ний относительно вектора узловых значений V 

Уравнение (14) всегда разрешимо относительно II и име- 
ет бесчисленное множество решений вида 


и = и + ѵ 


(!)■ 


(15) 


где V - одно из решений; ѵ - произвольное число. 

Минимум функционала (4) определяется с точностью до по- 
стоянной функции. Еоли в (14) задать одно из узловых значений, 
или к (14) добавить уравнение, определяющее функцию и 0 в 
заданной точке, то получившаяся система уравнений будет иметь 
единственное решение. 

Если конечные элементы (}. представляют собой симплек- 
сы, то согласно [4] можно использовать линейные интерполяцион- 
ные полиномы и‘ (X), .В этом случае отображения 

дМ- /дт не зависят от г и система (14) примет вид 

(ЫШ§1) и --ШѴ* • < к > 

где г г ^ = I с/х - объем области С { ; (17) 

Сі ш і іфйх . ( 18 ) 

' 6; 

Рассмотрим последний вариант более подробно. Поскольку 
0. является п -симплексом, то он является выпуклой обо- 
лочкой п+1 точек из пространства К/ 1 . Обозначим эти 

точки соответственно у/ , у ‘ , ... , у я ‘,, , тогда 

С,- [х/хш Т1у і Л. і \ П А;*', У-І,** 1 ]- ( І9 > 

3 соответствии с принятыми ранее обозначениями можно также 
полагать, что узловые значения х, , ..., Хр являются одновре- 
менно вершинами симплексов. 
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Линейный интерполяционный полином и '(т) напишем в век- 
торно-матричном виде 

и‘(х) - М 1 (х) , і - йн , (20) 

где М‘ - вектор-строка размерности п + 1 . Для записи функ- 
ции формы /V; (х) необходимо представить М‘ через уз- 
ловые значения управления. Нетрудно видеть, что 

і -/ 


\-\' и і У 

1 і Ул*і/ 


(21) 


( 22 ) 


(20) можно представить в другой форме: 

“‘«'-(“м 

\< \и(УІ,)) 

Таким образом, функция формы /Ѵ< (х) , выраженная через те 
"реальные" узловые значения х, , .. . , X, на п -симпѳк- 
се (?■ , имеют вид 


Ъ(х)-(1х) 




(23) 


где соответствующая С,- матрица Е- состоит из нулей и 
единиц; ее строки являются относительными координатами 
и(У})* ••• і и (у 1 п+}) 3 столбце V из уравнения (16): 


( и( *'Ч ) , гг 

\“(УІ+,)І ‘ 


Производная функции формы: 


№ г 

дх “ •' 


! ! і 


\~і 


'(оо„ ._0_ у 


і 

где І п - единичная матрица размерности п 

Объем п -мерного симплекса &■ с вершинами 
у} , . . • , , определяется формулой : 


(24) 


( 25 ) 
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( 26 ) 


Ѵ.-ІсІеіЬЧ - т . 

гдѳ Л‘ - матрица преобразования 

(27) 

При вычислении выражения (18) олѳдует учесть, что вектор- 
отроки обратной связи раоочитываются с помощью целого 

комплекса программ для решения задач линейного синтеза, поэ- 
тому точное вычисление интеграла (18) является времяѳмкой про- 
цедурой. С другой -тороны, если разбиение области О на 
элементы достаточно мелкое, то і(х) будет слабо 

измѳнятьоя в пределах &, , поэтому интеграл (18) можно вы- 

числять о помощью теоремы о орѳднѳм, фиксируя какое-либо зна- 
чение X* иа области (?. и полагая, что 

і в . "1*(Хщ) . (28) 

В качестве Х к щ рационально выбирать точку, близкую к "цент- 
ру" симплекса, например в центре тяжести симплекса. 

Более точные результаты можно получить, если ввести интер- 
поляционные полиномы для 1(Х) , иопольауя в качестве узло- 

вых вначений вершины симплекса. Если принять линейный закон 
интерполяции, то 

1(х)*$1 /.(и*)*- , х е 0. , (29) 

'Ч у 

где иі- - относительные координаты точки х в оимплѳкоѳ 
С , определяемые оогдаоно (19). После несложных вычислений 
получим 

ізо) 

Алгоритм раочѳта потенциального уравнения получается на 
основе соотношений (16), (25), (26), (27) и (30). Как уже от- 
мечалось выше, для однозначной разрешимооти оиотемы (16) необ- 
ходимо задаваться величиной управления и в произвольной 
точке х области 0- . Величину и удобно задать в точ- 

ке равновесия нелинейной сиотѳмы так, чтобы выполнялось усло- 
вие {(х * и(х*^)-0. 
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Для уменьшения объема вычислительной работы при решении 
системы линейных уравнений (16) необходимо специальным образом 
вводить нумерацию узловых значений (координат вершин симплек- 
са). При правильной нумерации вершин симплекса получится лен- 
точная структура матрицы системы линейных уравнений. Посколь- 
ку нумерация вершин симплексов имеет большое значение, рас- 
смотрение этого вс..роса выделим в следующий раздел. 

Разбиение многомерной области на конечные элементы 

В данном разделе рассматриваются области в виде многомер- 
ного параллелепипеда. Более сложные облаоти могут быть рассмот- 
рены с использованием комформного отображения на соответствую- 
щий многомерный параллелепипед. В аадачах управления обычно 
нет необходимости рассматривать сложные области, поскольку кон- 
фигурации областей определяются допустимыми режимами работы 
объектов управления в пространстве состояний. 

В п -мерном пространстве рассмотрим облаоть 5 в виде 
прямоугольного параллелепипеда. Предположим также, что исполь- 
зуются равномерные сетки разбиения по каждой переменной 
X , , . . . , Х„ о количеством узлов ^ , . . . , ф п соответствен- 
но. Для определенности считаем, что переменные оостояния Х П ...,Х Я 
пронумерованы так, что числа /] п расположены в 

порядке возрастания. 

Разбиение параллелепипеда 0 на симплексы и нумерация 
узлов осуществляются в два этапа. На первом втапѳ & раз- 
бивается на элементарные параллелепипеды меньшего объема и спе- 
циальным образом нумеруются узловые значения для получения ми- 
нимальной ширины ленты матрицы в (16). На втором втапѳ элемен- 
тарные параллелепипеды разбиваются на оимплѳкоы. 

Вершины элементарных параллелепипедов нумѳруютоя по слоям. 
Каждый ( п - I )— мерный слой получается путем сечения п -мер- 
ного параллелепипеда гиперплоскостью, ортогональной координате 
Х п с самым большим количеством узловых точек . Коор- 
динаты вершин нумеруются последовательно от первого слоя до 
последнего о номером ц п . В свою очередь нумерация каждого 
( п - І)-мѳрного слоя осуществляется по ( п - 2)-мѳрным слоям. 

Для этого проводятся сечения ( П - I )-мѳрного олоя гипѳрплос- 
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костями, ортогональными координате X П1 в соответствии 
с принятым разбиением по этой координате. 

Пример разбиения 
трехмерного параллелепи- 



педа на элементарные и 
нумерации вершин пред- 
ставлен на рас Л. 

После нзадерации вер- 
шин элементарных п -мер- 
ных параллелепипедов 
осуществляется их разбие- 
ние на п -мерные симп- 
лексы. Процедура разбие- 
ния одинакова для всех 
элементарных параллеле- 
пипедов. 

Проиллюстрируем 
сначала эту процедуру на 

примере разбиения 3-мѳрно т, о куба, представленного на рио.2. 

Из геометрических соображений яоно, что каждый симплѳко, ле- 
жащий в кубе, о общей вершиной в точке 
О может быть получен как выпуклая ком- 
бинация 4 точек, одна из которых 0, а 
три другие лежат на грани куба, не со- 
держащей точку 0. Вопрос состоит в том, 
чтобы поотроить симплексы, заполняющие 
весь объем куба без общих внутренних 
точек. 

Для разрешения этого вопроса ввѳ- 


РиоЛ. Пример 3-мерного парал- 
лелепипеда 



Рис. 2. Разбиение 
куба на симплексы 


дѳм обозначения 


в е 

і * 


для единич- 


ных ортов. Можно показать, что искомая 
оиотема симплексов получается как вы- 
пуклые оболочки точек (О, е і ,е { + , е, + е м ] , где А- 

произвольный набор различных индексов из трех чисел I, 2, 3 
с всевозможными перестановками индексов. Общее количество 
симплексов для трехмерного случая 3! 

Рассмотрим разбиение на симплексы п -мерного куба. По- 
лагаем для простоты, что одна из его вершин совпадает с нача- 
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лом координат. Единичные орты обозначим через е 1 , . . . , е п . 
Каждый симплекс является выпуклой оболочкой точек 

{°> е е «, * е Л г *— * е л„ 1 ' 

где о*,,... , а. п - произвольный набор различных целых чисел 
от I до п . Общее число таких наборов равно (л!) . Таким 
образом, п -мерный куб разбиваѳтоя на п! п -мерных симплек- 
сов, полностью заполняющих его объем с общей вершиной. 

Раэбиение параллелепипедов на симплексы проводится анало- 
гично. Рассмотренный выше алгоритм позволяет построить на ЭВМ 
процедуру разбиения многомерной облаоти на симплексы и нумера- 
цию вершин с наименьшей шириной ленты матрицы в уравнении (16). 

Заключение 

Процедура расчета потенциального управления достаточно 
сложна и ориентирована на реализацию в виде программного обес- 
печения для ЭВМ. Практическое использование авторами предлагае- 
мой методики расчета законов управления для нелинейных систем 
показало ее высокую эффективность. В [I] достаточно подробно 
рассмотрен расчет потенциального управления для модельного 
примера. 
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УДК 620.975 


A. М.Тарасьѳв, А. А. Успенский, 

B. Н. Ушаков (МММ УрО АН СССР) 


АЛГОРИТМ ПРИБЛИЖЕННОГО ПОСТРОЕНИЯ МНОЖЕСТВА 
ПОЗИЦИОННОГО ПОГЛОЩЕНИЯ В ЛИНЕЙНОЙ ЗАДАЧЕ СБЛИЖЕНИЯ 

С ВЫПУКЛОЙ ЦЕЛЬЮ 

Рассматривается задача построения множества позиционного 
поглощения в линейной дифференциальной игре сближения с выпук- 
лой целью. Считается, что момент окончания игры фиксирован. 
Приводится алгоритм приближенного построения множества позици- 
онного поглощения в случае фазового пространства Р т . Рабо- 
та согласуется о исследованиями [і-9] . 

I. Пуоть задана конфликтно-управляемая сиотема, движение 
которой описывается уравнением 

^--А(і)х + В(і)и + С( Огг, (І.І) 

здесь Т-т -мерный фазовый вектор системы; и, V - векторы уп- 
равляющих воздействий І-го и 2-го игроков, удовлетворяющие ог- 
раничениям 

иеРсР* , ѴеОс/? 9 , (1.2) 

ідѳ Р и а - компакты, Р р и - евклидовы пространст- 
ва; А(і), В (0, С(і) - непрерывные на заданном промежутка вре- 
мени [і 0) Ѵ] (і д < Ф) матрицы-функции. 

В задаче о сближении, стоящей перед І-м игроком,' требует- 
ся обеспечить попадание движения X [і] оистѳмы (І.І) в 
момент тТ на компактное выпуклое целевое множество М , 
содержащееся в фазовом пространстве Р т . Эту задачу можно 
решать с помощью экстремальной конструкции [1-3] , основным 
элементом которой является построение множества позиционного 
поглощения. Согласно определению [і] , множество позиционного 
поглощения еоть множество всех повиций [с о , і?1*. Р т , 

из которых разрешима задача сближения. Зная матрицы А СО, &((), 
С(і)' и целевое множество М , мы можем легко найти ту ог- 
раниченную область Л с [і 0 , іФ] х Р т , в которой со- 

держится множество позиционного поглощения. Область Л на- 
зовем пространством позиций ? «. (і г) игры. 
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Введем вспомогательные понятия и обозначения: 

І(і,х,и, ѵ) = А(і)х + В(і)и + С(і)ѵ ; 

у (і , х, С) “ /тюх тіп е'{( і, х , и, ѵ) ; 
иеР ѵеа 

Г(і,х) ш со[{ : 1 = І(і,х,и,ѵ), иеР, ѵер] ; 

Г' (*, Л - Г/ * *7>, *), Ц « Я^, х.П] ; 

5 =• [( : С е Р т } ЦЦ~1) , 

где символ Со[/] означает выпуклую оболочку множества 
[/) \ С - вектор из 5 ; Ц С// - евклидова норма векто- 
ра ( . 

Символом Х е (і*і , 1,) (і щ <і*) обозначим мно- 
жество всех тг-’ек х* е Р™ , в которые приходят в момент 
і* решения х[і] (І€ , ^ * ] ) дифференциального 

включения X е Г е (і, X) , х[І и ] - X* . 

Определение . Множеством позиционного поглощения Ѵ/° 
назовем максимальное из множеств IV в пространстве пози- 
одй В , удовлетворяющих условиям: 

1. К, сМ\ 

2. Для любых двух моментов , і*(І в « і щ < і * <і і?) 

любой точки х щ ( X* 6 ІѴ^ ) и любого вектора 

С ( I е 5) имеет место 

\Р^ПХ е (*'} , Т+)Ф Ф , 

здесь ѴѴ, = (х : X е Р т , (і х X) е V ] , 

Будем предполагать, что выполняется условие /Г : 

(і ш 9,г)еДк9, 

где 

Отметим, что в [2] вводятся множества 
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О Ц6Н- 


г г '(і, х)-г ѵ І,Х)П(І ;/е е'і 4 ?а, х,01, 


гдѳ Г (і , х) - некоторый шар в пространстве К 
тром в начале координат и достаточно большим радиусом. 
. Если обозначить при этом 

9М^,С) -тіа ?'} р'(*,х,1)= тох 

}*?*(*>*) 1 Г ’ /бГѴ*,*) ' 


( 1 . 3 ) 


и эаменить фигурирующие в условии Е функции р*(і, X , О, 
ф*(і,х,С) функциями (1.3), то условие Е будет заве- 
домо выполняться. 

Пусть выбрано разбиение Гш і , .. . , = і?} . 

Полагаем й ( - - - Е, . Определим ., числа е 0 - О , 

е 4 . ~со(а іш ,)+ (1 +-Л&і-і)Е і., , 1-7,... , /V , где 

- некоторая функция, монотонно убывающая, такая, что 


Л'/п 

4*0 




О ; ѵ? - /пах І/АСОІІ 


Для выбранного разбиения Г определим последователь- 
ность множеств IV*. , /• /V , И~ О . Полагаем 

у. « ~ [иг : и/е # т , тіп ІІѵ/*-гігЦ ** е.„\. 

с л/ 1 ѵ*тм • •’ 

Пусть множество V,. уже построено. Полагаем 

- {.га 6 Й т !щ іч пЪ і (і 1 „ ;і п ѵ) + Ф 

при любом Ее 5] , где (7*; ~ ы, + (7 *- Г е (і щ , ^ 

Определенной здесь системе множеств ѴѴ^ с Я т ,1=0, 
I,..., N можно поставить в соответствие множество 

Ѵ/(Г)= [((, ѵ) : ыв <і 4 Е. Ж] . 

Известно , что (ѵ/(г) , ѴѴ °) — 0 при сІіат(г)-~0 
[8] , где с/іа/п (V) - тот |^ +( -1,1. 

Символ рПѴ(Г),и' в ) означает хаусдорфово 

расстояние [іо] между множествами IV ^Г) и IV 0 

Цель настоящей статьи - описание алгоритма построения 
системы множеств ІѴ^ , і - О, 1 , . . . , Л' , которую 

назовем аппроксимирующей системой множеств, отвечающей разби- 
ению Г 
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Достроекиэ аппроксимирующей системы множеств { ѴѴ 4 . : 
і « А/ , Л- О] осуществляется согласно определению с 

помощью попятной процедуры от множества IV ^ , рассматривае- 

мого в конечный момент времени I? к множеству IV* , рас- 
сматриваемому в начальный момент времени і 0 . На каждом эта- 
пѳ [і ( , Ь 1 + 1 ] конструируем множество IV*. , исходя 

из построенного на предыдущем этапе [ і и , , ^ і+і ] множества 
IV*. . При этом на исходном этапе попятной процедуры пола- 
гаем . ' ' 

Для построения множества ІѴ*^ (і и Ж) обозначим 


<*Пф - Ч.„ пі(( ы : *і . *ЛЛ- 


Здесь 

Ь . *&7) - х [*і]+ & і г (Ь » *Л7) » 4 і - А 


Теорема I . 

К, - (х[іі]*# т : таг е й . А» Л А7 » О * . 


пип ('ѵг,есп и Л,. (*&:])+ Фі 

ѵ*м ім (*н# 

* ■» , если (х[і { ])-&. 

В силу линейности системы (І.І) и выпуклости цели Л? , 
множества ѴѴ<. , і - О,/, . . . , N выпуклы, и, следовательно, 
выпуклы множества М і . (х [і і }) при любых 

Х[іі ]е* т и /- О,..., N- 1 . 

Предполагаем, что наряду с условием Е * выполняется 
условие Н : множество IV* при любом і - О, . • N -1 
имеет непустую внутренность ШІ . 

Піи приведенных выше предположениях Е и И выполня- 
ется следующее утверждение. 
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Хад дама 2 . тах е а (I; , х[1 ( ], /) 

Ія $ 1 


- тах 


(^і і 


х [*і]> 1)-0, 


где 4 У*/*.;) - вамыканиѳ множества векторов І из 5 

внутренних нормалей к (х[Ь с ]) в точках множества 

л, №7) - М‘. иі (ЩТ)Пьл Х«и , ; *, , * А 7) • 

Здесь оимвол <?Ф оаначаѳт границу множества 0 
Заметим, что соотношение /7»7х г, (і , хГ{. 7, ^/7. 

ыі 0 (хЦ{і) • * ь * 7 

являщѳеоя критерием принадлежности точки х[і(] границе мно- 
жества к/*, , недостаточно удойно для конкретной реализации 

попятной процедуры на ЭВМ вследствие того, что максимум в ле- 
вой части этого соотношения вычисляется на множестве, завися- 
щем от искомой точки х [іі] , поэтому сделаем следующий 

яѳрѳход. 

Пусть • Поставим в соответствие точ- 

кэ х [<і*,] множество М* * ( х [*і+і]) ш Л * Г Д Ѳ 

0 г -[ Цх-х[і^,]Ц іг] (. г>0 ). 

При выборе числа Г полагаем: 

& і х [*і], і) - - * і 7 ~ ?)+?»<■ г*р иі] ) (Я ; 

4 */♦/ 4 * 

Аг(*ЬЛ-М°гПОМ^х[і иі ]\ і г (*[*іч]) - 

- с( (1:?е$, 1‘(х-у)<.0 для всех 


где ушА,(х[і и ,])]- 

Здесь оимвол означает замыкание множества в Л* 

Теорема 3 . Пуоть точка X ] связана о точкой 
х[і ІФІ ] е соотношением Мі \(х[і { ] )сіпіО г . 

Тогда олѳдупцие три утверждения эквивалентны: 

I. х[і { ]& дУ іх ; 
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г -і&*ьп е ‘‘ ( * ,,жМ,е) ’ 0 ' 

3 'е.™%ЬП 


(1.4) 


Соотношение (1.4) удобно тем, что максимум по 8 функ- 
ции с/ (і { , X [і- ] , 8) вычисляется на множестве 

і г (х [іц., ]) , не зависящем от искомой точки х[ (.] . 

Итак, в теореме 3 утверждается, в частности, что все гра- 
ничные точки х[{(] множества У ( . , и только они, 

удовлетворяют соотношению (1.4). Следовательно, перебирая точ- 
ки множества дУі +/ и решая уравне- 
ние (1.4) относительно х[і { ] , мы выделим границу дУ ( . 

множества , и, таким образом определим множество 

Приведем описание конкретной вычислительной процедуры, 
которая включает в себя и вычисление значений 

Предполагаем, что в задаче сближения множества Р , С/ , 

М - выпуклые многогранники с конечна* числом вершин. 

Ниже описывается попятная вычислительная процедура построе- 
ния системы множеств ( п*. : і ( е Г } , близкой к 

системе множеств ( У*. : 8, с Г } . Здесь близость множеств 
У*, и У і . ~ понимается в смысле хаусдорфовой метрики 
$ ( Ѵ/і , * ) ; Г - некоторое резбиѳниѳ отрезка 

[і в , ] , характер которого будет определен ниже. 

Отметим, что факт неточного вычисления множеств У}, , 

і і е Г связан с тем, что попятную процедуру мы будем вести 
не от множества М ем , а от множества М - выпуклого мно- 
гогранника. 

На начальном шаге процедуры полагаем «Л/ .В конеч- 
ном счете процедура сводится к построению для каждого момента 
і,- е Г границы дУі. множества на основе 

множества 

Р т . 

(I — 0 % А/) имеют 

которое обозначим символом К . 


границы дУі. множества 

информации об уже построенной границе 
г - выпуклого многогранника в пространстве 
Все многогранники Л 


ч 




одинаковое количество вершин. 
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Таким обравом, многогранники іѵ* ( (і = 0,1,... , /V) можно 
прѳдотавить в виде 

- (ю : Л/-ГЛ 1 [*<]> II Л" к* иК\- 

1 1 Н*1 

^Приведем рекуррентные формулы построения точек х <к) [і^е 

ед " и **&]-***[**,]+ г <м * т * 

* иРти*,)) і'(е*а 

Чиода Г и 4; выбираются согласно равенствам 


г - #й 1 х * х *<♦»]! • ' 7т/ » А '> ,е А** » 

4, - г Г, я ^ мвЛ ^ч>>^о« +см*пі 

При укв ванных соотношениях для нахождения величин Г и 
4 - множѳотво А*У* >+,]) прѳдотавляѳт собой 

многогранный конуо, сопряженный конуоу вовможных направлений 
в точке х (я) [ 1 І4 .,] к множеству IV*. 4) . Это повволяѳт 

вѳоти вычисление коэффициента ^ ы) по формуле 


•М. 


й: тах 


е , А(* І )і™[і і . І )+г , ва і )и;+ т^к ■ 


е'(Е+а к А(ф 


& 


хяе "1 Г (х ОЛ [і і .,])П8“П$Х і -конус 

внешних нормалей в у -й вершине многогранника Р \ 8„ - ко* 

нус внутренних нормалей в я -й вершине многогранника <7 . 
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УДК 62-50 С.Т.Завалищин, О.И.Стародумов 

(ИММ УрО АН СССР) 

ОПТИМАЛЬНОЕ ПО РАСХОДУ ЭНЕРГИИ УПРАВЛЕНИЕ ДВИЖЕНИЕМ 
МАНИПУЛЯТОРА В СРЕДЕ 

Проектирование машин, предназначенных для работы в экс- 
тремальных условиях, является ко.алексной проблемой [і] . Она 
предусматривает ряд специфических задач по управлению механиз- 
мами. Задача об управлении перемещениями манипулятора [2, 3] 


- 29 - 


в среде ставится как вариационная, в которой качество управле- 
ния оценивается по энергетическим затратам, требуемым на пере- 
мещения звеньев манипулятора. Поскольку эволюция манипулятора 
происходит в потенциальном поле силы тяготения, то такая зада- 
ча динамической оптимизации равносильна минимизации энергети- 
ческих затрат на преодоление силы сопротивления среды. 

Поставленная эадача имеет ряд особенностей. Во-первых, 
она нерегулярна [4] : уравнения Эйлера-Лагранжа не содержат 
управляющие моменты и, следовательно, не позволяют формально 
определить их оптимальные значения. Во-вторых, как выяснилось, 
энергосберегающіе программы изменения управляющих моментов 
имеют импульсные составляющіе, поэтому применение классичес- 
ких вариационных средств для решения задачи лишено логическо- 
го основания. Третья особен- 
ность вытекает из второй и со- 
стоит в определении подсчсчета 
энергетических эатрат. Оче- 
видно, для этого надо опреде- 
лить корректный способ умноже- 
ния импульсных управляющих мо- 
ментов на разрывные реализации 
угловых скоростей звеньев ма- 
нипулятора. 

Цель данной заметки заклю- 
чается в применении общего под- 
хода 15] , позволяющего преодолеть перечисленные трудности и, 
в конечном случае, решить обсуждаемую вадачу. 

I. Рассмотрим снвчѳлв однозвенный манипулятор, показанный 
на рис.І. Предполагается, что вся масса т сосредоточена 
в охвате манипулятора. Движение происходит в вязкой среде, в 
которой сила сопротивления пропорциональна скорости перемеще- 
ния массы т . Момент Ц , приложенный к звену, считается 
управляющим, а момент Г - возмущающим. Применяя метод Лаг- 
ранжа, получаем уравнение движения звена 

тг 1 Ѳ ’ъ-тдгсоіѲ *■ II* Г , (I) 



Рис.І. Одноавенный манипу- 
лятор 
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где О с - - аг*Ѳ - обобщенная сила сопротивления среды. 
Выражение для мощности имеет вид 

А=ѲѴ , А(ф-0. (2) 

Задача I. Найти закон перемещения звена из начального по- 
ложения Ѳ(0 ) - Ѳ 0 , Ѳ(О) = О в конечное Ѳ(і^)=Ѳ р . 

Ѳ(ір) <= О , обѳспѳчивапций наименьшее значение работы А({^. 

Дифференциальные связи (I), (2) в такой задаче составля- 
ют систему третьего порядка. Оказывается, что порядок этой 
системы можно понизить на единицу. Действительно, принимая 
обозначение 

ѳ = (з) 

разрешим уравнение (I) относительно управляющего момента 
II = тг 2 о> = тдгсозѲ+аг г 0 - Г. \ (4) 

Подстановка выражения (4) в уравнение (2) дает соотношение 

4 •= тг г ы<1) тдсОсозѲ ог*и>* -Гы. (5) 

Таким образом, выбор угловой скорости (3) в качестве уп- 
равления позволяет ограничиться дифференциальными связями (3), 
(5). 

Однако трудности, отмеченные во введении, сохраняются 
и для новой формулировки задачи. Так, первое слагаемое из пра- 
вой части уравнения (5) очевидно, теряет смысл для разрывной 
программы изменения угловой скорости, что соответствует им- 
пульсному управляющему моменту. 

В данной ситуации примем определение: тг г ыы - 
• Т)(іІ^тг г ы 1 ), где 2) - оператор обобщенного дифференцирования 
[6] . Для гладких функций ы(-) такое определение совпада- 
ет о классическим результатом. Теперь выражение для работы до- 
пускает представление: А » У + //^ тг*ш 1 , где величи- 
на У (•'і - решение уравнения (6). В результате получена 

стандартная задача оптимального управления 

У( — - тіп ; 

о) , Ѳ(0 ) •‘Ѳд ; 
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У(0) = 0 ; 


У = тдгсосозѲ + аг г ш г - Ги) , 

( 6 ) 

Ѳ(і р ) ■* Ѳ р , и)(0) - со(ір') = О. 

Решая последите задачу с помощью принципа Лагранжа [У ] , 
можно получить следующий результат. Оптимальная программа из- 
менения угловой скорости имеет вид 

Ір 

и„ Р і - (аг*У'(г+ ір(аг г (Ѳр Ѳ 0 ) - ] Шг . 

о 

Если возмущающий момент Г = 0 , то видно, что (О орі = 

= і~р(Ѳр-Ѳ„) , т.ѳ. угловая скорость на всем промежутке управ- 
ления постоянна. Таким образом, в ходе оптимального перемеще- 
ния обобщенная сила сопротивления среды сохраняет постоянное 
значение. Согласно первому слагаемому выражения (4) получает- 
ся, что оптимальная програша изменения управляющего момента 
имеет начальный и конечный импульсы с интенсивностью соответ- 
ственно 

тс г и 0 р1 ; ~тг г и) 9рі . 


Далее перейдем к вопросу о синтезе оптимального управляю- 
щего момента. Если следовать идее [2] построения синтезирующей 
процедуры, исходя из оптимальных программ, то формально синте- 
зирующая функция имеет вид (Г=0) 


ѵ(иѳ,ѳ) 


-тг г (д-(ір-і)(Ѳ р -Ѳ(і))д і , 
- тг , Ѳб і/> , і- і р . 


ОіЫір } 

( 7 ) 


"Бегущий" импульс из правой части этого равенства как обобщен- 
ная функния смысла не и.г.ѳѳт. Выход из зтого положения заключа- 
ется. в сужении оператора (7) на последовательность моментов: 

(Ь т е. і т *’) . Очередной корректирую- 

щий импульс осуществляет оброс фазового вектора объекта (I) на 
линию, описываемую уравнением (ір - і т )Ѳ + Ѳ = Ѳ р . 

С увеличением частоты коррекции управляемый процесс будет стре- 
миться к предельному режиму, задаваемому уравнением ( і р ~і)Ѳ + 

+ Ѳ ® Ѳр • Ясно, что для предельного режима 6-0 , т.ѳ. 
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и>~ сопзі «= и) орі . В итоге окаэалооь, что предельный 
режим не зависит от возмущающего момента, т. е. позиционное 
донтульоноѳ управление 
(7) компенсирует неиз- 
вестное возмущение 
(см.рио.2). 

2. На рио.З 
изображен матема- 
тический деухе вен- 
ный манипулятор. 

Как и выше предпо- 
лагается, что связи 
не обладают массой. Рио. 2. Отработка внешнего возмущения 
Движение происходит 

в среде, в которой силы сопротивления пропорциональны скорос- 
тям перемещения масс т ) , т г . Моменты ІЛ , , прило- 

женные к звеньям, считаются соответственно управляющими и 
возмущающими. Уравнения движения манипулятора имеют вид*. 



- 2рѲір - ^ір г +(т 1 -ип г )дг І соз&^ ( 6 ^ 
+ т г дг. соз(Ѳ* р) ■= ^ + Ѵ / 2 + 0 ( с - ^ - Г ж • 


,, , ^ ^ 
ы іг Ѳ + *п \ Р + № * т г9 Г г СО$(Ѳ* р) - Ѵ/ 2 +СІ г - Г 2 


Рис.З. Двухзвен- 
нкй манипулятор 
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Здесь приняты обозначения: 

<* н = (т, + т^с*+т г г*+2п> г г,г е созір ; 

Ы. 1г « т г г* + т г г 1 г і созір , а гг = т г г* ; 

/ - /п а г,г а $1пу ; 

в? ш г 8, ё , 8 І = (а, * а г )г/+ а г (г*+ 2г,г г соз ір ) ; 

4 “ а * г * » 

С 

где .(?. - обобщенные силы сопротивления среды. Выпишем вы- 

ражение для мощности: 

А = (Ц + 1/;)0 + 1/ л ф . ( 9 ) 

Задача 2 . Найти закон перемещения звеньев манипулятора 
из положения Ѳ(0)~Ѳ д , 0(О) = О, ір(0)=\р д , хр(0)=0 в положе- 
ние Ѳ(і^“Ѳ р , Ѳ(ір) = 0, у(і^~у>р, хр(і^)-0 о наименьшей затра- 
той работы Л/*/.) • 

Как и в п.І, здесь можно понизить порядок системы диффе- 
ренциальных уравнений, описывающих динамические связи. Принимая 
обозначения 

9*4), , у> ~ со г , (10) 

можно разрешить уравнения (8) относительно управляющих моментов 
и, подставляя найденные выражения в формулу для мощности, по- 
лучить уравнение 

А=в і со* + $ г 4 - * &>,)&), а) г +С/д, •>- Л7 г ) а), со$Ѳ + 

+ + + (Д ? й), +<* „ й),) Ш 1 +• (II ) 

+ (° ( іі 0) і +< *іг аі і)4 • 

Теперь роль правлений будут играть угловые скорости Д , и) г 
а дифферѳнпиальныѳ связи будут описываться уравнениями (10), 
(II). 
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В выражения (II) входят произведения а) • со- , которш 
надо придать смысл в случае импульсных управляющих моментов. 

Это достигается с помощью следующего приема. Для гладких про- 
грамм со- (•) верно равенство 

(<<„ а) ( + <Л іг іх)^)ЛоР 1 + (ы п а), + ы і2 сд^)Гяд 2 =-ІЖ- V , (12) 

где 

У+ і/2м 1і ^^о( іг со 1 а) г + 1/2 ы гг а)* . (13) 

Правая часть соотношения (12) имеет Смысл а для негладких 
функций со- (■ и, следовательно, может служить определением 
левой части в негладкой ситуации. 

Если искать выражение для работы Е виде (13), то относи- 
тельно величины у получится уравнение 

у = В 1 и)* * 8 г (ОІ + (т І +т^2г,и) 1 ео$е+ ц 4 ) 

+ т г дг ж (и), + и) г )сов(Ѳ+ ір) , у (О) =9. 

Поскольку и>і ( ір) - О , то обсуждаемая задача сво- 
дится к задаче минимизации величины у(і р ) при ограничени- 
ях (10), (14). Такую задачу уже можно решать о помощью класси- 
ческого вариационного исчисления и получить следующий резуль- 
тат. В течение оптимального процесса обобщенные силы сопротив- 
ления таковы, что 

67,* = сопзі , 8 1 , В І (й, с У + ('<?/) г * аті . ( 15 ) 

Тождественные соотношения (15) позволяют сформу даровать 
двухточечную краевую задачу 

Ѳ - *>, , Й, (ір) о), - С 1 , Ѳ(і 9 )~Ѳ 0 , &(і р ) =Ѳр ; {І6) 

2 ‘ С* 

У> -Ч Д ^ =С г , \?ау % , ір(ір) - V ѵ . 

которую надо решить подбором конотант С • . Получающиеся 

при этом программы 90) , ір(-) - оптимальное решение исход- 

ной задачи, сформулированной в момент і 0 

Пусть задача (16) решена для любого момента І е (& 6 ^ * Ір) 
и любого набора [Ѳ(і 0 ) , 4>(і 0 )} из рабочего пространства ма- 
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нипулятора. Тогда, действуя по схеме из п.І, можно построить 
позиционную импульсную процедуру, реализующую в виде импульс- 
но -скользяще го режима перемещение мавдпулятора из исходного 
положения в заданное. В отсутствии возмущения эта процедура 
обеспечивает оптимальнооть перехода по энергетическим затратам. 

Создано программное приложение к задаче, состоящее из 
двух модулей. Первый модуль РОБОТ производит следующие опера- 
ции. 

1. Построение квазиоптимальной программы движения с вы- 
числением констант Г, , С г . 

2. Расчет совокупной величины работ, совершаемых при пе- 
реориентации по построенной программе всеми внутренними и внеш- 
ними силами системы. 

3. Расчет совокупной величины работ, совершаемых всеми 
оилами при поочередных переориентациях звеньев, осуществляемых 
по оптимальным программам (делается это для демонстрации отно- 
сительной эффективности построенной программы движения). 

Тѳкот модуля написан на языке ГОКЕ К ОС ДОСЛАН. Вход- 
ными параметрами модуля РОЮТ явл-ются: 

1) лі, , т г - массы первого и второго звеньев, соответ- 
ственно; 

2) Г 1 , г г - длины первого и второго звеньев, соответ- 

3) Ѳ в , у> в - исходные обобщенные координаты оистѳмы; 

4) Ѳ р , \р„ - конечные целевые обобщеннее координаты системы; 

5) ііте - располагаемое время движения; 

6) а, , Од - коэффициенты сопротивления среды переме- 
щению первого и второго звеньев, соответственно. 

Выходные параметры модуля РОБОТ включают: 

1) С і С д - константы, определяющие программу движе- 

ния; 

2) совокупные работы всех сил пр: поочередном перемещении 
(ЗН 1 ) и при совместном перемещении ( ЭН1 ) звеньев по опти- 
мальным программам. 

Второй модуль СИНТЕЗ производит синтез алгоритма импульс- 
ного позиционного управления движением двухзвенного манипулято- 
ра как при движении без внешних возмущений, так и при возмуще- 
ниях (синусоидального, либо импульсного характера). Текст мо- 
дуля написан нѳ языке ГОКЕХ 
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Входными параметрами модуля СИНТЕЗ являются: входные 
параметры модуля РОБОТ, дополненные параметром РЕЖИМ, опреде- 
ляющим наличие или отсутствие ( О ) внешних возмущений. 

Выходными параметрами модуля СИНТЕЗ являются: 

1) фактические координаты системы в начальный и конечный 
моменты движения - Ѳ % , у д , Ѳ р , у> р ; 

2) ііте - фактическое время движения системы; 

3 ) ЭН - совокупный энергетический раоход. 


Данные теста 


/7?, 

= 30,00, 

кг 

т г = 

1000, 

кг 

Г 1 

= Д50, 

м 

Ъ = 

тро, 

м 

9 0 

= 3000, 

град 

% - 

-6ф0, 

град 

ѳ р 

= І5ДЮ 

град 

У>р « 

7Ц00 

град 

ііте 

ІІ 

О 

0 



РШИМ = I (без возмущений) 




а , 

= доо 

Выход (РОЮТ) 

а г = 

ДОО 


С. 

= 0.15 


С г - 

4.35 


ЭН 2 

—1.23 


ЭН1 = 

10.25 



Выход (СИНТЕЗ) 

Ѳ 0 * 30.00, град \р 0 = - 60.00, град 

Ѳ р =150.02, град хр р = 70. 00, град 

ііте = 2. 98 , о ' 

ЭН =2.21 
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УДК 517.938 Ю.И. Бердышев (Институт мате- 

матики и механики УрО АН 
СССР) 

О НЕОБХОДИМЫХ УСЛОВИЯХ ШШМАЛЬНОСТИ ПО БЫСТРОДЕЙСТВИЮ 
В ЗАДАЧЕ ПОСЛЕДОВАТЕЛЬНОГО ОБХОДА ГРУППЫ ТОЧЕК 

• 

Иоследуѳтоя задача оптимального по быстродействию управ- 
ления нелинейной системой о’ целью сближения ее в указанном по- 
рядке о заданной группой точек. Получены необходимые условия 
оптимальности управления и моментов сближения в форме принципа 
максимума I. С. Понтрягина Гі] , не использующие декомпозицию во 
времени. 

Пусть движение системы в п -мерном евклидовом пространст- 
ве Р п описывается векторным дифференциальным уравнением 

хеР", иеР (I) 

с начальным условием х(і 0 ) = X 0 . Здесь Ра Р г 

компакт ( г - размерность пространства управлений) , функция 
/ • [і 0 , <* *” [* Р п * Р -*■ Р непрерывна по совокупности аргумен- 
тов вместе со своей частной производной по г .В простран- 
стве "геометрических координвт" Р * задана совокупность 
точек 
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(2) 


{ и (і) е Р н , / е /, П) ], 

где к, т - натуральные числа, к 4 п . Обозначим через 
я оператор проектирования Р п в Р* .В качестве до- 
пустимых управлений выберем V - множество всех измеримых 
по Борелю функций 11‘- [і д ^^ [ — - Р . Каждому управлению- 
программе I] е. V соответствует ір ц 4 (ш (і) І і > і а ) - 
движение оиотемы (I), порожденное этим управлением из началь- 
ной позиции (І д , Хд) , й - равно по определению. 

Будем говорить, что управление II е У обеспечивает 
сближение системы (I) с совокупностью точек (2) за время 

в заданном порядке VI й { /, і? , . . . , /77 ] , если най- 

дутся такие моменты времени г, , . . . , г т , удовлетворяющие 
условиям г, < . .. и г ті т т ‘іІ, при которых имеют место соот- 
ношения 

\\*(ч> ѵ ('і))-ѵ т І = о, а*К"), о) 

где (| • II н - евклидова норма в Я * . 

Задача I . Определить управление 11 ° е V , обеспечиваю- 
щее сближение системы (I) с совокупностью точек (2) за наимѳнь 
шее время 1 9° - і д . 

Пусть 

Ті {г : г е Р т+< ,т в ~і в , г,„, -с г- Це />)) ; 

Г„ *{*:*« Г, г.-*}; 

5 т Ь[л:АяК п , а^О(іе і >), П - /}. 

Задача 2 . При фиксирован ном м оменте Т т ■= і} , глад 
ких функциях : Р п Я (іеі',т) определить пару 
(Т* 7 Ѵ*)е "Г# х У , доставляющую минимум функционалу 

Г* <4 V) “ П Ы і ( Уѵ >'*;>) (* е5 т). (4) 

Будем предполагать выполненными условия неограниченной 
продолжимости решений системы (I) 
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Ч«’*Л 6С(/+ІХІ Я ), С- соп$і , 


, п+1 


а также выпуклость в Р множестве обобщенных скоростей 

Ѵ^,х)- [!(і } т } и), иеР } 

при любых фиксированных (і,х) . Эти уоловия гарантируют 
([2] , о.284 ) существование оптимального управления II в 

палате I. 

Введем обозначения : (7) - \р ѵЛ (і) - движение, порож- 
денное управлением ІІ*; 5< ({, ^), фундаментальная мат- 

риц решений оиотемы: 

сЦбх) д/(і,х,и) д_ 

-гг V— " 

(5) 

* 4 {*:*«*", |г|, *Г], 

С* ~(дв;(Х)/дх)\ , 

' ІгТ-и^г*) 

Положим 

V» 6 /, л/ , ^ е , г,. ] ; 

(6) 

("штрих" означает транспонирование). 

Справедливо [з, 4] олѳдующѳѳ утверждение: 

Теорема I . Пусть пара (Г* 1/*)е Т# * V является реше- 
нием задачи 2. Тогда для каждого і е. Т~іп почти всюду на 
отрезке [г_ м Т. ] имеет место соотношение: 

ч (і), Ц*(і)) - пнп ФлЩіи ч ГО , «) , (V) 


а в точках х* (і е /, /т?-/) выполнены условия 

V Г Ѵ)/ГГД 4,(7^), и) - 

V < ‘*’ ) (г *)}(?■ , “)• 


№) 
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Заметим, что эту теорему можно использовать при функциях 


&: :/Г— Я, і е 


Г,т 


определяемых формулами 


&/т)йІ^(х)-ііг т І А 


(9) 


лишь тогда, когда 

Ѵі е І^ті : <3.(ЬГ Л (г*)) * О. (10) 

Это связано ^необходимостью вычисления производных д&^дх- 
(іе/,т , у'е /,х) в точках х^- 4 ( ?*) . В дальнейшем 
мы укажем один из способов выбора вектора <Л в выражении 
(3), при котором теорема I применима и в олучаѳ, когда функ- 
ции (іе. /,т) определяются формулой (9). 

Теперь при заданном моменте рассмотрим две сле- 

дующие задачи: 

(г ($ Т * V ' г ' гГ) ^ ш ' в,ві *' {Ш 

тах К], Ы,т, V) — тіп , (г, 1/)е Г# * V. (12) 

“* в 5 т 

Последняя из них тесно связана с задачей оптимального быстро- 
действия (задачей I). Действительно, в силу тоадеотва 



решеіше задачи I может быть получено как предел решения зада- 
чи (12) при стремлении момента & олѳва к наименьшему 
моменту времени 1?^. , при котором значение задачи (12) 

равно нулю. Теким образом, время - і„ .равно времени 
оптимального быстродействия 0 °~І о в задаче I. 

Можно показать, что вектор <ч 0 , являющийся решением 
задачи (II), обладает следующим свойством. Если (г* 1 , 1/ л ) - 
решение задачи 2 при и 
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( 14 ) 



то Ѵу е 1,т I М 0 : ы" - О . Использование вектора ы. 0 позво- 
ляет исключить из рассмотрения те точки , для которых 

<*і т О . Это дает возможность при решении задачи 2 и ы - «с 0 
применять теорему 2 в случае, когда функции Р(іеІ^т) 

определяются формулой (9), предварительно положив т~р (/> - 
число элементов множества М 0 ) и перенумеровав точки у іеМ 0 . 
Иначе говоря, при решении задачи (II) можно использовать тео- 
рему 2. При этом точки ъ/ (, \ і е 1,т /М в исключаются из 
рассмотрения, а сближение с ними происходит "попутно" при ре- 
шении задачи о сближении о точками V м , і е М 0 . 

Ив теории игр иэвеотно, что олтимумы задач (II), (12) - 
величины , определяемые формулами 


& рюх 


тіп 

(Т,Ѵ)еТ„. V 


Га Н» г « и )і 


в'* 


тіп 

(г ,Ѵ)*Т ѵ *Ѵ 


пах 

•** 5 П) 


удовлетворяют неравенству 

* Г* • (15) 

Очевидно, для достаточно малого *#) аначения 

положительны. Пусть , і? г - наименьшие моменты времени, * 
в которые зануляются соответственно величины 0® , • Ранее 

отмечалось, что ^ і д - время быстродействия в за- 
даче I). Теперь покажем, что і} в - . Это равенство позво- 

ляет искать г}° как наименьший момент времени т? , при 
котором величина обращается в нуль. Вначале ваметим, 

что соотношение 


г?# * . (16) 

вытекает непосредственно из неравенства (15). При гТ*. г^. 
найдется такой вектор ы. и соответствующая ему па- 

ра (т ы , и* ) - решение задачи 2 при <Л = <л° : 
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. ѴІ е /> : Ѳ- ( хр ѵ <е С гГ)) ' О. (17) 

Это связано с там, что лишь при уоловиях (17) величина 

Г*, & (пах &і(ѵ> л *( г/)) 

Г < е Гт 

обращается в нуль (как того требует выбор момента ). Иа 
(17) вытекает неравенство У° 6 і? г . Отовда о учетом (16) име- 
ем і? 9 = . 

Итак, каждому моменту «Я , меньшему , соответ- 

ствует определенный вектор «с 0 и разбиение множества /, т 
на чаоти М 0 (14) , 1,т /М 0 , при которых 

V' е М 0 . (18) 

поэтому следующие определения векторов ?; (&) , которые мо- 

гут быть использованы в формулировке теоремы 2 в качестве век- 
торов , являются корректными. Эти соотношения позволя- 

ют при определении оптимальною управления II* в задаче 2 
с функциями &■ (I е /, т) вида О) и моментом окончания 'У , 
мѳныд.м , использовать теорему I. Всякой последователь- 

ности ( моментов окончания, схсдящейоя 
слева к Л. отвечает последовательность управлений 
([/" , у» 1,2 , ...), удовлетворяющих необходимым условиям оп- 
тимальности (7), ( 8 ). Поскольку іУу - У ( 1} °- і 0 - втомя быстро- 
действия в задаче I), то при достаточно большом номере У уп- 
равление V; обеспечивает попадание оиотемн (I) в моменты 
времени (се 1,т) в сколь угодно малые б -окрестнос- 
ти точек Ы (і) (2). Поэтому оптимальное по быстродействию 
управление 1/° можно определять как предел (в смысле * - 
слабой сходимости) последовательности (Ѵ $ ,5*1,2,...) при 
У — іУ 0 управлений, каждое из которых удовлетворяет необ- 
ходимым условиям оптимальности в задаче 2. В случае, когда 
система (I) линейна по и , предельное управление выписыва- 
ется в явном виде. Б силу выбора функций €)■ (і е і,т) (9) 
их градиент в точках ( Т при условиях (10) имеет 

длину, равную единице. Поэтому определение функций (6) при 
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любом Ф < является содержатѳлькш, а соотношения, 

фигурирующие в теореме I, невырожденными . Обозначим через 
[ ( (і е /, т) предел I* (5) при і? 0 , а через 

У ^ ( і) - чму соответствующую функцию , Поло- 

жим 


V 1 * 3 » 


т 

с 


Ѵ-Лб ■ 


Теорема 2 . Пусть Ѵ° я. V - оптимальное управление в за- 
даче I, е* , . , , , г° - соответствующие моменты сближения 
системы (I) о совокупностью точек 12), цг(і )= у> ѵ „(і) . 

Тогда для каждого і е /, /я почти всюду на отрезке [Т 1 _ 1 , г ; ] 
имеет место соотношение 


Ф>)у & ??'» * л? ^ (і ' ѵ(і) » ") > 

а также V/ е /, т-І 

тіп иг( г/), «) - 

9 /* ^ (г *,ьх(\ гД м). 


Пшмер . Рассмотрим задачу об оптимальном по быстродейст- 
вию сближении управляемого объекта, описываемого нелинейной 
сиотемой 

X = С08 ір, 

у тешу», . (І9) 

ір-Кі/і <и| * 1 , Н • еоп5І 

о заданной системой точек (ы и) е /?*/«? /,/Л } . Система (19) 
является простейшей моделью летательного аппарата, имеющего 
аэродинамически!, рули и движущегося в гори зонта льной плоскос- 
ти. В рассматриваемом случае: 




1 0 <!о$ \р(і)~ео9ір(і д ) 

О 1 5іпу(і) - зіпір (ід) 

О О I 
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( 20 ) 


(созір (і) - сов ір ( * (апір(і)- ап у> (*,))'; 
и °(і)- 5ідп [ Е (еозу(і) - соир(і 0 )) + 

+ Лр(зіпу>(і)-зіпу>(І в ))] , іе[ г { .,,г„]. 



Можно показать, что из теоремы 2 следует справедливость соот- 
ношения 

^7 «к ѵ 1 ( О + Си 8іп і> (*') “ 0 ( 1 в й™) (21) 

на оптимальной траектории. Иэ [5] вытекает, что на отрезке 
[І' • г л» ] оптимальная траектория ооотоит из душ ок- 

ружности и отрезка прямой. При этом величина Ѵ д і ' п ' 1 ( і') должна 
быть положительной при любом ц>(т)е]\р(Т т )-2/г, 
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Отсюда имеем 

Л т1 - • - со$ч>( с т ) , Л т1 ш~ $сп ч>(г т ) . ( 22 ) 

в , ■ — 

Можно показать, что координаты л { , і »= 1,т вектора 
ы° - решения задачи (II) в рассматриваемом случае удовлет- 
воряют соотношениям (см. рисунок) 

2сі° т соа^ - Л; (і 6 /, т-І). (23) 

Кроме того, по аналогии с [5] можно показать, что оптимальная 
траектория состоит из дуг окружностей радиуса ^ и отрез- 
ков прямых, их соединяющих. Ввиду (21) выбор вектора ^ пол- 
ностью определяет дугу окружности, проходящей через точку ы (1 \ 
Заметим также, что в точках сопряжения дуги и отрезка соот- 
ветствующие выражения, стоящие в квадратных окобках в (20) , 
обращаются в нуль. Эти замечания и соотношения (21 )— (23) ис- 
пользуются при построении итерационного метода вычисления век- 
торов л а , (*, і е Г^т , определяющих оптимальное управление 
( 20 ). 

На рисунке изображена оптимальная траектория . і орі , по- 
лученная о использование ,і указанного метода при следующих на- 
чальных условиях: 

X, = 15, X, = О, Х,= 60, X, = 40, 

У 0 =-40, у і * 0, у і = 10, У} = 55, 

где х ; , у і - координаты точки ь/ (і) , / е /, 3 . Сравним тра- 
екторию І дрі и тра кторию , полученную в резуль- 

тате декомпозиции задачи с т целевыми множествами на т 
двухточечных задач. Из риоунка видно, что время движения по 
траектории существенно больше времени движения по 

оптимальной траектории. 
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УДК 517.938 В.П. Серов, А.Г. Ченцов (Уральский 

политехнический институт) 

КОНЕЧНО-АДДИТИВНЫЕ МЕРЫ И КОНСТРУКЦИИ РАСШИРЕНИЙ В 
ЛИНЕЙНОЙ ЗАДАЧЕ УПРАВЛЕНИЯ С ИНТЕГРАЛЬНЫМ 
ОГРАНИЧЕНИЕМ 

I. Введение . Рассматривается линейная система 

і(і)-А(І)х(і)+и(і)ё(1), (І.І) 

где х(і) е X к К п ; матрица-функция А ( О , определен- 
ная на Р , непрерывна по ^ ; управление иО) : Г— Р 

кусочно-постоянно и непрерывно оправа ( Т А [і , іХ 0 [ , где 

і 0 < і7 д ; Р - чиоловая прямая). Управление и (Л стесне- 
но ограничением ^ 

])и(і)І(И *с , съО. ( 1 . 2 ) 

Вектор-функция 6 (-) ' Т -» X не предполагается непре- 

рывной (конкретные условия, которым должна удовлетворять в (•) , 
будут сформулированы ниже). 

Полагая ІА[.і в ,^’ в ] = ТіІ (гТ„) , рассмотрим мно- 

жество X * 1 всевозможных функций / : / — X , наделенное 
тополо-иѳй поточечной сходимости [I, 2] . Будем рассматривать 
траектории системы (І.І) с заданным начальным уоловием (н. у.) 

Т ( і д ) => X 0 как элементы X 1 . Получающийся в резуль- 
тате пучок "обычных” решений не является, однако, компактнш 
множеством; управления и (•) , удовлетворяющие (1.2), также 
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нѳ образуют компакта в смысле какой-либо естественной тополо- 
гии. Более того, отображения вида 

'і — и(і)в(і): Г — X 

имеют "в пределе" омыол своеобразного произведения обобщеннной 
функции [3] на раврывную, а это представляет определенные за- 
труднения при использовании аппарата обобщенных функций для 
построения расширений. В связи с этим в данной работе будет 
использован другой аппарат, а именно, аппарат конечно-аддитив- 
ных мер (КАМ) [4, гл. III, ІУ] , что позволяет решить пробле- 
му существования оптимального управления в целом ряде экстре- 
мальных задач. 

Отметим, что исследования в области задач импулі'ного уп- 
равления линейными системами, базирующиеся на принципе двой- 
ственности, восходят к фундаментальным результатам Н.Н.Красов- 
ского, связанным о проблемой моментов [5, б] . Развитие подхо- 
да [б] получено в работе [7] ; см. также [8, 9] . Методы рас- 
ширения для задач оптимального управления и вариационного ис- 
числения подробно изложены в [10-12] . Идеология рассматривае- 
мых ниже конструкций расширения связана также о постановками 
[13 - 14] . 

2. Конечно- ~ шштивные меры . Еоли И -множество, то через 
Р(Н) обозначим семейство всех подмножеств Н . Пусть 

2 й й Р(Н)\ [<р] ; Уй [і ; 2; ... } ; 

: (Тйй {; : ,/е/У, 

Через Т обозначим семейство всех полуинтервалов [с, с( [ 
таких, что се Т , СІеІ . Тогда ( Г, Г) - измеримое про- 
странство о полуалгеброй [іб, с. 46] множеств. Множество всех 
положи темных КАМ ^ : Г — К обозначим через (ааісі')-* [Г], 
следуя [13, 16 ] , Пусть бе (асій) + [Г] - длина [13, 

с. Зб] как функция множеств, а именно : Ѵс е Г, е /с, 1^ ] 

С (ССу(іС) & (1-С . Через В а (Т 1 Г) обозначим множество 
всех ступенчатых отображений (Т,Т) в И? . Пусть, как и в 
[4, 13, 16] , В(Т,У~) есть замыкание В 0 (Т 1 1~) в смыс- 

ле метрики равномерной сходимости, Пространство В (Т , Г) 
с традиционной $ир -нормой [4, с.26і] является банаховым; 
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пространство, топологически сопряженное к В(Т,Г), обозна- 
чим через В*(Т } /у . В обычном смысле [4, гл.У] понимаем 
термин * -олабая топология 8 *( Г, Я : топология в 
порожденная элементами В (Т, .Г) . Пусть теперь 
А Г«Г) - множество всех мер у» - У , когда уи и V про- 
бегают (а(1<І)+ [.Г] . Наделяя линейное пространство А (У) 

* (КАМ ограниченной вариации) естественной нормой-вариацией , 
мы получаем банахово пространство, изометрически изоморфное 
8 Ц (Т , Я . Соответствующий изометрический изоморфизм 
ДГЯ на В ( Г, Я определяется как отображение, переводя- 
щее уи е А ( Я в функционал . 

д •— I ?(і)уи(с*і): В(Г,Г) ~ В ; (2.1) 

т 

интеграл понимается в смысле [іб, с. 75] . С учетом изометри- 
ческой изоморфности Д (*Г~) и В*(Т У <Г) наделим само 
А (&)- слабой топологией г, (У ) , вкладывая в это понятие 
смысл * -слабой топологии В *(Г, У) . Полагаем далее 

<С] ; (2.2) 

П а й{//:^е (асШ)+ [У] у ди (т)ла] (а* 0 ), (2.3) 

Народу с (2.1) рассмотрим неопределенный интеграл [16, о. 76 ] . 
Именно, Ѵ’/е А (Я,/ е В (Т, у) через ^ * ѵ обоз- 

начим меру из Д (ІГ) , сопоставляющую каждому множеству і е У 
уі/ -интеграл , ^ на множестве і . Пусть Ѵа е [о, *>[ 

М а й {$:]вВ*(Т,Г\ *.а), (2.4) 

г 

где В* (ту)- положительный конус Я в омыолѳ поточеч- 

ной упорядоченности. Между множествами (2.3)*и (2.4) существу- 
ет тесная связь [17] , сводящаяся к возможности аппроксимации 
мер из (2.3) неопределенными интегралами элементов (2. 4). Здесь 
мы учитываем, что і(, г ) + 0 при ГеУ\[ф] • Поэтому 
каждая мера из (2.3) является С -непрерывной в смысле [17]. 

Если ( Н, г) - топологическое пространство и Ае. У (Й), 
полагаем, что сі? (А , г) - замыкание Л в (М,т) . Тогда 

[17] Ѵа е [0, 
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(2.5) 


^ ГГ/ : } е ^а) і г « #1) • 


Обозначим черев 2/ множество всех 1/еВ а (т,^) таких, что 


1)Ѵ(і)\еыь* 


с. 


Легко видеть [13, гл.і] , что V - множество всех ступенча- 
тых отображений ( Г, & 1 ) в К , соблюдающих ограничение (1.2). 
Используя (2.5), можно показать справедливость следующего ут- 
верждения. 

Теорема 2. 1 . Множество 5 является * - олабым замыка- 
нием множества [V * 8 • I? е II ) : 


3-с4((1/*е : [ГеѴ),Т* (А). 


3. Обобщенная управляемая система . Рассмотрим систему 
(І.І). Будем предполагать, что Ѵіе/^п і - я компонента век- 
тор-функции 8(‘) является элементом В (Г, У) , т.ѳ. 

6- е В (Т у /*) .В этих условиях 8- (•) являются ограничен- 
ными борѳлевскими функциями, поэтому система (Т.І) имеет един- 
ственное решение (в смысле Каратэодори [10, гл.ІІ, п. 4 ] ) 

К ^ (Ѵи (*) е * ,Ѵ Ѵ , и это решение 

может быть эаписано посредством формулы Коши : для V 9 е. I 

V Ѵ(рР(Ѳ^Щ)ІЫ^). (З.І) 

[і в ,в( 

Здесь ѳоть фундаментальная матрица-функция решений * 

системы х •= А (і)х ; интегралы в статье вычисляются поком- 
понентно. Подчеркнем, что в (З.І) мы можем вернуться к С -ин- 
тегрированию [13, § б] , так как компоненты подынтегральных 
функций являются элементами В(Т } >?) ; см. по этому поводу 

[4, с. 182] . ПуотьѴ,д { у>„ : I/ е 2/] - пучок обычных ре- 
шений (З.І). Далее полагаем У^/и е 5 ,что \р^ й ( ^ (і), 

4 ( « і^) еХ 1 есть отображение 


Ѳ—Ф(е,і в )* а + | (<*$): 1-Х^, ( 3 . 2 ) 
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где интеграл по КАМ понимается в соответствии с [іЗ, глЛі] , 
[іб, о.7б] . Пучки V, и уі/е5“] рассматриваем 

как подмножества X 1 о топологией ? поточечной схо- 
димости. Отметим, что V, с V, , так как V- 1/ е V , 
д& Ѳ(Т, \Г) , ГеУ 

[дЮ ѴШШ) “ [§(*) (ѴѣШО , 

г л 

где 11*1 е. 5 . Таким обравом, "обычные" решения суть чаот- 

ный случай обобщенных. Ив (3.2) и определения Г, (/) следу- 
ет 

Лемма 3.1 . Отображение уи ’■ 3 ~ ■ X непрерыв- 

но как отображение ("5, г ^ (У)), Г* 1 ("./) - {5П6 •' С 6 г, (^)) , 

в (XV)- 

Теорема З.І . Множество Ѵ г и замыкание V, в ^Х г , г) 
совпадают. 

4. Задачи оптимизации с Фазовыми ограничениями . Пусть 
у': X 7 — >? - функционал, непрерывный относительно 

(Х г , г) , и пусть задано отображение і — -/V* : Г—^(Х) , для 
которого /Ѵу замкнуто в X при Г . Обычная зада- 
ча оптимизации имеет вид 

!г(у>ѵ)- іп /і ѴеѴ > УУ'*> еЛ і (4Л) 

Она, однако, может быть неустойчивой при изменении /V , в 
связи с чем целесообразна асимптотическая постановке в надлежа- 
щем классо приближенных решений. Рассмотрим, следуя [10,гл.Ш] , 
обобщенную задачу 

ЦѴ/и) — тіп »/ /6 5 < (* Л Я- < 4 -2) 

Теорема 4.1 . Пусть ограничения 8ѳдвчи (4:2) совместны: 
множество 

Г» (4-3) 


непусто; тогда Зуи 0 еГ. Ѵ^еГ 1'(Ч > ^/ 0 )* I ($//)■ 
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Рассмотрим теперь асимптотическую постановку задачи 
(4.1). Обозначим V е е I е I через /Ѵ/ с; замкнутую 

евклидову е -окрестность /Ѵ і , Л' сХ . Пусть Гіп (!) 
есть семейство всех непустых конечных подмножеств 1 
Условие 4.1 . Ѵс е ]0,^[ , Г е Гіп (1) : 


и е (Г) й[іЛѴ*Ѵ, (Ѵі е А '.у ѵ те /V* 1 )] + ф . 


Тогда семейство °Ц всех множеств II ^ (К) , когда 

е пробегает ]0, <*> [ и А' пробегает Гіп ( 1 ) , есть ба- 
8ио фильтра V . При этом [14, § 7] 


(и-аа)[г]лтр Ч г(ѵ„)- 

с* иі І/еЕ 


(4.4) 


$ир Іп/ Г(Ѵ> ѵ ), 

л,#)в]о,~[пПп(Л ѵши^к) 


где у е имеет вид {(1/)й X ( 1 Р Ѵ ') , есть асимпто- 

тическое значение задачи (4.1), достижимое в классе (прибли- 
женных) рѳшений-направлѳнноотѳй [і , 14, 17) . 

Тѳооема 4.2 . Величина (4.4) является значением задачи 
(4.2): 

(и-Шід)(Г)‘ ып. Г($*)- 


Решения уи в вадачи (4.2) позволяют построить оптималь- 
ные приближенные решения-направленности путем аппроксимации 
А г в (\(Г), г, ГЛ). 

5. Двойственность . Пусть т е/ ; /%, . . . , М т 
непустые выпуклые компакты в Г г (г е /~п) с опорными 
функциями р (• \М ,) ,... ,рО\М т ) [12, с. 57] { л і * 0 , .. ., ьО , 
и пусть заданы момелты времени 

^ О “ К і 6 ' * ' * Г />7 * ^ 0 • 

Раосмотрім функционал 


* «.о, 
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гдэ сі (■, М-) - функция эвклидова расстояния от точки (из 

Я г ) до множества М- , а [ у ) г - вектор первых г 
координат вектора уеХ . Тогда 



есть значение задачи (4.2) в условиях, когда /Ѵ^ * X. . Пусть 
2 & 5 оптимально в (5.1): {Орд) - Го • Покажем, что 
2, удовлетворяет принципу максимума Л. С. Понтрягина. Пусть 
і -замкнутый единичный евклидов шар в Р г с центром в 
"нуле". Для любого Л е і положим 

ѴШЁ'Ч (Л‘«)(<Р( г. , І 0 )Х 0 ) С - р (лмЩ)) - 

(5.2) 

- с тал ыр ІЕ <*;Л'(і){ф(г , , *)#0] г I ■ 

$еі,т и[г,. п т г [ *"* 

Теорема 5.1 . Пусть Л д е I. - максимизирующий вектор в 
8а даче У* (Л) — ■ тах , Деі . Тогда 

г: [ г * л' 0 с і ) [<р(п , о №}г й т = 

у - 1 і*$ 

= -с тах іир 

(еі,т 

В заключение отметим, что в данной работе рассмотрен для 
простоты случай скалярного управления; в [18] аналогичные ре- 
зультаты приведены для задач с векторным управлением в услови- 
ях, когда ограничена сумма полных импульсов координат управляю- 
щей функции. 
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УДК 517.977.5 А.Н.Сѳсѳкин (Институт мате- 

матики а механики УрО АН 
СССР) 

МИНШ "НАЦИЯ ФУНКЦИОНАЛА С ИНТЕГРАЛОМ ЛЕБЕГ А-СТИЕТЬЕСА 
ИМПУЛЬСШМ УПРАВЛЕНИЕМ С ОГРАНИЧЕННЫМ РЕСУРСОМ 

I Рассматривается задача минимизации негладкого интеграль- 

ного функционала вдоль траекторий линейной системы при ограни- 
ченном ресурсе управления. Получено выражение для определения 
оптимального значения функционала качества. Найдены необходи- 
мые и достаточные условия для определения оптимального управ- 
ления. Различные эадвчи оптимального управления рассматривались 
в [І-<*І . Задачи минимизации, аналогичные изучаемым здесь, но 
при мгновенных ограничениях, исследовались в [5] , а в дискрет- 
ной постановке - в [б] . 

Пусть движение объекта управления описывается линейным 
дифференциальным уравнением 

X - А(і)х + В(і)ѵ(і) , (I) 

где X - п -мерный фазовый вектор; V- т -мерное управление 
(ѵ&І) (-/а)) ; А(і) и В(і) - соответственно п*п 

и п * /и матрицы-функции о непрерывными на [{„ , ] элемен- 
тами; I/ (•/ а) - класс функций ограниченной вариации оо зна- 

чениями в А * т , удовлетворяющих неравенству: 

ѵаг ѵ ( •) = П ѵас V: (ѵ) 4 а (а >0) , 

в,.*.] *'•' ьлі ? 

где V- - і -я координата вектора ѵ . 

Производные в (I) понимаются в обобщенном смысле [7] . 
Согласно формуле Коши решения уравнения (I) можно пред- 
ставить в виде: 

і 

х[ѵ({)] *Х + 1 і)В(ЯСІѵ(ё), (2) 

где Х((,$) - фундаментальная матрица системы і ~А(і) X 
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В /?* , ГДѲ К - цѳлоѳ число, удовлѳ тв оряюще ѳ неравен- 
ству , для каждого і ^ [і 0 , & 0 ] задано не- 

пустое вшуклоѳ множество М (і) . Будем предполагать , что 
многозначное отображение ^ —■ М (і) удовлетворяет следующе- 
му условию: для любых і, и і г из отрезка [Ь д , & в ] , 

подчиненных неравенству 

/і (М(і,), М(і г )) 6 , (3) 

где ір(і) - неубывающая на [і д , & 0 ] ограниченная функция; 

Н (И, , М^) - хаусдорфово раостояниѳ мѳаду множествами 

[Ь] . 

Пусть С(і) - П * К - матрица-функция, элементы ко- 
торой - функции ограниченной вариации; ос (і) - неубывающая 
на [і д , & р ] ограниченная функция. 

Введем обозначение 

д*<о»,ѵС 0 )« I у(с(і)х[ѵ(-)] , М(і))с(ы.(і), (4) 

где <?(у, М) - функция расстояния от точки уе Я* до 
множества МеЯ* [8] . Интеграл в (4) понимается в смысле 
Лѳбѳга-Стилтьѳса [9] . 

Покажем, что подынтегральная функция в (4) есть функция 
ограниченной вариации. 

Согласно [8, с. 45] 

I -$>(**, Л/)| 4 |г, -2Г,|, (5) 

где 7, и 7 Л - точки, а М - множество из Я* . Из опре- 
деления хаусдорфова расстояния [8, с. 75] и расстояния от точки 
до множества [8, с. 45] следует справедливость неравенства 

(6) 

где 7 - точка, а М, и М г - множества из # к . Из цепочки 

неравенств: 

п , м(ь и ,)) - , м(Щ « 
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♦ {= ІР(*(Ь Л , "(*,))■ -9(»(*і), "(®\ < гѴ* 1+ ,№Ь 

+ П |гСО-*^)І * Ч> (*>,)- Іра^ ѵаг^ »(і), 

построенной с помощью (3), (5)-(6), следуе* ограниченность 
вариации подынтегральной функции в (4). 

Задача I . Требуется минимизировать функционал (4) по 
V е ѵ('/а) вдоль траекторий системы (I). 

Функционал (4) в силу определения функции расстояния й 
опорной функции [8] можно записать в виде: 

* 

ѵ-(о( , ѵ(-)) - Г тах Г, (Я, ѵ(-)) сіа(і) , (7) 

; #лі*< 

'б 

где 

Г,(Я,ѵоі)- Я г С(і)т[ѵ(-)}-?„ (і) (Я), 

' х[ѵ( •)] - решение уравнения (I); у?- * -мерный вектор, 

(у?) - опорная функция выпуклого множестве М . 

Операция взятия максимума в (4) порождает многовначное 
отображение А , ставящее в соответствие каждой тс же 
і е [і 0 , г? в ] некоторое множат гво значений Л ив оданич- 
нохо шара Ц ЯЦ *• 1 . Согласно [10, с. 175] отображение бу- 

дет ы. -измеримым. Тогда в силу [10, о.І7Э] существует Ы. - 
измеримая однозначная ветвь Я из Л , что 

р(С«)х(і), М(і1) ~Я т (і)С(і)х[ѵО)]-о т (Л«)) . 

Обозначим через і. множество <Х -измеримых веятор-функ- 
вдй Я (') со значениями в единичном шаре ЦЯ к 6 / 

После подотановш в (4) X [ѵ(-)] из (Я) задача I 
примет следующий вид: 

Задача I* . Требуемся вычислить: 
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I » м 


• /л/ 

і'т Ѵ('іа ) 


:ПОХ 

Л(-\ еі 


ѵЫ,Л } V ), 


где ^ 

V ^ , Л, V) */ *[Я (і)С(і)К (і, і в \ х *- р* л) ('Л <і *>)] Л* гл +■ 

*4 ; *' 

+ Г I Л т (і) С(і) Х(і у з)& (і)4ѵ($) <Ы(І) . 

*» <„ 

Поменяв порядок интегрирования в двойном интеграле соглас- 
но теореме ФуСини [9] , получим: 

і», 

ѵ{ы % Я, ѵ)-( [Я г (і)СШ(і. Ох а -? ЫІ) А(і))}аІ«(і)+ (8) 

V /' і'я т (і)с№(і % *)*<і№(Ъ<іг(і) . 

*• * 

Теперь с помощью теоремы Фанъ-Цзы [Ті] задаче I* поставим 
в ооотаетотвиѳ двойственную задачу, В идейном отношении рас- 
сматриваемая конструкция родственна [і, 2, 12] . 

Ив вогнутости иг (* , А , ѵ) по А л линейности 

по V следует равенство соглаоно [іі] : 


у в (и)*іп/ ,пах ѵг(ы,А ѵ) - шах /л/ 'ш(^,А,ѵ). (9) 

ѵлч(-/в)*чві. Л0)м М-\лі. 

■Ле мма X . Для любой си о темы функций 


^ О)еВѴ/’* 0 ,і> о ) (ік^Тп)- и іГ(-)еѴ(-/а) 

справедливо ровѳнотво 
т }*• ■ 

іп/ V. 1 1 (г)ау (в -~а таз сир //. 7 »/. 

ѴшѴ(-Іо) ‘ ^„, 9 ,] 


очевидно спрѳвѳдли- 


Докавательство . Для любого і е /п 
вы неравенства 

$ 

- 5ы/> //, ( фиг ^ (і) 6 I } (І)с/Ѵ-(І)<, сир //. ГЛ/ »аг Ц (і) . 

ГЬѴ г Ь.Ѵ [і оЛІ 
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Просуммировав последние по < , получим: 

т т $ 


- п 


5 Л,// Чг» , ‘ & * 


Пв .%1 


4 іі вир I/ - (і)\ ѵаг V; ('»). 

‘-Ъ.Ѵ П-М 

Заменим $ир ({, ( і)І на т ?.? Зш> !*.(*)! 

п/.ѵ ;•& ».кѵ н 

(неравенство о. втого лишь усилится). Учитывая, что в силу 
ограничения ѵ(і)г.ѵ(-Іа) П. ѵаг ѵ. (і) і, а . 

, Г. .11 ‘ 


имеем: 


-а 


^ »„« у 

и» ‘ѵ 

'»** Г«і /ЛГЙ, ‘й/^ 0Лг,г ' и 

•в 


«*/"«• “Р.ІАЮІ- 

‘*і,т Гі„Ѵ в ] 

Отсюда следует, что 

лі г% 

п/ И [ (і)гіѵ-(і)*-а тац зар І.(і} 

ѴшѵО/а) **' } 1 іаГ ” (І,М 


(Іо) 


Пусть экстремум в правой части неравенства (10) дости- 
гается при./-/* и на послѳдоватѳ явности і р — / * , т.ѳ. 

Сіт //,(•'/ *Ѵ" 5і - а //»('••''/. Положим е ({) * О 
р-*~ ‘ ’ 

при <>/*, Ч-, л (**»-« відпЦ.,(і 1 ,)І7С(1-і р '). 

Тогда С -«//,. ({,)}. • 

Переходя в последнем равенстве к пределу, убеждаемся, что 
л* Л 

Л/п Г ~ [ } (1)сІг ~-а$ир //.,('/)/. 

а—*"// 
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(II) 


Из (II) следует, что с участвует допустимая последовательность 
Ѵр(і) , реализующая в (ІО) точное равенство, а это и до- 
казывает справедливость леммы. 

ёШМВРир- Если г па$' 5ир достигает— 

ся на счетном мнокѳотвб точек ('*' ,(*• ,) * т °. как видно из 
доказательстве леммы, экстремум функционала достигается на мно- 
жестве отупенчатіи функций вида: 


а Л ,^,Егс„...гг^М,о, 


если при данном } лшт вир Ц м ( ОІ не достигается. 

Н 9 мИ /"4 іі ] ' . м 

В случае, когда ѵах §ир )і ( {) I достигается на 

' *•<> [ і 0 ,\> ,] 1 " 

множествах вида (а* , [ а ^* , ]) і то оптимальная 

функция Ѵ я (І) может оодѳ ржать и непрерывную составляющую. 

Из (8), (9) и лзммы I следует 

Теорема I . Для оптимального результата задачи I справед- 
ливо представление 


Ія/ г(и , ум) 
ѵ(‘)*Ѵ0/а) 


тах п(ы,Я(-і), 

ЛІ-)Ч < 


где 

* Іл/ ір(ы 1 /і(>) 1 ѵ ( ш Ь х 
Ф ѵО)сѵО/сі) 

- Г !'[л Т (*)ст«,*.)* в - р т (л«1>]<ью- 

** , ?' г (І 

- я ямѵ юпд I <И г I) ф)х (Ф, г)В(ё)\ йл (01], 

через ^ обозначена I -я координата вектора-строки 

Теорема 2 . для оптимальности управления ѵ в (і)е II (-/а) 
необходимо , чтобы »• 


тіп ѵ(*,у1 в ,ѵ)‘г4/(*,Л 0 ,ѵ в ), (із) 

ѵ(.)еѵ(./а) 
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где А„(і) - вектор-функция , доставляющая вкотремум з (12). 

Справ' тливость теоремы следует из свойств седловой точки А 0 , 
ѵ а • 

Теорема 3 . Для оптимальности управления ѵ 0 (-)я ѵ(-/а), 
удовлетворяющего условию (13) достаточно, чтобы вектор-функ- 
пия ограниченной вариации, доставляющая экстремум в (12), удов- 
летворяла условиям Л Ад (4) Ц 1 и на траектории х 0 (і) , 

соответствующей управлению ѵ 0 (і) , ы. почти всюду выполня- 
лось равенство 

?(са)х в ц),мц))-у1 0 г (*)с(()х 0 (*)~р м(4) (я/і)). (и) 


Доказательство . Предположим противное, т.е, существует 
управление ѵ(і) и соответствующая траектория системы (I) 
х(і) , для которой выполняется неравенство • 

Г К Ю) < Х(ы,ѵ„(-)). (15) 

Согласно (14), теорему 2 и формуле (2): 

гк ѵ в (-)) -/ р (с(і)х в а), т)сы(і ) = 

а *• 

- / О* ?Н(І) (*• 
і в і А 

* I Х(і,т)ІШч(*)]4ы({)* I [ЖМіЩіМ*'- 

*, і 4 

“ Рлм&М+Л'ММ) I * (і,т)ВШѵ(т)]сІ* (і) « 

А 

* / (л (ОФтМ- Рм(і) (А, (*>)] ы* (О Я 
* 

•, I ?(С(і)х(і),М ($<*«({). 

Сравнивая последнее неравенство с предположением (15) , 
убеждаемая в справедливости теоремы. 

Достаточные условия оптимальности в задачах импульсного 
управления рассматривались, например в [2, 4] , но в этих ра- 
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йотах ли о о функционал был другого вида, либо рассматривались 

другие ограничения на управление. 

% 
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УДК 519.6 


А. Г. Чѳ вдов (Уральский по- 
литехнический институт) 


РАСШИРЕНИЯ И ДВОЙСТВЕННЫЕ КОНСТРУКЦИИ Б ЗАДАЧАХ 
АСИМПТОТИЧЕСКОЙ ОПТИМИЗАЦИИ 

I. Введение . Рассматриваются задачи математического про- 
граммирования (МП) [I, 2] и оптимального управления (ОУ) [3, 

4) . Используются двойственные конструкции [б, б] ; формализа- 
ция приближенных решений соответствует [7] . Для одной задачи 
ОУ с фаэовыми ограничениями (ФО) получено выражение функции 
Веллмана. 

Пусть Е - непуотоѳ множество, & ( Е) - множество всех 

веще отвѳннозначных ограниченных функций на Е , з.е.В(Е), 
ф ф2сВ(П- ' 

Задаче МП 

5 д (X) — іп{ , хеЕ , со(х) «г<?. (и) е Я) ( 1 . 1 ) 

сопоставим семейство ее релаксаций 

<о(х) < е , (шеК), (1.2) 

где ее]0 , «= [ и М - конечное подмножѳотво Я . Подоб- 
но [7] вводим приближенные решения - направленности [8] в Е , 
соблюдающие каждое ограничение (1.2), начиная с некоторого 
момента [8, о.9б] , и обладающие свойством оходамооти соответ- 
ствующей направленности значений 5 " . Пусть X, - семей- 

ство всех множеств (х : X е. Е, (Фобе. К : со(х) < е)} ,гдѳ 

ее ]0 , С и И - конечное подмножество Я . В усло- 
виях Совместности релаксаций (1.2) имеем а качѳотве Х і не- 
четкое ограничение (НО) [7] і асимптотическое значение задачи 
(І.І) определяем, как (X, - тіп) [в^] [7, о. 19] . Смешива- 

нием (I . I ) назовем задачу: 
т 

У~ <л. 5 0 ( х ; ) — »■ Іл/ , пг - натуральное число, 
л • • * > * О у ^ <*і ■ ш 1 і 

х,еЕ , ..., х т е Е, ы. со (х^ * 0 (о) е Я ) . (1,3) 
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Задачи (І.І), (І.Э) некорректны [9, 10] , так что иопольвуем 
оерию релаксаций 
Аі 

П ы 1 сд(х { \ с с (и>*.К) , (і.4) 

где ее]0,»[ иА / - конечное подмножество 5? . Не ос- 
нове (1.4) также можно определить НО 3? 2 Г? Іи асимптотичес- 
кое вначение задачи (1.3) [7, с. 19] . Нсхоящѳние асимптотичес- 
ких значений составляет основную цѳль работы. 

^. Символика . Иопэльзуѳм кванторы, связки, символы^/ (по 
определению) и А (равно по определению)} через ф обоз- 
начаем пуотоѳ множество. Если ХІ и У - непустые множест- 
ва, тс черев У ѵ обозначим множество всех функций, действую- 
щих из V в V . Пусть Л' А [1\2 \ ... ] , Р - чи сда- 

вал прямая; 

Ѵѵ « /Ѵ\ й* А ( і : і е Ж, і в к] , Я к й р’'* , 
положи тельный конус . 

Множество Е наделяем полувлгеброй (см. [12, о.4б] , 
[іЗ, О.Э7, 39] ) X его подмножеств, для которой (обозначе- 
ния [13, о.60] , [14, с. 75] ) « г е 3(Е,Х) , ЯсВ(Е 

В (Е , X) - множество в„ѳх равномерных пределов ступенчатых 
отображений ( Е , X) в Р ; Е - полуалгобрв X , 
удовлетворяющая упомянутому условию на (е в , й?) , непременно 

существует. Пусть {аЫс()\ [X] есть о'е/ множество 
всех конечно-аддитивных мер на «Г [і4, о. 74] ; ом. также 
[15, 16] . Через Р д (X) обозначим выпуклую оболочку- [і7,о.ІЗ] 
множества всех мер Дирак;. 5, с (аХа)+ [Х ] , X в Е . 

Элементы Р 9 (X) формализуют эффект конечных взвесей в 
(1.3), (1.4). Пусть Р (X ) есть множество всех 

уи & ( ас/сО + рГ] , таких, что уи(Е)~1; 

Т(хи{уи^^Р(Х\(ѴішХ: (уи( і ) - О) V (фи (О -/])]. 

Через А(Х) обозначим множество всех мер л/-Ѵ , когда 
и V пробегают ( асісі )* [X] ; 4 (X) с нормой- 
вариацией изометрически изоморфно Е* (Е , X) (топологичес- 
кое сопряженное В ( Е , X) , так что наделяем А(Х) 

* -слабой топологией [15, гл.У] ; [16, с. 26 ] Г # (X) . При 
этом 
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Р(*)- р в (х) , Т'х)-{6 Х :хеЕ] , 

здесь черта сверху означает замыкание в (Л (X 1 ), т ѣ (Х)). 

3. Расширенія . Если (1.2) совместно при любом выборе 
€ е ]0, о° [ и конечного множества К, И с Я , то Я 

именуем слабо созместнім множеством; в этом случае % 
есть базис фильтра Е и асимптотическое значение задачи 
(І.І) определяется, как в I. Через 5 обозначим множество 
всех уи е Т ( X*) таких, что '■ / <о(т)уи (с/х) О ; 

интеграл понимается в соответствии с 14, с. 75 . 

Теорема 3, . . Множество 5? слабо совместно тогда и толь- 
ко тогда, когда 5 Ф ф . 

Теорема 3.2 . Пусть 5 Ф <р . Тогда (X, - ллл)[5„] 
является значением задачи 

Г $ 0 (х)уи (с/х) тіп (3. 1 ) 

Задача (3.1) разрешима (минимум достигается). 

Множество 5? назовем слабо совместным при смешивании, 
если (1.4) совместно при всяком выборе се ]0, <*> [ и конеч- 
ного множества А 7 , А' с: 5? . Через 5 обозначь ч 

множество всех мер уг/ е Р (X) таких, чтоѴ<*)е 52 : 

-$а)/х)уі і(с/х) 0. 

* Теорема 3.3 . Множество 5? слабо совместно при смеши- 
вании тогда и только тогда, когда § ^ ф . 

Пусть X г - семейство всех множеств: 

(у/ :// с р д (X ) , (Ѵй>е К . / а)(і е)/і (ах) < е)} , 

ж 

где В е ] О, * о[ и А 7 - конечно э подмножество Я .Ес- 

ли Я слабо совместно при смешивании, то Х г есть базис 
фильтра Р 0 (X) . т.о. НО [7] . Через П 0 обозначим функцию 
на (X) такую, что Ь 0 О* 7 ) » $ с 0 (х)/і (с/х). 

Если ? Ф ф , то (Х г - тіп) [іі а 7 е Я [7] 

есть асимптотическое значение задачи (1.3). 

Теорема 3.4. Пусть § Ф ф . Тогда ( Т г - тіл ')[/і (ѣ ] 

является значением задачи 
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Теорема 3.5 . Пусть ф ф . Тогда 



(3.2) 


(х г -Ш)[к ,] « 


= вир вир вир 



ІП} ( 5; (X) + 

Хее 


(3.3) 


4. Аппроксимативная выпуклость . Задачу (І.І) назовем аппро- 
ксимативно выпуклой (а.в.), если при всяком выборе Еб]0,«[ 
и конечного множества И , И с Я . имеем Ѵы. е [о, /] } 

X е Е , уе Е , что 3?с Е : Ѵи (в 0 ) и К 

и(з) 4 («и(х)+(/-ы)и(у))*е. 

Слабая совместность 52 влечет слабую совместность при 
смешивании. 

Теорема 4.1 . Пусть 5? слабо совместно, задача 
(І.І) а.в., тогда 

(X, - (%г-ПШ){Ьо\- 

Итак, если $ + ф и задача (І.І) а.в., то (ЦС, - /тл)[ 5 „] 

определяется (3.3). В І> (Е) вводим вир -норму 
[15, с. 26] . Если 52 слабо совместно и вполне ограничено 
в (В(Е\ )•/) [15, с. 34] , то семейство X 3 всех мно- 

жеств (х: X е Е , (Ѵь) е 52 : о)(г) «с е)] , ее]0, »[ 

есть НО [?) и (Х } ~ /пі/і)[5 р '] еЯ . 

Теорема 4.2 . Пусть 5? слабо совместно и вполне огра- 
ни ченр в (В (О, Гі) • Тогда 


(Х 1 - /7*Ѵ7)[5,] = (Х г тіп)[в д ] ; 


(4.1) 


если задача (І.І) а.в. , то (4.1) определяется максимином в пра- 
вой части (3.3). 

5. Задача СУ . Рассмотрим систему 


х = А(і)х + 0 (і,и) , «ер. 


(5.1) 


Здесь а? - п -мерный фазовый вектор (пс Ж) ; А (•) и 
& , •) - непрерывные функции, первая из которых матрично- 

значна; Р - непустой конечномерный компакт. Функциониро- 
вание (5.1) рассматриваем на отрезке Ій [і р , тІ д ] 

({„<!?„); Т й [і 0 , т? д Г . Через Ц обозначим множество 
всех кусочно-чостояньых и непрерывных справа на Т функций 
[/ер г . Фиксируем начальные условия г (І д ') •= а* (заданный 
п -мерный вектор) ; V V е V. через ір ѵ обозначим 
движение (5.1), по поддонное V из (і ді Х а ) - Тогда 

% - непрерывная вектор-функция, определенная на I , 

вектор уз у ( О имеет размеокость п . Ьыбером К е /~Л 
в качестве размерности пространства геометрических координат; 
Пи ость Ые/ вектор первых к координат п -мѳр- 
яого вектора х . В М -мерном пространстве определим 
систему множеств /1^ , і е I ; полагаем при этом относительно 
ФО //(,■) , что /Ѵ і , /Ѵ і * ф , ограничены в совокупности, вы- 

пуклы и замкнуты. Зависимость Д/, от времени полагаем не- 
прерывной (на Г ) в метрике Хаусдорфа [18, с.І7і] ; семо 
к -мерное пространство наделяем евклидовой нормой. Введем 
также непустой выпуклый компакт в к -мерном про- 

странстве. 

. Если 5 -непустой компакт в к -мерном пространстве, 
то через р(- , $* ) обозначаем функцию евклидова расстоя- 
ния до 5 * . Множество Е отождествляем с 21 ',$„(11) 

определяем, как <?([у и (& 0 )]и ■> ; Я полагаем 

состоящим из всевозможна отображений 

Рассмотрим получившуюся задачу (І.І), 

5, (Ѵ)~іп{, ІІе У, [ір ѵ [О], е /V, (+зГ), 

предполагая слабую совместность 5? , которое к тому же 

вполне ограничено в В(Е) = В (У) , так что справедливо 

(4.1). Наконец, тахсая задача (І.І) а. в. , так что для. нахожде- 
ния (4.1) можно использовать (3.3). 

Через Ф(* ? •) обозначим фундаментальную матрицу-функцию 
(решений) системы х = А (і) х • Определим систему мно- 
жеств Г % . і е 1 , полагая й У ^ при і с Т 
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И й N . Через И обозначим замкнутый единичный 

евклидов шар к -мерного л ростра нства с центром в начале ко- 
ординат. Если т е Л~ , полагаем, что Т т есть с/е { 

множество всех "наборов" в гр, 6 Т т таких, что 

V і 6 /,/77 , і < т имеет место + , . Пусть 

Ѵт е /Г: О^т = {о] ІІ Т^гп 
определяем условиями 


; если 6 Т т ,то 

I 


.(К Л *.)Ь (Цт К» : І- А і;)& (С, Л *.)• 

Аналогичным образом, Ѵт е А", (7Д- е Гот е пола- 

гаѳм^/). 6 у— у с таким, что ^ 4 <, при 

іе 1,т и й / . Наконец, если { & I и Л есть 

к -мерный вектор, то через р д (/]) (к) обозначаем 
значение опорной функции в точке к , т.е. максималь- 

ное (на ) значение скалярного произведения к' р , к] ; 
штрих здесь и ниже используем при обозначении скалярного про- 
изведения в качестве указателя транспонирования. 

Теорема 5.1 . Задача ОУ с ФО /Ѵ ( , (е Г , и критерием 

качества $ в (V) , I/ е V , обладает асимптотическим значе- 
нием: 

( Х і ~ ПМ)[ 5 о] - (% і -лм)[ 3 в] ■ 


$ир 

те/г 


$ир 

(*іЪе Гт еГ т 


зир 

ЮішІ?Ь е »а 


тах 

лес”*' 


&(ал<тс 


“ ?о (5.2) 

Равенство (5.1) устанавливается методами [19] , идейной ос- 
новой которых является общий принцип двойственности [к] . По- 
дробные обоснования всех утверждений настоящей статьи см. в 
[яо] . 
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У^СК 5ІЭ.854.2 Д.В.Зудихин (Институт 

математики и механики 
УрО АН СССР) 


ОНга-ИЗАШЯ дышл некого ПРОЦЕССА назначения 
I. Введение 

В статье рассматривается динамический процесс назначения 
заданий (работ, операций) некоторой совокупности исполнителей. 
Моменты начала выполнения задания предполагаются заранее за- 
фиксированными, а продолжительность выполнения неизвестна. Изу- 
чение таких задач представляет интерес по двум причинам: во-пер- 
вых, в связи с наличием широкого круга реальных ситуаций, где 
при распределении заданий важен именно момент начала их выпол- 
нения при неизвестной априори продолжительности, во-вторых, 
из-за недостаточной освещенности в литературе подобных динами- 
ческих процессов. В качестве примера рассмотрим следующую зада- 
чу. 

Задача диспетчера . Диспетчер, организующий работу И ис- 
полнитѳлзй, распределяет заказы на выполнение работ по обслужи- 
ванию М клиентов в желательное для каждого из них время. 
Каждый исполнитель должен обслужить не менее N клиентов 
(к/Ѵ ■» М) . Стоимость перемещения от одного клиента к друго- 
му может быть либо известна априори, если оговорены продолжи- 
тельность работы л вид , перемещения, либо неизвестна и опреде- 
лится например, заданными возле каждого клиента видами транс- 
порта и оценкой времени на выполнение заказа, проводимой по 
прибытию к клиенту. 

Диспетчеру требуется организовать оптимальное, с точки зре- 
ния заданного критерия качества, назначение заданий исполните- 
лям , если качество процесса оценивается взвешенной сушоіі вре- 
менных и материальных затрат на перемещение. 

Следующий пример является частным случаем приведенной выше 
задачи. 
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» а дача инкассатора . Бригада из К инкассаторов обслужи- 
вает М торговых точек (далее просто точек) по А/ штук 
каждый ( КА/ = М ) . Оговорено время с определенным для 

каждой точки периодом изъятия средств (денег) инкассатором в 
точке. Предполагаются следующие условия: 

- время накопления средотв в каждой точке существенно превос- 
ходит наибольший период времени обслуживания; 

- каждый инкассатор работает в условиях полной информации о 
действиях членов бригады; 

- апріори известен алгоритм вычисления затрат перехода от точ- 
ки к точке. 

Требуется организовать распределение торговых точек между 
инкассаторами: 

- обеспечивающее минимальные энергетические, затраты; 

- удовлетворяющее заданному ограничению на затраты; 

- обеспечивающее обслуживание точек за минимальное время. 

Прежде чем перейти к формальной постановке задачи, в -делим 
основные моменты, отличающие рассмотренные примеры от извест- 
ных в литературе [і, 2] задач о назначении, динамических пото- 
ках в сетях и теории расписаний. Для этого сопоставим задаче 
линейный ориентированный граф С -* (X , Л) , состоящий 

из совокупности вершин х в X и множества пар Л элемен- 
тов х е X , у еУ : ( х , у') е Л , называемых ребрами. 
Вершину в нашем случае отождествим с парой (р , Т р ) ~ х 
где параметр р определяет позицию точки (клиента), а Т р - 
момент начала обслуживания. Ориентация ребра задается момента- 
ми Ту , Тр , приписанными соответствующим вершинам у= (д, Ту) 
и X : при Т р <■ Т ч ребро ( х , у) ориентировано от г в 
у .Тогда упомянутое отличие состоит в следующем: 

- неизвестна заранее энергетическая стоимость ребрэ графа и оп- 
ределяется лишь в момент Т р , прихода в соответствующую на- 
чальную вершину X ; 

- каждой вершине соответствует заданный момент начала выпол- 
нения задания; 

- время выполнения задания априори неизвестно и определяется 
в момент Т р прихода в соответствующую начальную вершину х 
( в частности, может быть равно нулю); 

- позиция р может повторяться во времени в паре с соответ- 
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Сівувдим моментом Тр (і) : (р, . . . , (р, Т р <п} ) в случае, 

если ранее не произведено обслуживания. 

Приведенный отличия не позволяют использовать известные 
методы для решения задачи диспетчера. Однако в частном случав 
(задача инкассатора) при некоторых дополнительных условиях 
можно свести поставленную задачу к рассмотренной в теории рас- 
писаний задаче минимизации максимального штрафа за обслуживание 
п требований [г] т приборами, когда допускается преры- 
вание в обслуживании любого требования. При этом задано отно- 
шение строгого порядка, определяющее возможную последователь- 
ность обслуживания требований. 

В работе предложена численная процедура минимизации взве- 
шенной суммы затрат на выполнение задания. Алгоритм реализован 
в виде программы на ЭІМ. 

2. Постановка задачи . 

Рассмотрим абстрактное пространство позиций Р . Пусть 
на этом пространстве задано М многозначных отображений 
(заданий) 

А і : Р -Р , (I) 

где Р - 2 й , і е 1 , М . Затрата на выполнение задания 

Аі определим неотрицательной функцией 

В; : Р . ( 2 ) 

Назовем далее планом выполнение N заданий из М возмож- 
ных. Тогда процесс выполнения плана составит последовательность: 

1) , Ро еР 

_ Р-, > Р„ : 

2) Ы г е Г, м /(<*,) 

м) «*, е йм/(м'о К] .... и {*,_,}) 

Рл м е (р с * 

осуществляемую выбором номера о і { е /,М (Се /,/Ѵ) и пози- 
ціи лростргчства Р , на котором введен линейный предо о рядок: 
ДЛЯ ВСЯКИХ р И 9 либо р >9» либо 9 > р . Позиция р д 
в данном случае зафик сиро вана, а выбор элементов р ы . множест- 
ва (р^ У , I 6 /, N , подчинен условию 
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Ри іч < < Рт > 

где р т 6 Р - выбранная апріори конечная позиция. Случай 
р^>р т ( іе й # ) означает невыполнение плана. 

Таким образом, выбранная последовательность Р а1 ^(р <л ) 

определит величину затрат прі выполнении плана 

Ъ-Ц в. (р л , р ы ). О) 


Здесь символом ы, обозначена выбранная последовательность 
<* = (<*, , ■ • • , <**) ■ 

Задача I . Определить пару (ы.°, р ы ,) , минимизирующую вели- 
чину суммарных затрат на выполнение плана 


в °- л*а тіп гг Ь л> (р и . , , 
<* Р, и * •-> 


■Я-,)- (4) 

Задача 2 . Определить пару ('сі*, р^*) , обеспечивающую 


выполнение неравенства 
м 


БѴА*-, »А0 4А ’ 


(5) 


где & * О заданная константа. 

Задача 3 , Среди всех пар (V* , , решающих задачу 2, 

определить пару, параметр /> * которой минимален (о, * = 

✓ Л г л 

= (Ры* » • • • > А *;)) • 

В данной работе остановимся на исследовании задачи I , пос- 
тавленной для одного исполнителя при /V - ,Ѵ , предполагая 
производить дальнейшее обслуживание оставшихся клиентов испол- 
нителями последовательно. Предложенные далее алгоритмы решения 
зѳдачи I основываются на следующих понятиях и определениях. 

По смыслу задачи совокупность позиций 


Р ( Р- (р:р€А;(р д ), р <р т ) _ (6) 


при заданных р 0 и р т конечна. Перенумеруем их произвольным 
образом: 


Р 




К 6 /, Л’у , 


- 73 - 


Упорядочим далее множество позиций 


Г л 




/>Г 


><п 


ке1,к; ■, уе /,А} 


согласно заданному на Р прѳдпорядку. На полученном мно- 
жестве упорядоченных по возрастанию позиций Р^“[р ( '■ /,і] 

определим матрицу ^ задача назначения (матрицу графа С ) 
следующим образом. Элемент матрицы у ( и */Ѵ) зададим со- 
отношением: 


Д“ 


1 , Зх : Хе 1, к. 5 Р® - Р. ; 
О , V х : х е и~і<; »=»• Р® -*• Р , 


ІЩѲ < е /,Д ; /, А. 

Определение I . Ненулевые элементы /3- (іе /, и ;уе /,А) 
матрицы у назовем вершинами. Вершины р к( и р. . назовем 
несовпадающими , если і * у . Любую последовательность из А 
несовпадающих вершин 


А,«, » • •• ' 

назовем /V -цепью, есл и соответствующая последовательность 
позиций Ру ( (і е /, А) образует монотонную (в смысле 
введенного упорядочения) последовательность множества Р с : 


р. > о > 
V, 


> Р,- 

•я 


Таким образом, в задаче диспетчера всякой А -цепи бу- 
дет соответствовать возможный вариант обслуживания клиентов 
и наоборот. 

Определение 2 . Весом А -цепи , . . . , н 
назовем величину 


і1в- 5 вг л. < 7 > 

Следовательно, решение задачи I состоит в нахождении У - 
цепи минимального веса. Для численного нахождения А -цепи 
минимального веса предлагаются следующие алгоритмы. 

3. Алгоритм I . Рассмотрим А/ матриц ^ (хе Г, N !), 
полученных из матрицы У а» <р перестановкой столбцоч. 
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Рис.І. Блок-схема основной программы 


Определен:! о 3. Исходящей N -цепью матрицы (Ѵе /, /V') 
назовем всякую N -цепь Л- тЛу , А* ,...) при условии' 
г * * » </ * « для смі иьас элементов. 
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А лгор итм I состоит б определении для всякой матрице у 
(к ® ЛЛ'/) набора нисходящих /У -цепей, каждая из которых 
задается набором монотонно возрастающих чисел е 
/е ) , таких, что ^в іг ф 0 для всякого і е /,/Ѵ . 

Число таких наборов іуш матрицы Ц> к обозначим 1(.ч) . Тог- 
да величина X = II С (к") будет совпадать с числом всех 
возможных вариантов 'обслуживания клиентов. Выпишем эти вари- 
анты: 


К "Л--. К$) ■ *? е О, * , . • • , 1) \ К т < « X) ; 

А**® * 0 ; У в Л?ул) ; ѴУ/ } 

или, посла упорядочивания произвольна* способом: 

, Л»\ 




ѵ ) : і'е /,* }. 


Решение задачи I определится числом і* е / } Х , соот- 

ветствующим N -цепи минимального веса. Отметим некоторые 
особенности алгоритма I и приведем блок-схему программы (рио.І). 
Предложенная методика позволяет: 

- априори отсекать вершины физически нереализуемых вариантов 
обслуживания; 

- исключить по ходу дела вершину, не удовлетворяющую условиям 
или ограничениям; 

- формировать и запоминать вариант обслуживания \ клиентов 
с помощью одного натурального числа. 

К недостаткам алгоритма следует отнести возрастание объе- 
ма вычислений при больших N , Ь связи с этим рекомендуется 
использовать алгоритм I для N * 5 . Большее число клиентов 

можно обслужить последовательно, выделяя по 3-4 клиенте по не- 
которому приоритету. Основанием такой рекомендации служит 
быстрота счета в контрольных примерах оптимальных вариантов 
обслуживания при >Ѵ< 5. 

4. Алгоштм 2 . 

При создании следующего алгоритме сделана попытка для 
N Ь 5, сохранив достоинства алгоритма I, исключить некоторые 
недостатки, например, большой перебор возможных вариантов, а 
также многократный пересчет веса отрезков N -цепи, повто- 
ряющихся при обслуживании. 
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Для этого предполагается использовать память ЭВМ: во-нер- 
енх, запоминать неудачные (в смысле суммарных затрат) отрезки 
N -цепи о тем, чтобы исключить из рассмотрения /V -цепи с 
такими отрезками; во-вторых, запоминать стоимости переходов 
в позиции р е , ^ е /,/Ѵ , а также сами позиции. Это по- 

зволит снизить размерность задачи на единицу. 



Рио. 2. Блок Б 


Основные изменения внесены в БЛОК А блок-схемы (рис.І) и 
состоят в следующем: 

- вес N -цепи считается по мере ее формирования; 

- при невозможности перехода из вершин /<■ в вершину Н і + І 
запоминается отрезок N -цепи К, , # а , . . . , 

- создан блок апріорной отсечки вариантов, реализовать кото- 
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рые невозможно либо по определению матрицы (не /,/Ѵ!) 
либо по другим соображениям. 

Приведем БЛОК Б, который в алгоритме 2 заменит БЛОК А 
блок-схемы программы (рис. 2). 

Алгоритм 2 реализован в программе ОМА/? . Использованы 
стандартные комбинаторные программы РЕП М11 , 50ПТ2Ѵ. 
Параметры ввода: 

МА - число клиентов; 

ПО - максимальное число вершин матрицы У ; 
в - ограничение на затраты по обслуживанию; 

_ позиции клиентов; 

ТО - начальный момент обслуживания; 

Г - конечный момент обслуживания. 

Параметры вывода: 

ІА (А/) - номера столбцов матрицы У (Л/ е. /, МА) ; 

У5 - суммарные затраты на обслуживание; 

М1,...,ММА - оптимальный вариант обслуживания. 

1. Твнаѳв В.С. , Гордон В.С., Шафранский Я.М. Теория рас- 
писаний. Одностадийные -истемы. М. , І9Ь4. 382 с. 

2. Форд Л. Р. , Фвлкѳрсон Д.Р. Потоки в сетях. М. , 1966. 
276 с. 


УДК 517.977 В.П.Серов (Уральский политех- 

нический институт) 

К ЗАДАЧЕ ПОСЛЕДОВАТЕЛЬНОЙ ОПТШПЗАЦШ 

системы отклонений 

Для собственно линейной системы рассматривается задача 
программного управления на минимум системы рассогласований. 
Задача на экстромум в функциональном пространстве сеодптся 
,к конечномерным экстремальным процедурам на ЭВМ. ' Исполь- 
зуемые конструкции математического программирования следуют 
подходу [Т], предполагающему совместное исследование пары экс- 
тремальных задач, находящихся в двойственности. Работа отно- 
сится к кругу задач, рассмотренных в [2, 3] . 


1. Рассмотрим управляемую систему 

х-’А(4)х + 0(і,и) + В(і)Пс4 і ^(і-і е ),х(і в )^х в , 

ХеР п , иеР, Р*, і [Ь д ,$] . 

Здесь х - фазовый вектор; Р - ограниченное замкнутое 
множество из Р с ; 6(4) - дельта-функция Дирака; п>п - мат- 
рица А(4 ) и лкр -матрица 3(4) непрерывны; векгор-функ- 
пия &(•,■): Р *■ Р - 1 - Р п непрерывна по і, и . Для оп- 
ределенности считаем х(і ) непрерывной слева в точках 

. Назовем множеством всех допустимых программных 
управлений ІІ множеству всех измеримых (по Лебегу) функ- 
ций Ь ■ Г— Р _ 

Моменты і ь- подачи импульсов /, я) пронумеруем в 

порядке их возрастания 

і 0 6 І 1 6 і л < .. . 6 і м 4 1?; , ( 1 . 2 ) 

суммарный ресурс импульсного управления равен уи > О : 

ИК/4/ (1.3) 

5 »/ 

( І'І - евклидова норма). Движение х(і ) , іеТ системы (І.І) 
из начальной позиции ( і д , х^) определяется формулой Коши 

і 

где К[і,т] - фундаментальная матрица-функция решений систе- 
мы х - А (Ох. 

Пусть задана система непустых выпуклых компактов /і ; (і е.[^т) 
из Р* и, кроме того, 

1 Л {йеР Щ суи , оі. е Р р , і е 1~к ), 

Уа{ге^ да |^-Г <; 6Г^ ... 4Г Я 4 I?}. 

2. В этом пункте рассмотрим задачу на минимум взвешенной 
суммы расстояний от траектории системы до множеств М- , т. е. 
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качество управления будем оценивать функционалом 

^ а, р) - г «; а (р. х(г) ,м.х (2.і) 

где V; (іе1,т) - эаданиыѳ /ип -матрицы; <і 1 >С(іе %т)~ 
ва данные числа; а'- 4 ; символом ‘ 4(у,М) обоз- 
начено евклидово расстояние от точки у до множества М 
Пусть У - множество векторов р'- и,,..., , координаты 

которых удовлетворяют условию (1.2). 

Задача 2,1 . При Зиксированных Ге Г , ре У вычислить 

е (г, р) - /я/ (г (Г, Ѵ,с/,р)І([/,сОс 1/*2]. 

Пусть і" 1 - множество всех векторов-строк А та- 

ких, что 

Л «(Л Л (л»)), |Л(ОІ < /, лѴЛе#*, і в . 

С помощью соотношений 

тах и'у с -р <?/*<)) - 

ё «іілщ-еМ'Ѵ] 

получаем для величины (2.1) вырежѳетѳ 


Дг, о; <*, р) - тах ,Л, г, 1Г } </, р)| Л « І т ] , 


где 


са^в.л.г, Ѵ,4,$-Г*і [АЩХ($ -р(Л(П \ *)] , (2. 2) 

Щ)-Чь *•]*.♦ ( +Пх[ъ,* в ]в(Ца, , 
л/(4-(8ш?ряи { <ь). 

Аналогично [2] имеет место равенство 

ГМ,р)-< Л /[ Г (г ,^,р)|(/ е 1/} - 


(2.3) 
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- тот (Н(«с,Л,г)+І1ы { А (00; ЕГ X [г іЛ В(іАЫ. ] , (2.4) 

Лбі* 4 і‘І іе#(/) 1 1 


тле 

И(сі, А, т)~П а,- [Л (Щ X [г ( , і а ] т в -р(Л <7)| А/-)] ♦ 

*,? і‘ (ё4 л щ*ь,міс<’щ.4</і ■ 

С Л' 

Функция б фигурных скобках в (2.4) выпукло-вогнута по (СІ,А ), 
поэтому по тѳорѳмѳ о шіьлмаксѳ [4] : 


Е (г, р) - пчп [у (г , тах {/?и,Л,г) + 

т Ае.І." 1 ^ 

тіп ЕГ <*,• Л(іЩ ЕГ Х[г,ФЛ &(і 5 )<У, } . 

ае 2 ».» іел'О) } 

Пусть теперь К-т и , кроме того , моменты подачи импуль- 
сов удовлетворяют условиям і $ е /г 5 _, , г 5 ) , зеТ^т. 

Зафиксируем вектор се Г и определим множество (?'=(*,, ... 

■ * Ч ^7»)) 


К- = * і в , бе^т). 


Тогда 


е: лгг.-.^х -± х/г . 

$ШЛц) $ші 

(2.6) 

гл і т т 

Учитывая равенство Щ ЕГ = ЕГГ ЕГ б 1 ,', 

1»? г»/ *»/ і»« 

, получаем 

ЕГ <*;УЦП2). ЕГ X /г. , ]&и г )4{ = 

І»і $*< 

(2. 7) 

= Е II «і Л(сЩХ[ г, , і,]В(1,)е(, . 

$*1 1*5 


‘ Из (2.5)-(2. 7) вытекает справедливость, равенства 


&/г,я) * ядах (#і~*,Л,гѴ ; 

Леі'” 

(2.Ѳ) 

777 

I Г ^ лгоА х Гс„ *, ] а (Оі) • 



Перейдем теперь к нахождению оптимального управления (7° . 

При этом предполагаем, что по (7 из і/ в эадаче 
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2.1 достигается. В противном случае следует расширить класс 
допустимых управлений до обойденных управлений -мер [5, 6] ; 
в атом к ль „се управлений наше задача на минимум всегда имеет 
решение, и оптимальное обобщенное управление удовлетворяет 
соответствующей модификации принципа максимума Л. С. Понтрягина. 

Справедлива (при фиксированных ) следующая 

теорема. 

Теорема 2.1 . (принцип максимума). Если (1Г° с!")е Л* 2 
оптимальное решение в задаче 2.1, в А 0 - вектор, доставляю- 
щий максимум в (2.6), то 

а) для любого І7т почти вс оду на [г 4/ , г с ") выполня- 
ется соотношение 



Задача 2 , 2 . При фиксированном ге Т вычислить і(т) = 


шІл/{е(т } р) ІреѴ ; ). 


Задача 2.3 . Вычислить е°» [с(г)|ге Г ). 

Ниже, в п. 4, задача 2.3 нахождения оптимального вектора г* 
решена для случая так называемых "простых" движений. 

3. В этом пункте качество управления будем оценивать сте- 
пенью максимального отклонения траектории системы от множеств 
Мі . \ именно, будем минимизировать функционал 


тах сі ( Ъ. х(г г ) , Ц). 


іе 1,/п 


Имеет место равенство 



где о) задается соотношениями (2.2), (2. 3) ,^а 
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вычас- 


Пуоть т и е , г 4 ] , 5 & /, т . 

Задача 3,1 . При фиксированных г е 7', р е У^ 

лить 

е (г, р) - { рСг, II, <і, р)\(Ѵ, 4) е I/ * 2]. 

Задача 3.3 . Пря фиксированном Т е Т вычислить 

2 (У) « і/г/ [ 2 (V, р)І р е ѵ г } . 

Задача 3.3 . Вычислить 2“’« гл//"2^г))ге Г}. 

Справедливо следующее равенство (см. [?] ): 

2 С г, п\ - /пах тах Л , г) - 

(3.1) 

/таг | ДвС,АаЩХ(*і,#,]в^}. 

5ъ1уПІ *■- 

Кзк и в л. 2, предполагаем, что гл/ по 17 из 77 ь зада- 
че 3.1 достигается. 

Теорема 3 .1 (принцип максимума). Если (17\е(*)е 17*2- 
оптамальное решение в задаче ЗЛ, а ес* г Л* - векторы, до- 
ставляющие максимум в (3.1), то справедливо следующее утверж- 
дение: 

а) для любого зе1,т почти всюду на Г*Ѵ-/« */) вшол- 
няется соотношение: 

Т1ы‘Л а (1Щ Х[г : , і)0(і , а У О) - 

- лп-л с <Л7< )2>. X /г, , і]С (г ,и): 

б) 21 Пы'Л 0 (<ЩХ[г с , - 

= -л/ л*тх | II ы*Л*(ЩХ[Ъ , і'І&аХі- 
$6/, т і‘$ 

4. Рассмотрим решение задач ЗЛ - 2.3 для случая "простых" 
движений. Система (І.І) имеет вид: 

Х*и+ ВТІа! х(і^’Х 0 , (4.1) 
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и е Р, с( е 2 ; матрица В постоянна. Отметим, что для систе- 

мы (4.1) 5 - А импульс .'ложно подавать в любой момент времени 
е Г г *-» » Г» ) - вектор (2.3) при этом не меняется, так 

что задача 2.*, в данном случае выроэдаѳтся: г(т } = е.(г). 

Функция б(т') согласно (2.8) имеет вид с(т') = 

- тах[Я(г,Л)\Л в Г], цце 

XI <*; [Л(Щ х в -р6иО|А/;)] + 


"Ц , т . т 

+ Х_ (<7 -с,.,') лкя ( ЕІ <=<.- -Л Г02?- и) -/и тах ІГ а..К(Сф-Ь\. 

і т ! иеР'і.] ' $ефп ‘ 

Пусть 

б - {<5е *" | (5 ; і 47, і е , С 6- 4 г?- ^}. 

Тогда (хе# 1 *) 

Ѵ-, (Ѵѵ-'И “ Т іп &‘* я ( *-<•) *»> л*; • 

г еГ/«» ' ' 4 ИеД .«/рй ' 

Функция 5? ("г , УГ) линейна по г и вогнута по Л 
множества Т к і ~ выпуклые компакты, поэтому 

Е°*тіп тах Я(ТуЛ)^тах гтп Я (г, у|\ = 
г вг ЛаІ а Лес» теГ 

"лес* ( Е ^ (01>і *' ' * 


(4.2) 


+ тіп ( тіп П и ; Л(іЩ и)-уа тах Ц И <*■, Л(і)й ■ 31 ] . 

у'» V» итР і.; 4 ее/^т "і-і ‘ 11 

Пусть 


Г^{г*Г\е(т)-*°], Я>(г,Л)-е°]. 

Множества 7^ , непусты в силу эамкнутости Т , і т и не- 
прерывности 5? , • ) . 

Теорема 4.1 . Пуоть г°е Т ^ - оптимальное решение в задаче 

2.3; (1/° , а °) - оптимальное решение в задаче 2.1 (при 

С - г° ); Л°е I* . Тогда 
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а) для любого 56 и т почти всюду на {V*., , г') 
вшолняетоя соотношение 

П Ы;Л Ф (іЩѴ'(і) - Я*> (І. «іЛ 7«)і). «) ; 

4-5 Ыв/> «'•* 

б) вектор аС’ удовлетворяет условию 

И ЕГ <*,• Л7‘)Ц йаі" ‘•-'/і тах I ЕГ <*,- Л’(7)А 3 I; 

ІЧ І-і $еГ^п чі 1 

в) координаты вектора с ® удовлетворяют условиям г/ - г“_,+ 

<5, , 5 е , где г/ ■* г„ , а вектор ^ ° находит- 

ся из условия 6'(Л°)6 в ш/піп?(Л 0 )6 .где 6.(Лі*тіп(Ііы.Л(іф.и\ /7*. 

ЛеД і/еЯ 4 

Таким образом, схема решения заключается в следующем. Из 
(4.2) находим Л” , затем по Л° находим необходимые уо- 

ловия для V" и й° , далее находим вектор' 6(Л°) , потом 
по 8(Л°) находим вектор 6° и, наконец, по 6° находим 
г • . Как и в [2] , имеет место следующее условие выравнива- 

ния. 

Теорема 4.2 . Пусть существует невырожденный вектор г°& Г %1 
т. е. вектор, координаты которого удовлетворяют соотношениям 
Т 0 < г, с г г <: . . . < ѵ т . Тогда для любых А° е I. # 

ме , 5 е Ц~т справедливо равенство б я (А а )»б і (Л л ) 

5. Применим вышеизложенную теорию для решения прикладной 
задачи из области медицины. А именно, рассмотрим задачу выбо- 
ра оптимального режима дозирования лекарственного препарата 
Іл.п. ) в организме больного. Пусть с(і ) - концентрация л. п. 

в организме, который представлен п. -камерной фармакокинети- 
ческой моделью (см., например, [8 - 26] ). Производится непре- 
рывное и импульсное введение л.п. с целью создания его нужной 
| концентрации. Поскольку мы лекарство только вводим, то на ско- 
ростъ введения и разовые дозы л.п. сразу следует ограничение: 
они неотрицательны. Время выведения инъекции или таблетки из 
1 организма на несколько порядков выше времени их поступления в 
организм, поэтому для описания введения разовых доз л.п. в. урав- 
нении модели процесса можно применить <5 -функции. Кинетику 
всех процессов переноса и элиминации л.п. предаолах-аем линей- 
ной относительно с(і ) 


Скорость изменения ( д; ( і)) концентрации л.п. в і -й 

камера описывается уравнением (і е (,п ) 

С;«) — П <*;.(()€;(*) + П а^сЛО-к.^Ка)* 

/-Ѵ« ' *",*»' 1 (5.1) 

+ Я (0</;(()+Л { а)11 а я 6( і-і'), с : (Ь^~с я , 

где а у (і)^0,((і > ^')е х Г^п , і+/) - коэффициент скорости 

переноса л.п. из і- й камеры в у - ю; к- (і) іО (і>і 0 ,іе^п) 
коэффициент скорости элиминации (т.е. выведения л.п. из I - й 
камеры во внешнюю среду за пределы а -камерной модели) для 
I - й камеры; р.(Я*0, 8;(()>0({в%д)} с/^ д' ад ). 

Все скалярные функции , к 1 , д і , 8- непрерывны. Отме- 
тим. что в общам случае а + а^, . Таким образом, л.п. 

в I - ю камеру вводится непрерывно со скоростью о. 

$ разовыми дозами 8- ( I 5І в моменты времени . Отре- 
зок времени, в течение которого вводится л.п., фиксирован. 
Число инъекций щ также задано, а моменты инъекций 
і п . и величины доз могут выбирать- 

ся. Обозначим 

а и (0 - И а-, (і) > О (іе І^гі); 

4 Ш Ч* 1 

с(і)-(с,(і),...,с„(Я) , к'(і)- (к, к„ (і)) ; 

і'(0 - (і (*\- • • , &(*1 , е'(і) - (І, (О, ...,3 Й (О) ; 


Л(і)~ 



а ц(*> ••• 

\ 

<>»($ 

- а 2і (0 • • • 

а і л & 1 


' • 

-«*«) 1 


Тогда система (5.1) запишется в виде 

, ( 5 . 2 ) 


- 06 - 


где Е - единичная л «/7 -матрица; 6 (і)-с/іад {0 (і), ... 

* В(*)”С/іар (в, (і) , ... , 8 п (і )] . Убедимся в 

том, что модель сохраняет физическое свойство исходной ре- 
альной системы (организма), заключающееся в неотрицательнос- 
ти с(і ) . Известно, что условие Пц(і)>0 (і+/) явля- 
ется достаточным для неотрицательности решения х( 1 ) > Ь , 

Уравнения х(і) - А/(1)х(і) , х(і 0 )*0 (п г (0 - элементы 

лк л -матрицы Л' (О ). Отсюда с учетом неотрицательности 
координат векторов д(і) , в(і) , а п ) , с / 5 , с„ вытекает, 
что решения уравнения (5.2) неотрицательны при любом выборе 
с 9 ь О и допустимых управлений и(і), СІ $ . Вид матрицы 

Ек'({)) (в частности, число и расположение ее нуле- 
вых элементов) определяется, очевидно, структурой данной л -ка- 
мерной модели. 

Ограничения сверху и (і) е Р , Ые 2 .где 
Р-[иеВ"\Очі 1 . <е/ГЗ)'; 

г,{«,е""\й Гсі„ 

' 9*1 1 = / 


на скорость введения л.п. и на величины доз обусловлены сооб- 
ражениями клинической практики и имеют своей целью оградить 
организм от нежелательных побочных явлений. 

Ущерб от отклонения реальной концентрации с(і) от жела- 
. іельной будем изменять функционалом (см.п.З) 

%(і,и,4,п)~тах Цс(і-.)-л? ( / , 

і е 7 ,т 

где і д = & Т і 4 , ... і г т - заданные моменты времени; 

координаты векторов т і (і <= /, т) неотрицательны. Таким об- 
разом, желательный уровень с(і) задается величинами 
т п ...,т т в моменты времени Т/, ... , Г т . Вели, например, 

\ т. = О , то это означает желательность того, чтобы в момент 
времени г. л.п. вообще не было в организме, а если і - я. 
координата вектора т '■ равна нулю, то в момент Г. л.п. 

не должно быть в і -й камере. Имеем 


П\ п 


Г Гълтщ , -гг а і$ й , (5.3) 
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где а и . ("у е /, п ) - координаты вѳктора-строки а і% » 

-*іЛ(ПХ [Т±, (■;]. 

Из (5.3) следует равенство 


т т 

тіп ЕІ ЕГ А( ОХ/г., 1АВ(*Лс1, - А /ям 

<*«2 *->«-* 1 1 


ШОГа ). 


Функция (3.1)' имеет вид 

Ь(г,?,с,)-тах тах яе(а,Л ,т, п, с в ), (5.4) 

Л«і* о»«5 


где 


аг (<н , Л , г, ? , с/) - Е<*,-Л (о[* [г , („)с,-/л.] + 

V* 

вектор-строка ?у(^) > ?/*(?)) определяется ра- 

венством 

?,• (?) - Е <*,-Л(У)Х /г,,?] , уе^. 

Условия, которым удовлетворяет оптимальный режим лекарствен- 
ной терапии 11 °, й° запишутся согласно теореме 3. 1: 

в) для любого 5 е /^/л почти всюду на [г , вы- 

полняются соотношения (къіТп): 

9х(^и ( і)и * СО ш (*)/»>» [ 0 ] 9 , ( 0? Л се) , 

где С^) “ Е <*?Л'(І)Х[г. ,1] : 

6 ) вѳктор сі 9 удовлетворяет условию 


т т л Л 

Г Г а $ В(^ - 

л» >і *• 

где л* - ос/ЛѴОХ [г. 

Аналогичные соотношения можно получить и для случая, ког- 
да эффективность терапии оценивается величиной (см.п.2) 

ЕГ*; - т с \ . При этом функция (2. Б) имеет вид 

€ (Т } р, с д ) - ^/тюг ар С*,-Л ,г, 9 } Сд'). (5.5) 


А ЛѴЯ _ (б(і.)И о"..), 


Б (5.4), (5.5) величина С в введена в состав переменных 
функций значения задач 2.1 а 3.1, цоскольку на результат зада- 
чи установления нужней концентрации л.п. можно влиять и выбо- 
ром начальной концентрации с с (т.з. выбором С а (і « ) 

во всех камерах). В связи с этим возникает следующая задача: 
Пусть мы имеем возможность варьиѵювать Сд в некоторых пре- 
делах: с в е С , где С 9 - непустой вьлуклый компакт, 
содержащийся в неотрицательном ортантѳ пространства Я" , 
Задача 5Л . Найти (при фиксированных г, р ) вектор 
Сд .доставляющий равенство с' ( г, р , с ° ) «• гЫп {&(?, с д )] с 0 еС*}- 
Задача 3.2 . ;дйти (при фиксированных г, г/ ) вектор с/ , 
доставляющий равенство ё (с, і?, с 0 °) - гт.п {& (і, <Ѵ) | с » 45 )• 

Отметим, что функции С (Г, ?, і„) и ё (Г, Сд) выпуклы 
по с 0 и достигают своего минимума на компакте С’* . 

Рассмотрим решение задачи 5.1. пункция X вогнута по 
Л и линейна по С 0 , поэтому цо теореме о минимаксе 
получаем: 


тіп е (г, о, с а ) - пни тох ге(а т Л,т. п, с в ) “ 

с 0 *С* * с 0 « С Л с і, т 


- />юх(х, (л ,л , г, ?)|Л е І т ] , 

Х,(*,Л,г П«. ; Л(Ол>; «- 

У - ц л (Л К {(< ,(,]■ 

Отсюда находим условие, которому удовлетворяет С? 


(5.6) 

(5.7) 
1 

(5.8) 


уСд - тіп(уСд\Сд е с в } , (5.9 ) 

где вектор у /задается равенством (5.8) при Л* А* , а 
Л* - вектор, на котором достигается максимум в (5.6). 

Обратимся теперь к решению тадяпз 5.2. Фу„осі;, я сС (2.2), 
(2.3) для системы (5.2) линейна по переменной 1Г*(1Г,с/,с^)е. 1/*2*С* 
Можно показать справедливость равенства 


ё (г, о) г- тах тах эе, (<*, Л , г, /у) . (5.10) 

‘ леі т ые-5 

где X, (<*,Л, г, $>) задается (5.7). Оптимальна» начальная кон- 
центрация Сд* удовлетворяет условиям (6. 9), (5.8) при 
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оі‘ Ы.° , А « Л , ГДѲ аі° у Л - вѳкторн , нв ко- 

торых достигается максимум в (5.10). 
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УДК С2.50 


Г.И.Найфѳльд (Уральский по- 
литехнический институт) 


ПК'ЕКТИРОВАШЕ ЗАКОНОВ УП ЯВЛЕНИЯ САМОЛЕТОМ 
МЕТОДОМ МОДАЛЬНОГО СИНТЕЗА 

Введение 

При проектировании си стам управления полетом принято вы- 
делять задачи синтеза внутренних и внешних контуров управле- 
ния для продольного и бокового движения самолета. Критерии ка- 
чества внутренних контуров традиционно формируются в виде 
требований к переходным процессам, расположи кию корней или ви- 
ду частотных характеристик. Внешние контуры управления исполь- 
зуются для реализации программных траекторий по отдельным сос- 
тавляющим вектора состояния: по углу тангажа и высоте в про- 
дольном канале и по углам крена и рыскания в боковое. Качест- 
во внешних контуров оценивается по точности отработки програм- 
мных воздействий, отсутствию автоколебаний и г.д. Расчет зако- 
нов управления обычно осуществляется частотнши и корневыми 
методами, что объясняется простой физической интерпретацией 
зтих методов. 

В данной работе проектирование законов управления движе- 
нием самолета осуществляется методами модального синтеза, что 
позволяет получить единообразную форму закона управления в ви- 
де наращиваемой по сложности структуры [і] . Для облегчения 
изложения методика проектирования рассмотрена на численном при- 
мера. 

I. Математическая модель самолета 

Пространственное движение самолета для фиксированного ре- 
жима полета описывается системой уравнений [2] : 

оі-д), -ідр*и> х -0,667*<*/сО5/~0,Н*6 $ /соі/) +0,047]* С05 О . со5/‘)ссг р ; 

/3 - О х * 5ІПЛ +<Оу * СОМ * 0,047} 0)5 0* 5і/г/~- 0,23 */ - 0, 046 * б н і 

С0 3 ~-0, 09* и) г - 4,8* - 8,7* 6 в ; ® 

сЬ х --2,6»и) х - 0,25*07 у -38,0*/ 17,0 *б 3 - 7,о*$ н ; 
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со у - 0,075 *сО х - 0 , 27 * со у - Ч,Ч*р + 0,82*6, - 5,2*0, ; 

0 - СО у* ііпр -1-1 0;* С0${- ; 

(2) 

1 ш 0> х -ідО* (сОу * с05/г-и) г * зіл/г) ; 

Ц/ » ((Оу *СО$Х -СО г * $ілу)/ С05 О ; 

И - 207,2 *(со$ы.* со$р*со$ц/*со$'д+зіі*і*соі^*(со$ц/ *зіпі7* 

* С0$І"+5іП У*5 іП]г) + ІіЛ/ *(С05ІІ/*5ІЛ^*5іп^+ 5,ЛЦ * СОЗу'І). 

Чп сланные значения аэродинамических коэффициентов описывают 
самолет Р-С при Л/ = 0,6 , Н - 6000 Ом [3] . В полный вектор 
состояния входят: <* -угол атаки, ув -угол скольжения, (О х , 
(0у,(2>ь - угловые скорости, іУ - угол тангажа, )Г - угол 
крена, цг - угол рыскания и высота полета Н . 

Приводы руля высоты, элеронов и руля направления имеют 
ограничения го положению: 

-26,0 <6 Й < 6,0° ; | <5, | < 30,0° ; )б н | < 6,0°. 

Моделирование пространственного движения с помощью систе- 
мы уравнений (I) не учитывает изменение скорости полета и аэро 
динамических коэффициентов при изменении составляющих вектора 
состояния Модель органов управления не учитывает динамику 
приводов, ограничения скорости перемещения управляющих поверх- 
ностей и нелинейности типа "люфт" и "зона нечувствительности". 
Кроме того, считается, что все составляющие вектора состояния 
доступны для измерения и могут быть использованы в законе уп- 
равления. Эти упрощения не имеют принципиального характера с 
точки зрения применимости методики синтеза законов управления. 

2. Метод модального синтеза 

При проектировании систем управления самолетом традицион- 
но принято разделять продольное и боковое движение. Это разде- 
ление возможно, если перекрестные связи между продольным и бо- 
ковым движением сравнительно малы. С точки зрения методики 
синтеза законов управления различия продольного и бокового ка- 
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налов управления определяются размерностью вектора управле- 
ния и необходимость раззязки движений в боковом канале уп- 
равления. ' 

Рассмотрим объект управления в виде 


(3) 


х - А чх ■+ В*и , 


где X- п -мерный вектор состояния, и- /п-мѳрннй вектор 
управления. 

Управлониѳ и форлируется в виде линейной комбинации 
командных сигналов и вектора состояния 


( 4 ) 


І/г-Ьчх + С*ѵ. 


Матрица обратной связи рѳосчитываетсл с помощью соотноше- 
ний [3] : 


і «у/* Г ' 1 ; 

А*Г - ГѵЭ - В *// , 


( 5 ) 


где 0 - желаемая диагональная матрица собственных чисел 


замкнутой матрицы динамики З і , (і » , л) ; Г - матри- 


ца собственных векторов замкнутой системы; уи - матрица 
вспомогательных векторов. 

Задача линейного синтеза решается, если будут заданы мат- 
рицы 0 и Г. Однако определение этих матриц не являет- 
ся очевидным даже при наличии требований к переходным процес- 
сам или виду чзстотных .характеристик, которые обычно фигуриру- 
ют в техническом задании на проектирование, Фактически, тре- 
бования технического задания позволяют определить область на 
комплексной плоокости, в которой должны размещаться нули и 
полюса звена второго порядка, имеющего заданный вид переход- 
ных процессов. В реальной системе управления качество будет 
определяться не только обратной связью, но и степенью влияния 
нелинейностей исполнительных механизмов, ограничениями ресур- 
са управления и т.д. 

■'ак показано в [і] , матрица коэффициентов обратной связи 
для скалярного управления рассчитывается по формуле: 



_л[_ 

в 
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(С) 


п 

гді-і СІ (Л) -П (Л - б' ( ) - жэлаемый характеристический 

многочлен А ; <5,- желаемые собственные числа замкнутой 

системы управления; Н~' » (А* } . . . , Л+ ) т - матрица 
двойственных векторов; Я і - собственные числа матрица А . 
Очевидно, что минимальные значения коэффициентов обратной свя- 
зи будут получены при перемещении корней Я і матрицы А 
в точки з і которые должны находиться на границе желаемой 
области размещения и быть ближайшими к Я ; . 

Для сиотем, имеющих скалярное управление, линейная обрат- 
ная связь позволяет разместить собственные числа замкнутой 
системы управления, т. ѳ. удовлетворить треб отава я к динамике 
переходных процессов. Для векторного управления может быть ре- 
шена задача развязки движений по составляющим вектора состоя- 
ния. В этом случае должны быть заданы собственные числа и соб- 
ственные векторы замкнутой системы. Матраца. собственных век- 
торов может быть задана в виде 


Г 


а . 

. к 

0 . . 0 

0 . 

. 0 

* . . * 

0 . 

. 0 

0 . . 0 


ш 



Я (І г 


где у. - количество корней, определяющих движение по коор- 
динате х- ; ? - произвольные числа. 

Далее, из выражения 

(Яі*і-А) - {. (Ь) 

определяются //. и рассчитывается матрица обратной связи 
і. . Следует отметить, что. полная развязка свободных движе- 
ний возможна лишь при равенстве размерностей вектора состоя- 
ния и вектора управления. 

3. Внутренние ко,: уры управления 

Требования к качеству ручного управления формируются на 
основе анализа и формализации понятия "хороший самолет" в тер- 
минах теории автоматического управления [4] . Известны неоколь 
ко критериев, нормирующих отдельные параметры системы унравле- 
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кия, однако, принципиальным требованием является апериодичес- 
кий или слабоколебательный характер переходного процесса с 
заданный временем регулирования по наиболее важным для летчи- 
ка координатам. 3 продольном канале тагами координатами явля- 
ются перегрузка и угол атаки, а в боковом - угол скольжения 
и угловая скорость крена. 

Для расчета ручного управления используются уравнения: 

<*■ «й>, -0,667 *оі- 0,11* 6 й ; 

- 9,6 *:* -б, 7*6 В , 

с'і, * - 2.6 %<.д х -0,25* у - Х,0*ув - 17,0*6 } -7,0* 6 н ; 

«й * - 0,С75*сй х - 0,27*0, у -9,9*р + 0,62*6 3 -8,2*6 щ \ 

/ - 0,076* (Л х * 0,027* О у - 0,23*/і - 0,096*6, ■ 

Собствешшѳ значения матрицы динамики продольного кана.па: 

Я 13 *- -0,56 +/2,2. Переходный процесс разомкнутой системы 

по углу атаки имеет 
колебательный харак- 
тер о перерегулиро- 
ванием 405? и време- 
нем регулирования, 
превышающим 5 с. 

Предположим, 
что желаемый переход- 
ный процесс должен 
иметь перерегулирова- 
ние но оолеѳ 10% л 
время регулирования 
нѳ более I с. Эти 
требования определяют область желаемого расположения корней, 
представленную на рис Л, 

Ближайшие к собственным значениям разомкнутой матрицы ди- 
намики точки этой области определяют’ желаомыѳ корни замкнутой 
системы 5 2 2 = ± / 2,0 и закон управления 

5 & -0,568* <*. *■ 0,699 *иі і ■<- Х р , 
где Хр - сигнал датчика перемещения ручки летчика. 
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йгс.І. Желаемая область расположения 
корней системы ручного управления 
продольным каналом 


Разомкнутая матрица да на лига бокового движения имеет три 
собственных значения: Я , * -2,25, Я }1 = -0,44 + ^ 2,56. При 
наличии двух органов управления [5) линейная обратная связь 
позволяет разместить собственные числа замкнутой системы в 
желаемой области и разделить свободное движение но ц> Л от 
свободных движений по иіу и / . 

Предположим, что желаемая облаоть размещения корней опре- 
деляется соотношениями 

$е[С,5, /], 

которые определяют аперіодический процесс с Г = 0,2Б по а) х 
и слабоколебате "ъкьеі процессы с Т *= 0,38 , Ц -- С, 52 по&^Д 
Матраца собственных векторов замкнутой системы 




1 1 

0 

0 ^ 


Г я» 


/г (О 

ш 



\ 0 

I 

I / 

где ) 


- функции собственных значений 5, .. . 5 , 


Линейный звкон управления, рассчитанный согласно [5] 



-0,0496 0,257 2,26 


<*>х 


У, 


* 




-0,0077 -0,589 -0,061 


/» 


? 


где Х 3 


- сигнал датчика отклонения ручки летчика. 


4. Внешние контуры управления 


В отличие от ствтичѳских внутренних контуров внешние кон- 
туры управления проектируются астатическими, что позволяет 
реализовать отработку заданных значений и стабилизацию отдель- 
ных составляющих векторе состояния. При этом считается, что 
объект управления включает в себя ранее спроектированный кон- 
тур, корни которого не требуют коррекции. 

Типовая структура внешнего контура, включающая интегратор, 
нелинейные звенья, внутренний контур управления и обратные свя- 
зи внешнего контура, представлена на рис. 2. Величины ограничений 
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Рио. 2. Структурная схема внешнего контура управления 


определяются при численном моделиоовании процессов отработки 
больших входных воздействий. Следует отметить, что при таком 
подходе внешние контуры, относящиеся к одному каналу управле- 
ния, оказываются несовместными. Так, например, нельзя совмес- 
тить режимы стабилизации угла ангажа и высоты в продольном 
канале В боковом канале несовместимы режимы стабилизации уг- 
ла крена и угла рыскания. 

Система уравнений, включающая интегратор и внутренний контур 
<к -0, 93<?5 * О, 7295** -0,11*6 й ; 

'і>, - - 5,27 *о), - 9,7М*<<- в,7*д ь ; 

І7 \ 

эе = і?. 

где ді - переменная интегратора. 

Желаемая динамика определяется корнями : 

* м - -3,0 + у2,0, *= -0,85, «= -0,75 , 

что приводит к строке козіфіциѳнтов обратной связи 

і = [1,035, -0,191, -2,064, — 1 , 4Е ] 
и ограничениям в интеграторе |гё| 6 I, \эе\ 
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3. 


Ресчет контуре стабилизации высоты осуществляется с по- 
мощью системы уравнен. й 

Я - 0,93954,10 * - 0,7295* ы - і7, 1і . <5, ; 

<о г = - V 7 * а), - 0 1 7ЧЧ * < -8, 7*б й ; 

О - ; 

Н ■= 207,0*0 ; 

<к * /У 

Для размещения корней замкнутой матрицы динамики в точках 
5 і<2 = -ДО ±у'-,0 . $ 3 г “0.45, 5 у = -0,35, 5> « -0,3 

необходимы коэффициенты обратной связи: 

і = [1,056, -О, И, -І.'ДЗ, -0,0041, -0,00056 } . 

Ограничения интегратора: |ае|* 300,0 м/с, |эе| « 40000,0 к. 

Для расчета контура стабилизации угла крэка используется 
система уравнений 

- - 0,05 Ч* о)_ -і?0'3* р - /7,0* 6 2 ; 

и) х - Д 366 *сО у - 0, т *р+ 0,82* д } ; 

/3 = 0,0т - * о) х - с^7 * а) у - О, ?і*/і ; 

І =С0 Х - 0,078* о) у ; 

* “ Г- 

Желаемая динамика контура управления определяется корнями: 

5, =-3,5, « ^ 3 --1,?5±^2,25, $ ч -~1,25, 5 У «~/,3. 

Строка коэффициентов обратной связи: 

/ = [-0,166, -0,127, -0,283, -0,639, -0,342]. 

Ограничения интегратора составляют: ( і Т , | ре | 3. 

Для расчета контура стабилизации угла рыскания нолсльвуат- 
ся система уравнений: 

й) х - - <7,05? > ■ 0, 45 *,3 - /7,0*6, ; 

о) у - *(О х - 2,306* оду - 6, ЧЧ5*р+0,В2* д 3 ; 

р - ц) х ■*- <7,^7* - 0,73 */» ; 

Г = “>* “ 407* * ^ ; 
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ѵ^у * 

эе « ф . 

Желаемая динамика определяется корнями: 

5, - - 3,5, з гі - -1,35, 2,25 , ^ = - 1,25, з 5 — /,3, 5 6 --0,15. 
Строка коэффициентов обратной связи: 

/. « [-0,153 , 0,316,-7,71,-0,082 , 8,32, 1,15] . 
Ограничения интегратора составляют: } эё | 6 0,15, )эе| 4 40. 

Заключение 

Рассмотрены примеры синтеза основных контуров управления 
движением самолета методам модального синтеза. 



Рис.З. Совместная отработка входных воздействий в про- 
дольном и боковом каналах управления 


Проверка качества управления проводилась путем численного 
решения системы уравнений (I) на ЭВМ "БЭСМ-6". Примеры процѳс- 
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сов отработки программных значений по углам кренп, тангажа, 
рыскания и высоте представлены на рис. 3. 

Полученный опыт подтверждает возможность использования 
программного обеспечения метода модального синтеза для созда- 
ния САПР систем управления движением самолета. 
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УДК 681.3 Е.Г. Беленький (ЦНИИГАиК) 

МОДЕЛИРОВАНИЕ РАЗРЕШШХ СИСТЕМ НОРМАЛЬНЫХ УРАВНЕНИЙ 
РАЗЛИЧНОЙ ОБУСЛОВЛЕННОСТИ 

Решение систем нормальных уравнений является одним из наи- 
более трудоемких вычислительных процессов в уравнительных вы- 
числениях, особенно при большом количестве определяемых пара- 
метров. При решении больших систем нормальных уравнений матри- 
ца системы, К8К правило, является разреженной. Для эффектив- 
ной организации решения требуются разработка и применение спе- 
циальных методов, алгоритмов и математического обеспечения, 
использующих структуру разреженной матрицы и минимизирующих 
количество вычислений и объем необходимой памяти. Эти вопроси 
являются весьма актуальными сегодня как в практике уравнитель- 
ных вычислений, так и в развитии численных методов линейной 
алгебры. 
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.Для решения систем нормальных уравнений большого поряд- 
'-і , возникающих при уравнивании геодезических, и при 

аналитической блочной фототривнгуляі>ии , разработано ма- 
тематическое обеспечение, которое позволяет решать разрежен- 
ные системы нормальных уравнений с матрице», «и профильной и 
профильно-блочной структуры как в оперативной памяти, так и 
при использовании внешней памяти на ВД. Математическое обеспе- 
пение носит универсальный характер и может быть использовано 
при решении любых задач, приводящих к разреженным системам 
нереальных уравнений денной структуры. 

Профильно-блочная структура матрицы выбрана в связи с 
тем, что зто наиболее общая структура матрицы нормальных урав- 
нений, которая возникает при уравниваніи: блокоанѳлитической 
ф о то три а і і гу ляци и и геодезической сети с нерегулярной структу- 
рой. Отметим, что регулярная структура сети приводит к ленточ- 
но- о'лсчной матрице, что является частным случаем профильной 
структуры. Ь математическом обеспечении предусмотрена возмож- 
ность получения и сокращения профиля матрицы с помощью плио- 
ритма локальной минимизации профиля [і ] . 

В дальнейшем будем предполагать, что в памяти ОШ хранит- 
ся по строкам блок за блоком нижний треугольник матрицы нор- 
мальных уравнений, причем нулевые блоки вне профиля не хра- 
нятся и не участвуют в вычислениях. 

Для решения системы ’ 


АХ = Б 


(I) 


используется известный вариант метода Холецкого, модифицирован- 
ный для обработки блоков, в котором выполняется треугольное 
разложение матрицы 



( 2 ) 


где X — нижняя треугольная матрица с единицами на главной 
диагс элн; I) - днаюнальнвя матрица. Вектор решения опреде- 
ляется в результате последовательного решения систем 


і 2-5 , ВУ~2, Х Г Х= V, 


(3) 


где решение первых двух систем реализует прямой ход, а после- 
дней - обратный ход. Для этого метода характерны высокая точ- 
ность и устойчивость численного процесса. Алгоритм решения 
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системы при использовании внешней памяти и основные формулы 
приведены в [2] . 

При разработке математического обеспечения для решения 
разреженных систем норме -ьных уравнений для тестирования, 
оценки быстродействия алгоритма вычислений, оценки устойчи- 
вости метода решения необходимо уметь моделировать тестовые 
системы нормальных уравнений различного порядка, разреженнос- 
ти н обусловленности. Генератор тестовых задач является очень 
полезный инструментом при таких разработках [3] . 

Здесь будут изложены алгоритм и возможности генератора 
тестовых систем нормальных уравнений, даны некоторые оценки 
числа обусловленности матрацы, а также приведены результаты 
численного моделирования на ЗИЛ. 

Под тестовой системой нормальных уравнений будем понимать 
систему нормальных уравнений с известным решением. Обычно тес- 
товые системы получают, исходя из решения модельных задач, 
приближенных к реальным. Однако для всестороннего тестирования 
математического обеспечения и анализа его эффективности , осо- 
бенно, если математическое обеспечение имеет широкое назначе- 
ние, наиболее целесообразным является непосредственное модели- 
рование систем нормальных уравнений произвольного порядка, 
структуры и обусловленности. 


л і о С • 



Рио.І. Треугольное разложѳі"" 3 профильной матрицы 

Такой подход используется в предлагаемом генераторе тесто- 
вых систем. Он позволяет получать в оперативной памяти ЭШ 
| системы нормальных уравнений с профильной и профильно-блочной, 
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8 также с ленточной и ленточно-блочной структурами матриц 
различного порядка и обусловленности. 

При моделировании тестовой системы (I) необходимо вычис- 
лить положительно определенную матрицу А заданной или моде- 
лируемой структуры и вектор правой части системы В. Матрицу 
А о заданным профилем можно получить из формулы (2), исходя 
из того, что структура профиля и ширина ленты у матрица А 
и А одинаковы (рис.І). Задавая профиль (ширину ленты) мат- 
рицы А , а также элементы матриц А и І> , можно вычис- 
лить матрицу А заданной структуры. Таким образом, алгоритм 
моделирования разреженных систем нормальных уравнений заклю- 
чается в следующем. 

I. Моделирование профиля матрицы 

Если профиль или ширина ленты задаются, то этот пункт сле- 
дует пропустить. При моделировании длина і - й строга про- 
филя определяется по формуле 

е ; - тіп{« а (і-1)Е 1 + /,<], (4) 

где к 0 - коэффициент разреженности матрицы; &■ - псевдо- 

случайная величина, равномерно распределенная на интервала 
(0,1) ( Е- <г Я (0,1)). Уменьшение к 0 < I увеличивает разрежен- 
ность матрицы. При к ѵ >> I моделируемая матрица нормальных 
уравнений будет заполненной. 

С помощью формулы (4) определяется структурный массив про- 
филя Л'РЛ : 

НРЛ(і)-1, /№>(()• ЮФ(1- Г) + Ъ, (5) 

где величина Л'РЛ(і) указывает номер диагонального элемен- 
та С -й строки в одномерном массиве , в котором хранится 
профиль. 

2. Моделирование треугольной матрицы А 

Элементы 6ц (і »/) также моделируются с помощью датчи- 
ке псевдослучайных чисел с равномерным распределением, причем 
опродѳліяются толіько элементы, лежащіе в пределах профиля по 

формуле 

6ц - (в~а)е ц , 

- І п 4 - 


( 6 ) 


где 8 и а - заданные величины {8>а ); е г е# (0,1). 

Из (6) следует, что 8 ■■ е. 9 (а, 8). 

3. Моделирование диагональной матрицы 2) 

Элементы аі; вычисляются но формуле 

Ч , (7) 

где сі п сі п - наименьший и наибольший элементы диагональной 
матрицы (сі п і й, *0), е ( еВ(0,І), (/,• в к(аІ п с/^). 

4. Вычисление матрицы А по формуле (2) 

5. Определение вектора правой части В 

Вектор правой части вычислялся по формуле В=АХ Т , где Х т - 
точный вектор решения, все элементы которого принимались ров- 
ню® I. По изложенному алгоритму разработан комплекс программ 
на языке ФОРТРАН -ІУ, который реализует генератор тестовых 
систем, позволяющий моделировать профильные (ленточные), в 
также профильно (лэнточно) - блочные системы нормальных урав- 
нений с одинарной и двойной точностью. 

Профиль матрицы может моделироваться по формулам (4) и 
(5) или задаваться непосредственно. При моделировании о оди- 
нарной точностью все вычисления выполняются с двойной точностью. 
В качестве датчика псевдослучайных чисел с равномерным распре- 
делением использовался датчик, рекомендуемый в [4] . 

Для исследования обусловленности получаемых систем нормаль- 
ных уравнений и для определения интервалов моделирования мат- 
риц і и 2) для вычисления матрицы А заданной обусловлен- 
ности был выполнен ряд численных экспериментов. Прежде чем 
перейти к изложению результатов, рассмотрим некоторые оценки 
обусловленности системы нормальных уравнений. Количественной 
мерой обусловленности систем нормальных уравнений является 
число обусловленности матрицы А 

НА)-ІА{ г ІА-%-А,/Я„ , ( 8 ) 

где А 1 и А„ - соответственно наибольшее и наименьшее соб- 
ственные значения матрицы А. Однако вычисление к(А) по фор- 
муле (В) требует значительных вычислительных затрат, поэтому 
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для определения величины к (А) использовались более "дешевые" 
оценки. Хорошей оценкой к ( А) является следующая величина: 

к, (А) « а тах а пах 6 к (А), (9) 

где а тх 0 ^/пах максимальные диагональные элементы матриц 
А и А' 1 соответственно. 

Для другой оценки к (А) введем следующие определения. От- 
носительной точностью решения системы будем называть величину 

9-ІЛХ-Ві/ійІ, • (Ю) 

Л 

где X - вектор неизвестных, полученный в результате реше- 
ния системы. Относительной ошибкой решения будем называть ве- 
личину 

?*3х-х г і/іх г /. (и) 

Известно [5] , что величины ^ и ^ связаны между собо* 
следующим соотношением: 

ц* к (А)у. (12) 

Таким образом, 

/<(А)>ц/у. (13) 

Величины ^ ир можно определить в результате решения 
системы, а величину 

к г (а) «■ р/р . (іі) 

использовать в качестве нижней оценки числа обусловленности. 
При решении тестовой системы оценку (14) нетрудно вычислить, 
так как точное решение системы известно. 

С помощью генератора тестовых систем моделировались сис- 
темы нормальных уравнений о ленточными матрицами порядка N 
и шириной ленты уЗ ( а ^ = 0, если |і-/| * /і ), которые 
затем решались по формулам (2) и (3). Моделирование и решение 
систем выполнялись на ЭВМ ВС-1045. 

В табл. I приведены результаты моделирования и решения сис- 
тем порядка N = 250 и шириной ленты уЗ = 50 в зависимости 
от длины интервала моделирования элементов матрицы Ь л ( =в-а 
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Результаты моделирования 

Таблица I 


N 

1 

Ч 

^тіп 

^ тах 


1 

к г (А) 

I 

0 

0,5 

1,0 

1,0 

4,0 

І,0ЛСГ іа 

24,9 

2 

0 

0,75 

1,0 

І.с 

2,8 

6,4-ПГ 11 

2, ЗЛО 4 

3 

0 

1,0 

1.0 

1.0 

2,3 

9,0-ІСГ 6 

3.9Л0 7 

4 

0 

1,25 

1.0 

1,0 

3,5 

1,1 • КГ 2 

3,1 ЛО 12 

5 

0,5 

0.5 

1.0 

1,0 

5,1 

3,6ЛСГ П 

7,0 ЛО 3 

6 

1,0 

0,5 

1,0 

1,0 

3,1 

3,5 -І0“ 9 

І.І-ІО 6 

7 

1,5 

0,5 

1,0 

1,0 

5,4 

7.0-І0* 4 

І.ЗЛО 11 

8 

0 

0,5 

1,0 

10000,0 

3,5 

7,9 -ІО -13 

229,3 

9 

0 

1.0 

1,0 

10000,0 

3,7 

1,4- ІО -6 

3,7- ІО 8 


Таблица 2 


N 

N 


ё 




к г ( А') 

I 

500 

50 

0,0 

0,5 

4.0 

3,6 л о~ 13 

09,8 

2 

1000 

50 

0,0 

0,5 

4.0 

2.0Л0" 12 

507,3 

3 

1000 ' 

20 

0,0 

0,5 

1,4 

І.6Л0 -15 

1,3 

4 

500 

50 

0,0 

0,8 

2,6 

7,3 ЛСГ® 

2.8Л0 9 

5 

250 

25 

0,0 

0,5 

2,5 

4,3 Л О" 15 

1,7 

6 

250 

100 

0,0 

0,5 

5.9 

2, 9 Л О -13 

49,3 


и их среднего значения і *(а+8)/2 , а также от диапазона 

изменения элементов матрицы Ь . Как видно из таблицы, гене- 
ратор тестовых систем позволяет моделировать матрицы практи- 
чески любой обусловленности от очень хоросей (примеры I, 8) до 
очень плохой (примеры 4, 7). Из табл Л следует также, что уве- 
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Таблица 3 



лнчение интервала моделирования Л ( или смешение значения I 
относительно нуля приводит к ухудшение обусловленности мат- 
рицы. В то ке время увеличение диапазона изменения элементов 
матрицы 2) влияет на изменение обусловленности незначительно 
(для сравнения моино взять прим. 8,9 и 3). 

Выбор параметров моделирования для получения матриц заданной 
обусловленности зависит от порядка матрицы и размера профиля 
или ширины ленты. Для анализа этой зависимости было проведено 
моделирование оистем различного порядка и ширины ленты при 
фиксированных параметрах моделирования. Результаты представ- 
лены в табл. 2. На обусловленность матриц с узкой лентой(см. 
прим. 3,3 из табл. 2) сильнее влияет ширина ленты, чем порядок 
матрицы. Для матриц о сравнительно широкой лентой (см. при». I 
из табл.1 и 6 иэтабл. 2) эта зависимость не так ярко выражена. 
Следует отметать, что скорость изменения обусловленности сис- 
темы от изменения порядка матрицы (ширины ленты или размеров 
профиля) зависит от выбора параметров моделирования'. Так, напри- 
ример, использование параметров моделирования из прим.З табл.І 
для матрицы порядка N =500 приводит к настолько плохо обуслов- 
лов ленной системе, что получить ее решение оказалось невозмож- 
ным. Полученные результаты подтвердили высокую устойчивость ме- 
тода решения по формулам (2)и(3). Как видно из табл. 1,2 величи- 
на у , характеризующая количество верных знаков в векторе пра- 
вой части, получающегося при подстановке вычисленного решения 
в систему, имеет порядок 10 и лежит на границе разрядной сетки 
при вычислениях о двойной точностью. 
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этом величина р практически не зависит ни от порядка и 
ширины ленты матрицы, ни от обусловленности системы. 



Рио. 2. Графики накопления ошибок 

Для сравнительного анализа устойчивости было выполнено 
моделирование и решение систем Нѵ_ лальных уравнений различной 
обусловленности с помощью формул (2), (3) и при использовании 
обращения матрицы. Результаты приводятся в табл.З. Для всех 
примеров было взято N * 100, = 20 . Из табл.З видно, что 
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решение с помощью обращения матрицы становится неустойчивым 
при ухудшении обусловленности системы, хотя ошибки решения име- 
ют сопоставимые величины. Из табл.3 следует также, что вели- 
чины к, (А) и к : (А) хорошо согласуются между собой. 

Анализ величин у и оценок числа обусловленности в 
табл.І - 3 подтверждает эмпирический закон, согласно которому 
из-за ошибок округления в процессе решения теряется Рдк(А) 
десятичных знаков [6] , а именно: 



(15) 


Где вместо к (А) нужно подставить полученные оценки. 

По результатам моделирования был выполнен анализ накопле- 
ния ошибок в процессе решения в зависимости от обусловленнос- 
ти системы. На рис. 2 представлены зависимости ошибок решения 
от номера элемента в векторе решения для примеров 1-4 из 
табл.І. Значения модулей ошибок приведены в логарифмическом 
масштабе. Данные зависимости строились по точкам через каждые 
12-13 номеров, поэтому полученные графики являются сглаженны- 
ми по сравнению с истинными колебаниями модуля ошибки, однако 
характер накопления не зависит от этих колебаний. Из рис. 2 
видно, что чем хуже обусловлена матрица, тем сильнее влияние 
накопления ошибок, причем имеет место тенденция возрастания 
ошибки от последнего элемента вектора решения к первому. Харак- 
тер накопления ошибок связан с выполнением обратного хода. 

Тагам образом, реализованный алгоритм моделирования раз- 
реженных матриц позволяет получать положительно определенные 
матрицы практически любой обусловленности. Результаты модели- 
рования показали, что увеличение длины интервала моделирова- 
ния и смещение среднего значения элементов матрицы относитель- 
но нуля ухудшают обусловленность моделируемой системы и нао- 
борот. 

С- помощью данного генератора тестовых систем был выполнен 
сравнительный анализ устойчивости при решении системы с помощью 
треугольного разложения а обращения матрицы и показана числен- 
ная неустойчивость последнего способа, а также выполнен анализ 
накопления ошибок в векторе решения при вычислениях по форму- 
лам (2) и (3). 
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Данный генератор можно попользовать также в задачах ста- 
тистического оценивания для моделированія априорных ковариа- 
ционных матриц параметров, реальные ковариационные матрицы 
которых неизвестны. 
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УДК 62-50 К. А. Селиванова (Уральский по- 

литехнический институт) 

АНАЛИЗ ДВИЖЕНИЙ НЕЛИНЕЙНОЙ ЦИФРОВОЙ СИСТЕМЪ: МЕТОДОМ 
СЕКУЩЕЙ ОБЛАСТИ 

Рассматривается система стабилизации, объектом которой яв- 
ляется интегрирующее звено второго порядка. Управление осущест- 
вляется цифровым регулятором. Анализ динамики системы прово- 
дится методом секущей области и базируется на работах [і, 2] , 
где алгоритм стабилизации основывается на измерении выходной 
координаты объекта. В данной работе для обеспечения стабили- 
зации САР предусмотрено измерение первой производной коорди- 
наты аналоговым датчиком. 
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I. Структура системы стабилизации и ооновныѳ расчетные 
соотношения 


Структурная схема системы стабилизации представлена на 
рисунке. Здесь С - коэффициент, характеризующий эффектив- 
ность исполнительного органа; Т 0 - период квантования по 
времени; о (і) и &*(і) - соответственно функция управления 
и приведенная функция управления; й[к] и и [и] - управ- 

ляющее воздействие и приведенное управляющее воздействие; а д 
и а, - параметры управления; с(€) - выходная координата. 


Переключение реле происходит при выполнении следующих 
условий: 


!«• (0 - > С<. 

'0 

(I) 

Уравнение движения объекта 

кмѳѳт вид: 


т- 

\СТ„ , 


0 , 

(2) 


-ст . . 


Функция управления: 



в(і)-а в г(і)+а,г(і). 

(3) 


Как видно из риоуяка, и [к] принимает дискретные значения 
-І.О.+І. 

Для анализа динамики системы, следуя (I), запишем конеч- 
но-разностные уравнения системы: 


х,[н+1]-х,[к]+и[х] ; (4) 

т г [«*>]-т<[*] + х г (к]. (5) 

Это нетрудно сделать, записав выражение для импульсной пере- 
даточной функции приведенной непрерывной части системы и вво- 
дя обозначения: 




и(з ) 
г-і 




Переменные Х 1 [к] , І г [к] 


и (2) _ х,(7) 

(7-1)* ~ 2 - ' 
связаны с координатами 


ѵ(і) « г (О соотношениями : 


ѵ[к]-Т 0 Сх,[к] ; 

(*Ъ [*]+*, Р*і]- 
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3(0 И 

(6) 

(V) 


м/7 1 м/7 7ГТ г ' м.?- 

* і 



Введение координат 
х, и т г позволя- 
ет упростить про- 
цесс расчета движе- 
ний в системе за 
счет того, что при 
изменении дискрет- 
ного времени к на 
единицу, х, изме- 
няется нс один и 
тот же постоянный 
по величине квант, 
равный +1. 

Легко показать, 
что 


+ 2а 0 (8) 

Вводя обозначения 


*а 


А .ІЙі 

А г ш 2а 0 ; 

2ы. 

Тт7 


« і 


€'ш 


( 10 ) 

(И) 

( 12 ’ 


представим функдаю 
управления в виде: 

9[*) -*,*,[*}* 

+ Л,*,[м], (із) 


- из 


а условия переключения 


Ів[и]І > е. 

(14) 

Из соотношений (І3)-(І4) определим выражения для 
ключѳния (ЛІ): 

линий перѳ- 

ЛИ х г .-±-Ліх,і 

. л* 

Ца) 

лі2 <?/*; - е , х г - - А х, . 

( 16 ) 

Согласно [I] уравнения движения объекта при постоянном 


значении приведенного управляющего воздействия ы[к] опреде- 
ляются формулами: 
при и Ф О 

х, [м] • х, [О] + к идя и[*-(] ; (17) 

х г [к] •т і [о]*к$ідли[к-і] ^-+нх,[о] ; (іе) 

при и - О 

х,[к] - х,[0 ] ; (19) 

х г [и] - их, [О] *х г [0]. (20) 

Уравнения фазовых кривых системы соответственно при и+0 
(21) и и- О (22): 

х г ~ С 0 ]' } *,[<>](*, [0]*У»и - ^ Х?зідпи-^Х, ; (21) 

Х,п х,[0]. (22) 

Из (21) следует, что точки (х,[и ] , х г [и]) 'располагаются 


на параболах, симметричных относительно прямой х, - ~ $сдп и ■ 

2. Функции последования 

Для исследования движений нелинейной цифровой системы ста- 
билизации попользуется метод секущей области, описанный в [і] . 

Для установления характерных типов движений в системе до- 
статочно провести исследование траекторій, начало которых 
находится в полосе 5 : 

{*,*(- 0,5, 0, 5), т г 

Координаты произвольной точки в области 5 обозначим: 6= х, , 


х г • Область 5 выбрана таким образом, что все тра»»., 
тори движения системы, начинающиеся в произвольно заданной 
точке (х п х 2 ) плоскости (х и х 2 ) , попадают в эту область и 
покидают оѳ за один шаг. Такой способ выбора области $ при- 
водит к тому, что вместо исследования поведения системы при 
произвольно заданных начальных условиях достаточно исследо- 
вать поведение системы при начальных условиях X, [О] - б , 

Х 2 [°] “<Г- 

Назначая в области 5 произвольную точку М( у, б~) , 

можно определить о гображенпѳ М(р,д) этой точки в область 
5 и, следовательно, получить фу нкцию последования М - { (м), 
М, М е 5 ; 

л- Л Гг. Л 

Решения системы (4), (5), (13), (14) таковы, что если на- 
чальная точка ( т і [0], х г [О] ) принадлежит линиям 
іб = [і 1 [0] = б 1 бе (-0,5, + 0,5)} , то любая траекторія с 

ростом времени бесконечное число раз возвращается на ось . 
Это дает основание свести анализ динамики системы к изучению 
функций последования вида р ? б-сопзі-. 

Укажем характерные виды траекторий движения точка 
-Т*> (т,[н] ,х г [к]'). Начало траекторіи находится в точке т,[о]- 

= б, х г [ б ] - г- 

1. Точка Т в течение т тактов дискретного времени пе- 
ремещается под действием управления и = -I, проскакивает за 
один такт зону нечувствительности и далее за т тактов о 

и - л-1 достигает линии / у . 

2. Точка Т движется по траектории т тактов с и = -I, 

Ч тактов с и = 0 (в зоне нечувствительности) и т тактов о 
и = +1, возвращаясь на линию . 

3. Точка Т достигает линии в режиме скольжения от- 

носительно линии переключения. 

Можно показать Гі ] , что функция пооладсвания для второго 
из указанных случаев (у 1 0) определяется выражением 

р ~^*т(2б-т) + у(б-т). (23) 

Если в зоне нечувствительности точка Т не задерживается 
(? = 0), то 
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X = Х* т (2& ~ т )- 

функция последования является нечетной функцией обоих аргумен- 
тов, т.ѳ. 

/ (Г, & -<*)■ ( 24 ) 

Это обстоятельство позволяет проводить исследование движений, 
зарождающихся на линии при б е (О, 0,5] . 

Как показано в [2] , функция последования не обладает 
свойством непрерывности на всем учаоткѳ изменения координаты 
. Существуют точки разрыва функции последования, т. е. 
точки, которым будут соответствовать изменения числа шагов до 
моментов перемены знака и величины управляющего воздействия. 

Точки разрыва функций последования можно разделить на 2 
группы. К первой группе относятся то^ки а ^ • 

В указанной записи знак - (или +) указывает начальное значение 
и[0] • Точки обладают следующим свойством: при 

Г = І~/т) ~ О точка Т движется о и[к]--1 до изме- 
нения величины управляющего воздействия в течение т тактов; 
при Х я Тѵ*> * 0 изменение управляющего воздействия 
происходит через т + I такт. Ко второй группе относятся точ- 
ки ' определяющие изменение поведения точки 

Т в зоне нечувствительности после т тактов движения при 
и[к] ~ - 1(и[к] ~І) и обладающие свойство.'.!: при Г” 0 

точка Т будет двигаться т тактов с и [к] - - / , попадает 

е зону нечувствительности, а затем д тактов будет двигать- 
ся по жнии х, • б -т ; при у» у™ + 0 в зоне не- 

чувствительности точка Т будет двигаться (д + <) такт. 

На основании соотношений (І5)-(22) легко получить выраже- 
ния для точек разрыва функции последования: 




‘ (м) 

< * 

И + > ~ (д*пі)-тб ж 

1 <*» А г л ' 


- 

'"М> Л г А г 1 " 2 


- - 4 1 (б-т)-тб+ 


т(т~ 1 ) . 
2 * 

т(т-і) 


- д(б-т); 




(25) 

(26) 

(27) 

(28) 
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Таким образом, для построения функции последования, соот- 
ветствующей фиксированным значениям параметров а, , а , , <л , 
необходимо: 

1. Вычислить точки разрыва и по формуле 

(25) начиная с т = 0. 

2. Последовательно вычислять точки 0,1,2,..., 

отбрасывая все точки )Г^'. ^ . которые не' принадлежат отрез- 
ку ГГ ( Ц > ііт-к) ] • 3 Результате будет получен упоря- 
доченный набор точек Г = /V - . У**' . „ . .. 

^ 1 ѵ т,т*і II (Л > I іт+іур > 1 » 

)(/»»/) '• 

3. Для всех точек ре да,/ по формуле (23) определя- 
ются значения р справа и слева от р . Значения р сле- 
ва от точки Р г, 

р - С + (/и+ 1\2д -т- /) + у(д-т-І), 


оправа: 

р = ^ + (т+1Х2д-т-1)+((і+1)(д-т-і). 

Значения р слева от точки у - у^ при т + О : 

У= Г+ т(2д-т)*(ц^ + 1)(д-ю), 

где ^ - максимальный индекс ^ из упорядоченного набора 

- Значения р справа от р. Г( ^ : 

р • у*(т*1)(2д-т-і)*()^ (д-т-і). 

Здесь ц р - минимальный индекс из упорядоченного набора 
Г 

При т = 0 значение р слева от точки равно: 

4. Па плоскости (р, р) в точках разрыва отображения 
точки р соединяются вертикальной прямой, а между точка- 
ми разрыва - наклонными линиями. 

5. Для значений ^ е значения р вычисляются на 

основании свойства (24) функций последования. 


3. Анализ процессов в САР 

Прежде всего целесообразно установить предельную ограни- 
ченность процессов в САР. 
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Условие продельной ограниченности процессов в САР может 
бить записано в следующем виде [ I ] : 

я *2 (29) 

где - координаты точек разрыва функции последования; 

р + - значение функции последования справа от точек разры- 

В8 ! • 

г-г(г<% ♦«)• 

Условие предельной ограниченности движений в системе с 
учетом (10)— (12) , (23), (25) может быть преобразовано к виду: 



Выполнение условия (30) влечет за собой [I] существование для 
данной системы ограниченного асимптотически устойчивого мно- 
жества значений ( х , , х г ) , к которому стремятся в течение 
времени , либо входят в него и больше его не покидают все тра- 
ектории системы (4)-(5). 

Выделение области значений координаты ре ід , в кото- 
рую сходятся движения системы (4)-(5), может быть осуществле- 
но на основании анализа функций последования. Располагая соот- 
ветствующими программами для проведения расчетов на ЦВМ, мож- 
но получить семейство функций последования при различных пара- 
метрах управления о 0 , а , . Анализ полученных зависимостей 

позволяет сделать выводы о характере возможных движений в 
системе. 

При выполнении условия (29) для любого к величина 
ІГМІ . Усилим условие (29), выбирая некоторое конеч- 

ное значение т , что приведет к уменьшению размеров облас- 
ти 8 х , ограничивающей при к -» все движения систе- 

мы. . 

Б качестве примера рассмотрим расчет области параметров 
управления а 0 , а 1 для которых условие 

ГѴгад >-і- 

выполняется для всех , При этом все движения системы 

будут сходиться в область , в которой максимальное ко- 

личество шагов с ненулевым управлением при перемещении точки 
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Т с линии ід на линию равно 4: 2 шага о 

и = , д шагов с и- 0, 2 шага с и- + 1(и~-{). 

Расчет области параметров управления проведем по следую- 
щей методике: 

I. Рассмотрим упорядоченный набор Г- г т при т = 3: 


Г 4 г *' 0 гѴ у'->] 

т-1,т I і (т-і) 1 о(/п,у) ’ *(т.<і+і) * « пл) 1 ' 


»?> 

2. Последовательно назначим д = 0,1,2 полагая д - 

минимальным индексом точки разрыва отображения в 

наборе Г" , т . В зависимости от значения д ’ определим 
область параметров А ^ 

3. В соответствии с (і] область параметров определяется 
выполнением следующих соотношений: 


Гш.ѵ > Г, "-л і (30 

Г, ‘ Г,'1„ ■ (32) 

г*-гад..,-«і (зз) 

Г-)Чг<іѴ«- (3« 

4. С учетом найденных соотношений для точек разрыва функ- 
ций последования, по формулам (31)-(34) определим параметры 
О д , а , , для которых эти соотношения являются верными. 

Рассмотрим, например, определение границ области А 0 
В соответствии с (31), (33), (34) должны выполняться условия: 


Га») > ГЙ і (35) 

Г и) > ~ Г /о ; (36) 

У,й > - Г&) • (37) 

Учитывая выражѳни (23) для функции последования и выраже- 
ния (25), (27) для точек разрыва отображения, получим соответ- 
ственно следующие выражения: 

Е < -зг + а ѵ (2 1 5-б) • 

'о 

е > (6-2)+ а „(5- 6)-, 

'о 
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С < — уг- (д-3) - іа с • 

'О 

При фиксированном значении е - / (ы) эти условия определя- 
ют на плоскооти параметров (а 0 , а ,) область А „ 

Аналогичным образом определяются области А ? при д + 0 
(при фиксированном /п ). Объединение областей А ^ на плос- 
кости (а в1 а,) дает область А , представляющую собой 
множество параметров а 0 , а , , которые гарантируют сходи- 
мость всех движений в область , в которой максималь- 

ное количество шагсэ с ненулевым управлением равно 4. 

Границы каждой области А ^ зависят от значения д 
Область, являющаяся пересечением областей Ац при 6 *= С 
и 6 «■ 0,5, является областью параметров а „ , а , , удовлет- 
воряющих неравенствам (31)-(34) при всех де [0,0.5 ] . 

4. Заключение 

В статье рассмотрены вопросы построения функций последо- 
вательности для релейно-импульсной системы стабилизации, 
объект которой описывается передаточной функцией Ѵѵ7/>) - ~ г . 
Управление объектом осуществляется на основе измерения выход- 
ной координаты объекта и ее производной. Описано построение 
областей А * на плоскости параметров (а 0 , а,) . Выбирая 

параметры управления из этих областей, можно обеспечить опре- 
деленный тип установившегося движения в нлосі ,сти ( т 1 , . 

1. Косякин А. А., Шамриков Б.М. Колебания в цифровых авто- 
матических системах. М. , 1963. 336 с. 

2. Шамриков Б.М. Качественное исследование нелинейных им- 
пульсных систем методом точечных отображений // Изв. АН СССР. 
Техническая кибернетика, й 3. 1972. С.І5-20. 
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РАЗРАБОТКА ИЖТРУМЕНТАЛЫЮГО СРЕДСТВА ПОСТРОЕНИЯ 
ДИАЛОГОВОЙ СИСТЕМЫ, ОРИЕНТИРОВАННОЙ НА ПОЛЬЗОВАТЕЛЯ- 
НЕПРОГРАММИСТА 

В настоящее время не вызывает сомнений необходимость соз- 
дания и широкого внедрения человеко-машинных диалоговых систем 
(ДС) , благодаря которым эффективно используются средства вычис- 
лительной техники для решения многих практических задач. Осо- 
бый интерес [I] представляет разработка диалоговых систем, ко- 
торые требуют от пользователя в первую очередь специальной 
(профессиональной) подготовки и во вторую - подготовки в об- 
ласти программирования. Как средство диалога для неподготов- 
ленного пользователя широкое распространение получили меню- 
системы. В меню-системах у пользователя создается иллюзия 
взаимодействия с ЭВМ на своем профессиональном языке. Считает- 
ся целесообразным создавать инструментальные средства построе- 
ния меню-систем, позволяющие настраиваться на любую предмет- 
ную область. Примером инструментального средства может послу- 
жить меню-система подсистемы настройки на заданную объектно- 
технологическую среду робота, входящая в состав "АКИ проекти- 
ровщика программ управления робототехничзским комплексам 
ИЛ-01 \ 

Цель разработки подсистемы настройки заключается в постро- 
ении диалога между пользователем и ЭВМ, ориентированного на 
пользователя-нѳпрограммиста и обеспечивающего настройку АРМа 
на заданную объектно-технологическую с;реду, а также подстрой- 
ку структуры диалога в процессе работы при переходе на новую 
объектно-технологическую среду. Подсистема обеспечивает ввод 
информации, описывающей расположение технологических зон, гео- 
метрическую форму стационарных и перемещаемых объектов в рабо- 
чем пространстве робота, в том числе и геометрию самого робо- 
та, что составляет предметную область системы. Требование мо- 
дификации архитектуры ДС, начиная от изменения сценария диа- 
лога и реплик диалога и кончая подключением новых функций , оп- 
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рѳделяло характер и организацию разработки программного обес- 
печения ДС подсистемы настройки. 

Основу построения ДС, удовлетворяющей перечисленным тре- 
бованиям, составляют выделение и объединение функций управле- 
ния диалоговым процессом в отдельный блок - программу управ- 
ления диалогом (ПУД). Жесткий формат кадра диалога и построе- 
ние сценария диалога на основе использования структуры фрей- 
мов обеспечивают независимость ПУД от сценария диалога и па- 
кета прикладных программ пользователя. Введем определения 
понятий, необходимых для формализованного описания диалогово- 
го взаимодействия человека с ЭВМ и использующихся в этой ра- 
боте. Под диалогом следует понимать процесс структурированно- 
го обмена сообщениями между человеком и ЭВМ, направленный на 
обеспечение решения конкретной прик адной задачи. Шаг диало- 
га- это часть диалога, включающая следующие действия: вывод 
сообщения из ЭВМ, анализ сообщения пользователем, ввод чело- 
веком сообщения в ЭВМ, обработка введенного сообщения маши- 
ной. Диалоговое сообщение - это алфавитно-цифровая или графи- 
ческая информация, которой партеры обмениваются на каждом 
шаге диалога. Кадр - это формат сообщения при его выдаче на 
экран дисплея. Информация, содержащаяся в кадре, может быть 
дополнена правилами реакции ЭВМ на те или иные ответы пользо- 
вателя на данном шаге диалоги. Фрейм - это кадр, нополнѳнньй 
праг іламв поведения ЭМ на данном шаге диалога. Сценарий диа- 
лога - описание на специализированном языке полного набора 
фреймов и структуры диалога [2] . 

Для описания сценария диалога подсистемы настройки исполь- 
зуется понятие кадра. Кадр имеет жесткую структуру. Стандар- 
тизация расположения различных видов информации в кадре умень- 
шает затраты времени на ее анализ человеком. В кадре выделе- 
но три поля (рис I): поле вопроса или предложения выполнить 
действия по формированию ответа; поле меню, которое содержит 
правило формирования ответа или список допустимых ответов; 
поле ответа . 

Примеры кадров сценария диалога выглядят следующим обра- 
зом: 

I Введите имя технологического поля 

Имя - последовательность букв и цифр, начинающаяся с 
буквы, длиной 6 символов. 
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Введите имя: 

2. Введите тип технологической операции 

I. Взять деталь 

2 Отдать деталь 

3 Ввдержка 
Введите номер ответа: 

Описание 
сценария диало- 
га подсистемы 
настройки осу- 
ществляется с 
помощью кад- 
ров. Каждому 
кадру в полном 
наборе присва- 
ивается поряд- 
ковый номер. 

Структура диа- 
лога представ- Рис.І. Формат кадра диалога 

ляется в виде 

ориентированного графа, вершинами которого являются номера 
кадров диалога. Дуги графа, выходящие аз одной вершины, отра- 
жают возможные пути продолжения диалога в зависимости от от- 
вета пользователя Набор кадров и граф диалога используется 
при заполнении набора данных фреймов диалога. 

Вершине на графе вместе со всеми выходящими из нее дугами 
соответствует фрейм. Структура фрейма следующая: 


I 

2 

3 

4 

5 

- ^ ■ — ■ ... ^ 


Информационная часть Командная часть 


I - кадр, выдаваемый на екран дисплея; 2 - поле допустимых 
ответов; 3 - поле ответа пользователя; 4 - поле указателей на 
выбор следующего фрейма; 5 - имя подпрограммы обработки ответ! 

Заполнение фреймов осуществляется специально разработанной 
программой в диалоге о пользователем. Укрупненная блок-схема 
программного обеспечения ДС подсистемы настройки приведена на 
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рис 2. Формирование набора данных фреймов выполняет блок ре- 
дактора фреймов. Для этого предусмотрен специальный администра- 
тивный режим работы подсистемы настройки. В административном 



Рио. 2. Блок-охѳма программного обеспечения диалоговой 


системы 
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режиме осуществляется и корректировка сценария диалога. Вве- 
дение в командную часть фрейма имени подпрограммы обработки 
вводимой информации обеспечивает независимость диалога от 
пакета прикладных программ пользователя. 

В состав блок-схемы программного обеспечения ДС входят 
пять блоков. ПУД являотоя организующим блоком и выполняет 
следующие функции: считывает фрейм из набора данных фреймов 
диалога, высвечивает кадр на экране дисплея, считывает и 
анализирует ответы пользователя. В первую очередь выясняется, 
не является ли ответ высокоприоритетной директивой, например 
"Просмотр протокола диалога" или "Конец диалога". В этом слу- 
чае директива интерпретируется и управление передается соот- 
ветствующим блокам - блоку протоколирования диалога или бло- 
ку обработки результатов диалога. В других случаях управление 
передается интерпретатору командных строк. Интерпретатор 
осуществляет сопоставление ответа пользователя с возможными 
для данного состояния ответами. Если при анализе не установ- 
лена допустимость ответа, система классифицирует ответ как 
неправильный, формирует диагностическое сообщение и передает 
управление ПУД. ПУД выдает диагностическое сообщение на эк- 
ран и предлагает повторить действия. В случае правильного от- 
вета интерпретатор выполняет дешифровку командной части фрей- 
ма и возвращает управление организующему блоку. Разделение 
управляющей и исполнительной функций .диалогового процесса по- 
зволяет достичь мобильности программного обеспечения, возмож- 
ности гибко управлять диалоговым процессом, а также эффектив- 
ности и удобства разработки программного обеспечения. 

Предложенный вариант диалоговой системы является экспери- 
ментальным образцом При проектировании использовалась кон- 
цепция "быстрого прототипа" [3} . Дальнейшее развитие ДС мс- 

І жет в первую очередь коснуться интерпретатора командных строк 
Новые функции интерпретатора появятся в том случае, если поя- 
вятся новые типы связей между элементами предметной области. 

В представленном варианте ДС присутствуют два типа связей: по- 
исковая (задает статистические сьізи); действие без выдачи 
какого-либо ответа (ввод информации). Блок редактора фреймов 
и блок обработки результатов предлагается отроить с использо- 
ванием СУЩ. 
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УДК 621.865.8:681.3.06 Э. Л. Пономарев (Уральский поли- 

технический институт) 

ОБЗОР АЛГОРИТМОВ ПЛАНИРОВАНИЯ ТРАЕКТОРИЙ 
ДВИЖЕНИЯ РОБОТА 

Одним из важных направлений исследований в области робото- 
техники является создание средств программированіи промышлен- 
ных роботов на уровне заданий, а не команд При этом одно из 
центральных места этих исследованиях занимает проблема авто- 
матического планирования траекторий движения робота в рабочей 
среде с препятствиями. Рассмотрению некоторых алгоритмов 
планирования траекторий движения робота и будет посвящен дан- 
ный обзор. 

Часто трудность учета всех деталей рабочей среды, возник- 
новение непредусмотренных препятствий вызывают необходимость 
использования методов планирования движения робота в незнако- 
мой среде только на основе информации, получаемой роботом от 
своих сенсорных устройств. Если раачѳры окрестности, в кото- 
рой робот может получать информацию, соизмеримы с размерами 
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зрепятствий, то прокладка маршрутов может быть осуществлена 
о помощью алгоритмов, являющихся дальнейшим развитием мето- 
дов, разработанных для полностью известной рабочей среды. 
Однако в том случае, когда размеры информативной окрестности 
малы, необходим принципиально иной подход. Основная особен- 
ность заключается в том, что после обнаружения препятствия 
робот, оставаясь на прежнем месте, не может оценить возмож- 
ные направления дальнейшего движения и выбрать из них наи- 
лучшее; он может только попытаться обойти препятствие "всле- 
пую". Для решения этой проблемы предлагались эвристические 
алгоритмы, не гарантирующие достижения цели. Созданы также ал- 
горитмы обхода препятствий "вслепую", гарантирующие достиже- 
ние цели. Эти алгоритмы могут быть отнесены к двум различным 
классам: алгоритмов, не создающих модели среды в памяти ро- 
бота, и алгоритмов, создающих и использующих проблемно-ориен- 
тированную модель. 

Алгоритм, относящийся к первой группе, предложен в [I] . 
Рассматривается круглый робот, перемещающийся по поверхности. 
Предлагается планировать перемещение некоторой вспомогатель- 
ной точки (ВТ) , не связанной жестко с роботом (например, точ- 
ки контакта тактильного датчика о препятствием) и задавать 
способ преобразования маршрута ВТ в маршрут робота. Предлагае- 
мый алгоритм позволяет строить маршрут ВТ на плоскости, содер- 
жащей заранее неизвестные зоны, имеющие произвольную форму и 
расположение в рабочей области. В работе приведены доказатель- 
ство сходимости алгоритма и результаты экспериментальной про- 
верки алгоритма. Недостатком работы является то, что алгоритм 
не дает оптимального маршрута. 

В [2] рассмотрен пример алгоритма, относящегося ко второй 
группе. Алгоритм основывается на описании задачи о помощью 
графа. Граф отроится по клеточному представлению чаотично не- 
известной среды с использованием не двухцветной, а трехцвет- 
ной карты местности, где наряду о "белыми" (свободными) и 
"черными" (непроходимыми) клетками входят в рассмотрение "се- 
рые" клетки, о которых робот покы не имеет достоверных сведе- 
ний. По мере попадания "серых" клеток в информативную окрест- 
ность движущегося робота, они раскрашиваются либо в черный, 
либо в белый цвет. Неявное задание трехцветного графа и его 
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последовательное наращивэние по локальной (и даже контактной) 
информации не вызывает принципиальны:: затруднений. Недостат- 
ком этого алгоритма является то, что о увеличением подробнос- 
ти карты местности экспоненциально возрастает вычислительная 
сложность. 

В большинстве случаев движение манипулятора с большим чис- 
лом степеней подвижности можно разделить на два типа: позици- 
онирование (или движение "в большом") и ориентацию охвата. 
Предлагаемые в [3] 'алгоритмы служат для формирования движения 
”в большом". 

В предлагаемых алгоритмах рабочее пространство представ- 
ляет собой плоскость, границы которой определяются кинемати- 
ческой схемой манипулятора. Начальное и целевое положения 
концевой точки манипулятора заданы в рабочем пространстве. Ра- 
бочее пространство может содержать конечное число препятствий, 
ограниченных простыми замкнутыми кривыми конечной длины. Сог- 
ласно используемому подходу вводился в рассмотрение некото- 
рое [-( -пространство, являющееся отображением рабочего про- 
странства (в І-а пространстве манипулятор отображается точ- 
кой). В работе также рассматривается возможность использования 
данного алгоритма для пята типов плоских манипуляторов. 

Один из наиболее общих подходов к планированию маршрута 
движения мобильных роботов основывается на сжатии движущегося 
объекта в точку и увеличении препятствий для компенсации раз- 
меров действительного объекта. Карта безопасных прямолинейных 
траекторий между раскрытыми вершинами препятствий, представлен- 
ных многоугольниками , * -зываѳтся графом видимости. Траектория 
минимальной длины может быть найдена с использованием поиска 
по графу видимости [4] . 

Преимущество этого метода заключается в том, что легче 
рассчитать пересечение точки о набором объектов, чем вычислить 
пересечение объектов. Недостатком этого метода в применении 
к манипуляторам, является то, что размеры перемещаемого объзк- 
та должны быть известны при разработке алгоритма и на могут 
изменяться. Другой недостаток заключается в том, что этот ал- 
горитм дает траекторию, которая стремится быть очень близкой 
к препятствиям. 
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Похожий метод был применен к задаче в трехмерном простран- 
стве с пѳремешатѳльными и вращательными .движениями объектов. 

В этом случае вычислительная эффективность очень низка, не 
существует гарантии, что полученная траектория будЬт оптималь- 
ной [5] . 

В (б] описан алгоритм нахождения пути точечного робота 
с минимальным числом поворотов в свободной области, ограничен- 
ной многоугольником (не обязательно выпуклым) из данной началь- 
ной точки в любую заданную целевую точку. При формализации за- 
дачи предполыэется, что многоугольник "простой", т.ѳ. огра- 
ничен простой замкнутой кривой, которая делит плоокость на 
две области: неограниченную - внешнюю и ограниченную - внут- 
реннюю. Точка плоскости считается точкой многоугольника, если 
она ’нѳ находится во внешней области. 

В [7] предложен алгоритм, основанный на точном описании 
свободного пространства. Свободное пространство делится на 
три части: I) основное рабочее пространство - пространство, 
определенное набором конструктивных ограничений на диапазоны 
перемещения каждого звена; 2) частично свободное - простран- 
ство, где возможно столкновение захватного устройства о пре- 
пятствием; 3) полностью свободное' пространство , в котором 
гарантируется отсутствие столкновений яри всех допустимых 
ориентациях. В работе обосновывается переход от планирования 
перемещений всей руки к планированию перемещений захватного 
устройства. Рабочее пространство представляется в виде кле- 
точной структуры: вое смежные клетки, принадлежащие простран- 
ству одного класса , помечаются одинаковыми метками. 

Другое точное описание свободного пространства предлагает 
в своих работах Брукс Р.А. В [в], [э] вадача решается черев 
описание свободного пространства о помощью обобщенных цилинд- 
ров. Автор пытается найти траекторию так, чтобы она являлась 
центром свободного пространства. Хотя полученная траектория 
может быть не оптимальной, она бѳаопаона. Дополнительным пре- 
имуществом является то, что пространство зедаѳтоя независимо 
от перемещающихся объектов. Найденная траектория является бо- 
лее сложной, так как перемещается объект, а не точка. 

В {ІО] свободное пространство представляется в терминах 
двух элементарных форм: обобщенные конусы и выпуклые много- 
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угольники. Обобщенные конусы хорошо описывают узкие участки 
свободного пространства между двумя препятствиями. С другой 
стороны выпуклые многоугольники удобны как раз яри списании 
больших областей, а в случае узких коридоров описание оказы- 
вается слишком громоздким. 

Предлагаемый подход к планированию маршрута состоит в 
том, чтобы попарно проанализировать пространственные соотно- 
шения между соседними препятствиями и построить граф "оосед- 
ства" препятствий. После этого по дуальному (по отношению к 
построенному) графу из обобщенных конусов и выпуклых много- 
угольников можно сформировать соответственно каналы и прохо- 
димые области. 

В [іі] авторы подходили к задаче с помощью представления 
отношений между движущимся объектом ь препятствиями в виде 
вращательного планирующего графа ( гвШт тарріпд дгар/і). 

Задача нахождения траекторий, свободных от столкновений, таким 
образом, трансформировалась в рассмотрение связности графа. 

В литературе встречаются еще некоторые способы представ- 
ления свободного пространства. Н тример, в [12] свободное про- 
странство представляется следующим образом. Через стороны 
препятствий, представленных в виде многоугольников, проводят- 
ся прямые линии. Всевозможные пересечения этих прямых линий 
образуют множество соприкасающихся многоугольников, которым 
в соответствие ставятся узлы графа свободного пространства. 

В другой работе препятствия представляются прямоугольниками 
оо сторонами, параллельными осям координат [13] . Прямые, про* 
веденные через стороны препятствий, разбивают рабочее прост- 
ранство на нерегулярную решетку. Затем методом Оиіпе Мс СРу&кау 
из множества свободных прямоугольных областей выделяются мак- 
симальные, пересекающиеся, прямоугольные , выпуклые , свободные 
области, которым в соответствие ставятся узлы графа. Планиро- 
вание траекторий движения робота в этих случаях происходит на 
основе поиска на графе свободного пространства. 

Среди моделей свободного пространства, большинство из ко- 
торых базируется на формальном описании задачи с помощью тех 
или иных графовых структур, несколько особняком стоит представ- 
ление рабочей среды в виде поля потенциалов [2] . Однако метод 
потенциалов может быть применим к задачам обхода лишь доволь- 
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но простых (желательно выпуклых) и расположенных на достаточ- 
ном расстоянии одно от другого , препятствий. 

В связи с проблемой обхода препятствий роботами, пред- 
ставляет интерес задача перемещения многоугольного, не обя- 
зательно выпуклого, предмета А в стационарной среде с препят- 
ствиями на плоскости, известная как задача о ''перемещении 
кресла через дверь". Предполагается, что граница объекта явля- 
ется простой замкнутой кривой. Одно из оѳшѳний такой задачи 
приведено в статье [14] . 

Как можно было заметить, характеристика любого метода 
планирования траекторий движения робота во многом определяет- 
ся принятым сп-ѵ-'Хюм формального описания задачи , в котором 
можно выделать два уровня: первичное представление карты мест- 
ности и построенная по нему абстрактная модель рабочей среды, 
непосредственно используемая для планирования маршрута робота. 
Основными способами представления карты местности являются: 
представление в виде клеточных, контурных и других моделей; 
во учет реальных габаритов робота требует определенного про- 
межуточного преобразования рабочего пространства для сведе- 
ния задачи к планированию движения фиктивного, точечного робо- 
та с обходом поѳвдопрапятотвий , расширенных, по сравнению с 
истинными препятствиями, на величину некоторой зоны безопас- 
ности. Решение общей задачи перемещения объекта затрудняется 
при необходимости учета ориентации перемещаемого тела. Одна 
из наиболее интересных идей преодоления этой проблемы состоит 
в переходе в конфигурационное пространство. 

В [15] приведен алгоритм планирования траекторий движения 
робота-манипулятора, использующий конфигурационное простран- 
ство. Конфигурация движущегося объекта - это набор параметров, 
полностью определяющих положение любой точка объекте. Например, 
набор значений шарнирных углов манипулятор* определяет его 
конфигурацию. Конфигурационное пространство движущегося объек- 
та это пространство его конфигураций. Следовательно, простран- 
ство шарнирных углов робота есть эхо конфигурационное прост- 
ранство. Напротив, декартовы координаты охвата робота не опре- 
деляют конфигурация робота ввиду неоднозначности решения об- 
ратной задачи кинематики. К недостаткам описанного метода от- 
носится экспоненциальный рост сложности алгоритма от числа 
степеней подвижности. 


Еще один алгоритм, использующий конфигурационное простран- 
ство, предложен в [16] , В большинстве известных работ исполь- 
зуются алгоритмы автоматического планирования траектории, но 
в ряде случаев предпочтительнее работать в интерактивном ре- 
жиме. Интерактивный режим может быть обеспечен путем удобно- 
го представления в графическом виде конфигурационного прост- 
ранства, например, на экране дисплея, что позволяет оператору 
прокладывать траекторию с помощью курсора. Такой алгоритм и 
рассматривается в [16] . 

Из обзора видно, что алгоритм планирования траекторий дви- 
жения робота-манипулятора о п -степенями подвижности облада- 
ет значительной вычислительной сложностью. Но существует воз- 
можность перехода в зависимости от решаемой задачи . к планиро- 
ванию траектории движения только первых трех звеньев или толь- 
ко одного захватного устройства. В последнем случае можно 
испольаовэть алгоритмы, разработанные для мобильных роботов. 

В некоторых случаях возможен переход от планирования в трех- 
мерном пространстве к планированию траекторий движения робота 
в двухмерном пространстве. 
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